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SUMMARY 


A numerical study to determine the temperature distribution 
in the guide vane blades of a radial inflow turbine is pre- 
sented. A computer program has been developed which permits 
the temperature distribution to be calculated when this blade 
is cooled internally using a coit±)ination of impingement and 
film cooling techniques. The study is based on the use of the 
finite difference method in a two dimensional heat conduction 
problem. The results are then compared to determine the best 
cooling configuration for a certain coolant to primary mass 
flow ratio. 
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INTRODUCTION 


The power output of a gas turbine is directly proportional 
to its inlet gas temperature. The inlet' temperature of a gas 
turbine is controlled however by metallurgical limitations. 
Hence, effective cooling of the turbine blades would consi- 
derably improve the power output by allowing higher inlet gas 
temperatures. Thermal stresses can be minimized by reducing 
temperature gradients in the blade. This can be accomplished 
by matching the coolant distribution and internal heat transfer 
coefficients with the external heat flux distribution to achieve 
the desired temperature distribution. The cooling of a turbine 
blade therefore serves two purposes: 

(a) Reducing the temperature level in the blade 
material. 

(b) Equalizing the temperature differences in the 
blade section. 

There have been several experimental and theoretical 
studies dealing with turbine cooling. A good discussion of 
the different cooling techniques, namely convection, trans- 
piration and film cooling is given in Reference [1] . References 
[2, 3, 4, 5]p-resent a review of the present state of the art 
for the internal cooling of turbine nozzles in aircraft appli- 
cations. Okapuu and Calvert [3] have presented design 
features of the turbine and the results of thermal and stress 
analyses of its components. Kuhl [6] reported temperature 
’measurements obtained within the cooled and uncooled turbine 
blades, Hamed, Baskharone ^and Tabakoff [7] investigated the 
temperature distribution in the rotor of a radial inflow 
turbine. Tabakoff and Clevenger [8] carried out an experimental 
investigation of heat transfer characteristics for various 
configurations of air jets impinging on a concave surface 
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representing the leading edge of a cooled blade. Metzger, 
Yamashita and Jenkins [9] have experimentally investigated 
the heat transfer characteristics between single lines of 
circular jets and concave cylindrical surfaces. Kercher and 
Tabakoff [10] studied the heat transfer by a square .array of 
round jets impinging perpendicular to a flat surface. Chupp, 
Helms , McFadden and Brown [11 ] were concerned with impinging 
characteristics of a concave surface sized to match the 
leading edge of a typical turbine blade. 

Experimental development work for blade cooling techniques 
is very expensive. It also does not provide a deeper insight 
into the local physical interrelations ovzing to the complex 
measuring techniques involved. Computational methods can 
therefore be developed to provide valuable information about 
the temperature distributions resulting from the different 
cooling techniques. In the present study a numerical method 
is developed to determine the temperature field in a radial 
turbine guide vane blade. This analysis is based on the finite 
difference method in the two dimensional conduction problem. 

A computer program has been developed for calculating the 
temperature distribution throughout the blade section. Some of 
the available experimental test results are drawn on to back 
up the theoretical procedure. 

The blade profile used in the present investigation is 
shown in Figure 1 . The computations are carried out to deter- 
mine the resulting temperature distribution v>7hen this blade 
is cooled internally using a combination of impingement and 
film cooling techniques. The impingement tube configuration 
and the position of coolant discharge slot were varied to 
study their effects on the blade temperature distribution. 

The resulting cooling effectiveness and temperature distri- 
butions are reported . 
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ANALYSIS 


A radial turbine guide vane blade does not usually 
experience high axial temperature gradients at its inlet. 
The heat conducted to the two guide vane back plates is 
also small. The use of a two dimensional solution to 
determine the blade- temperature distribution is therefore 
justified. 


Governing Equations 

Under the steady state conditions and assuming that the 
conductivity of the material is constant, the differential 
equation governing the temperature distribution in the turbin 
blade is given by; 


2 2 
3 T ^ 3 T 

7 2 

3x 3y 
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( 1 ) 


There are three different types of boundary conditions 
which may be assigned to the different blade surfaces. They 
are ; 

i) The boundary conditions on the gas side surface 
at points far from being affected by coolant 
discharge . 

ii) The boundary conditions on the gas side surface 
at points downstream of the slot. 

iii) The boundary conditions on the coolant side 
surface. 

The general convective boundary condition is written as 
follows : 
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is the thermal conductivity of the blade material. 

is the temperature gradient normal to the surface. 

represents the local convective heat transfer 
coefficient evaluated at the reference stream 
temperature, T^. . 

is the local temperature on the blade surface. 

The value of the reference stream temperature T^, depends upon 

the type of the boundary point. In the case of a point on 

the hot gas surface which is not affected by the coolant 

discharge, the reference stream temperature takes the value 

of the free stream gas temperature, T . If the point happens 

to be on the coolant surface, then assumes the value of 

the temperature of the coolant, T . On the other hand, if the 

o 

point is downstream of the point of coolant discharge, the 
value assigned to T is the adiabatic wall temperature, T , 

0 0W 

Finite Difference Equations 

In this section, the differential equation (1) and 
the boundary condition, equation (2) , are recast into a 
finite difference form which will be solved using an iterative 
successive over relaxation technique. Referring to Figure 2, 
for the mesh point P(i,j), its neighboring points A and C 
are situated at and fractions of DX and points B and 
D are situated at 5^ and §2 fractions of DY. The following 
finite-difference equation is obtained by expressing the 
temperatures at A, B, C and D in terms of the temperature at P 

■ -2-^ - T 

DX^ (5^+52) ^2 ^2 

+ — 2 " T. . (^ + ^)+ = .0 

DY ^( 6 ^+ 62 ) ^1 ^2 *^2 

(3) 
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The above equation can be written in the follov/ing more con- 
venient form: 



T. , . 


•j+i,i 




2 T 

r iri-1 



(4) 


(5) 


For equally spaced grid lines the factors 5^/ ^ 2 ' ^2 
be equal to unity and the equations (4) and (5) will reduce 
to the following form: 


2" ^^^^^^i+1 i Vl DX^(T. 

2 (DX^ + DY*^) i+J-O 1 1/j ifD+1 


Referring to Figure 3, the convective boundary condition, 
represented by equation (2) , is used to evaluate the temperature 
at point D on the blade surface. If NDN' is the normal at 
the boundary point D in Figure 3, then the finite-difference 
approximation of the normal temperature derivative at D is 
given by : 


T - T 

9T _ D N 
8n DN 


where is the temperature at point D 

T-, is the temperature at point N 

DN is the distance between points D and N„ 

Equation (2) can therefore be written in the following finite- 
difference form; 
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COOLING CONFIGUPATION 

The purpose of the present investigation is to predict the 
temperature distribution in a radial turbine guide vane blade, 
using different cooling techniques. The present chapter is 
intended to give the description of the various cooling con- 
figurations considered in the present study and their basic 
terminology. 

The nozzle blade, due to its small size, presents a formidable 
design challenge. The internal cooling passage configuration 
producing a satisfactory vane temperature distribution is one of 
considerable complexity. Various investigators have pointed 
out that transpiration cooling is the most effective means 
of cooling turbine blades. Unfortunately, such cooling system 
is very complicated from the mechanical point of view. One way 
to avoid these difficulties and in the meantime retain a high 
cooling efficiency, is offered by the simultaneous application 
of impingement and film cooling {Figure 4a) . This is accomplished 
by placing a tube inside the hollow airfoil and orienting the 
tube such* that a series of holes machined into the tube are 
opposite the inner surface of the blade wall. Cooling air intro- 
duced into the tube (Figure 4b) exits through the holes as a 
series of jets directed towards the inner surface. The cooling 
air can later be ejected through a slot near the blade trailing 
edge, to provide film cooling. 

Figures 5A and 5B show schematically ' the four different 
impinging jet geometries which were investigated in the present 
study. A schematic representation of the three different slot 
geometries considered for coolant discharge is shown in Figure 6. 
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The impingement tubes are identified by a number Cl? 2, 3 or 4) 
and the slot positions are identified by a letter (A, B or D) . 
In a particular cooling configuration one of the coolant 
discharge slot geometries is used in conjunction with one of 
the impingement tubes. A particular cooling configuration will 
therefore be identified by a combination of a letter and a 

I 

number. Table 1 gives a description of the seven cooling 
configurations investigated. In all cases the mass flow rate 
of coolant constituted 3% of the main flow rate. 

Once the blade geometry and its cooling configuration are 
defined, the external and internal flow parameters can be 
determined according to the procedure described in Chapter 4. 


TABLE 1 


Slot 

Position 

A 


Impingement 
Tube No . 


Cooling 

System 

A1 


A2 


A 


A3 


B 


B2 


B 


B3 


B' 


B4 


C3 


Description of the Cooling System 


The coolant is injected through a 
spanwise slot in the insert and is 
discharged on the lower surface. 

The coolant is injected through a 
spanwise slot in the insert and is 
discharged on the lower surface. 

The coolant is injected through a 
spanwise slot and an array of holes 
in the insert and is discharged 
on lower surface. 

The coolant is injected through a 
spanwise slot in the insert and is 
discharged on the lower surface near 
trailing edge. 

The coolant is injected through a 
spanwise slot and an array of holes 
in the insert and is discharged on 
the lower surface near trailing edge 

The coolant is injected through a 
spanwise slot and an array of holes 
(opposite suction side) in the 
insert and is discharged on the 
lower surface near trailing edge. 

The coolant is injected through a 
spanwise slot and an array of holes 
in the insert and is discharged on 
the upper surface near trailing edge 
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PRELIMINARY CALCULATIONS 


In addition to the blade geometry, the values of the flow 
velocity, reference stream temperature, heat transfer coefficient 
over the entire blade surface, and the coolant mass flov; dis- 
tribution are required in the blade temperature computation. 

In the following the procedure for determining these parameters 
will be described. 

For an essentially subsonic flow, the nozzle vane velocity 
distribution can be obtained using the computer program of 
Katsanis [12], Knowing the velocity, the gas temperature can 
then be calculated using the energy equation. In the present 
analysis, the nozzle vane geometry and the main flow conditions 
are taken similar to the experimental turbine of Reference 
[5] . The flow channel between the two neighboring nozzle 
vanes is shown in Figure 7. The flow accelerates on both the 
pressure and suction surfaces, up to the throat. The dotted 
line represents the sonic line for the operating conditions 
under investigation. The flow experiences further acceleration 
on the blade suction side up to the trailing edge. The gas 
velocity, the gas temperature and the heat transfer coefficient 
distributions over the nozzle vane surface under investigation 
are shown in Figures 8, 9 and 10 respectively. This data was 
available in the experimental turbine study of Reference [5 ] . 

The heat transfer coefficients of Figure 10 v?ere obtained 
without taking into consideration the effect of coolant discharge. 
It has been shown in References [14 ] and [15 ] that the heat 
transfer coefficient does not change appreciably in the presence 
of film cooling, for low coolant mass flow rates. Figure 11 
presents the ratio of the heat transfer coefficient h, obtained 
with film cooling to the heat transfer coefficient h^, without 
film cooling plotted over the dimensionless distance d/w 
downstream from the slot. This figure was plotted in Reference 
[14] using the data of Reference [15] and is reproduced here. 

The parameter on the curves is the ratio of the density p times 
the exit velocity u^ from the slot to the external mass velocity 
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p u in the mainstream outside the boundary layer. It can be 

observed that starting with a distance d/w = 22, the heat 

transfer coefficient differs by less than 10 percent as long 

as the ratio p u /p u is smaller than 1. In the present 

investigation the mass flow rate of coolant constituted 3% of 

the main flow rate and the ratio p u /p u was less than 0.-4. 

c e g g 

This justifies the use of the same values of the convective 
heat transfer coefficients of Figure 10, downstream of the 
discharge slot. At these locations, the adiabatic wall 
temperature, was determined using the empirical relations 

of Reference [16] . 


Tip 




( 10 ) 


1,9 


Pr2/3 


^ Cp y 0.2 
0 329 Re — 2- C— ) 

^c c 


pud 0.8 
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m 



u 
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-5 ^c 

S = 1 + 5 X 10 Re sina 

s y 

g 

is the film cooling effectiveness, 
is the adiabatic wall temperature, 
is the temperature of the hot gas, 

is the temperature on the blade surface, 
is the specific heat at constant pressure, 
is the coolant mass flow per unit length of the slot, 
is the Prandtl number, 

= m /y is the slot Reynolds number, 
c c 

is the velocity, 

is the distance of the point from the point of 
coolant discharge, 
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a is the ejection angle, 

u is the viscosity, 

p is the density, 

g subscript for the gas side parameter, 

c subscript for the coolant' side parameter- 

The physical properties of the coolant such as p, Cp, y depend 
on the film coolant temperature T^. Initially this temperature 
is assiomed to be constant and equal to the temperature of the 
coolant in the impingement tube, T^. Subsequently it is 
corrected in the computer program and the temperature of the 
coolant which increases from leading edge to trailing edge is 
also obtained using the heat balance equations. 

The convective heat transfer coefficient over the internal 
blade surface, depends on the geometry of the impingement 
tube as can be seen from Figures 12 through 15. Impingement 
tubes niimber one and two provide cooling only near the leading 
edge with the rest of the internal passage being mainly cooled 
by convection due to the coolant flow through the passage 
between the impingement tube and the blade inner surface. 
Impingement tubes number three and four with the arrays of 
■romd jets over the tube surface, in addition to the leading 
edge slot provide impingement cooling over a larger portion 
of the blade surface. It therefore results in higher values 
of film coefficient Chi over the intermediate portion of the 
internal cooling passage as shown in Figures 14 and 15, The 
experimental results of Reference [8J were used for the 
calculation of h at the leading edge in all cases. For 
impingement ttibes number three and four, the experimental data 
of Reference [10] for jet arrays impingement in the coolant cross 
flow was used. The forced convective heat transfer coefficients 
due to the coolant flow either betv/een impingement number two 
and the blade inner surface or between the tv;o internal blade 
surfaces downstream of the tubes were calculated using the 
empirical relation for flow between parallel plates. 
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0,014 He 


Nu. = 


0.81 


C12) 


The mean hydraulic diameter V7as used as the characteristic length 
in the calculation of Nusselt and Reynolds numbers. The 
above equation is based on the experimental work at' Rolls Royce 
and gives a value slightly higher (up -to 5%) than the 
traditional equation of Reference [19 ] which is 

Nu = 0.023 Pr'^*'^ Re°’® (13) 


All the preliminary data necessary in the determination 
of the temperature distribution in the blade cross-section 
is now available. In the follov/ing chapter, the computer 
program for blade temperature computations and its input 
variables, will be described. 

PROGRAM DESCRIPTION 

In the following, a description of the computer program, its 
input and output and other information regarding its use will 
be given. The program listing is given in Appendix A. 

Main Program and Subroutines 

There are two separate main programs listed in Appendix A. 

The first program is concerned with the determination of temperature 
distribution throughout the blade section having coolant discharge 
on the blade outer lov/er surface (e.g., configurations Al, A2 , A.3, 

B2 , B3, B4) . The second main program does the same job for the 
blade section having coolant discharge on the blade outer upper 
surface (e.g., configuration C3) . Many of the computational patterns 
are similar in the two main programs and it might seem that the two 
could be efficiently combined. However, the differences are funda- 
mental enough so that the cause of clarity is best served by keeping 
them separate. One of the two main programs is used along with all 
the stibroutines described later in this chapter, to obtain tem- 
perature distribution throughout the blade section. 

The main program begins by calling subroutine DATA. This 
instruction reads in all of the input data required. The 
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appropriate formats and description of the input are discussed 
in sections (5.2) and (5.3). As an initial guess the temperatures 
at all grid points are set to any convenient value. The program 
then calls subroutine TEMAD which calculates the adiabatic wall 
temperatures at points downstream from the slot using equation 
(11) , At points downstream from the slot, the gas temperatures 
will be replaced by the corresponding adiabatic wall temperatures 
in the main program. 

The Gauss-Seidel iteration procedure is used to solve 

equations (4) and (9) , giving the better values of the temperatures 

at the grid points. The convective boundary condition (equation 9) 

is used to evaluate the temperature on the blade surface point 

such as D in Figure 3. Referring to this figure, it is necessary 

to determine the normal DN at point D, in order to use equation 

(9) . Subroutine CUR finds coordinates of points F and C which 

are adjacent to point D. Subroutine SLO finds an equation of 

2 

the curve y = a + bx + cx passing through points F, D and C. 

The constants b and c in the parabolic equation are returned to 
the main program to get the slope of the tangent to the surface 
and hence slope of the normal to the curve at point D. If 
■the adjacent grid points on the blade surface happen to be 
very close to each other, the value of the slope of the tangent 
obtained by subroutine SLO sometimes is not correct. In such 
cases, the correct value of the slope must be directly specified 
by inserting a card in the main program. The point N where the 
normal DN intersects one of the closest mesh lines inside the 
blade surface can thus be found. The temperature T^^ at this 
point can be determined by linear interpolation using the 
temperatures at points A and P. Subroutine MESH calculates 
the temperature T(i,j) at any interior mesh point P (see 
Figure 2) . It determines $2 ^2 Figure 2 

and knowing the temperatures at neighboring grid points 
proceeds to find the new temperature T(i,j) using equation (4). 

The temperature of the cooling stream increases along its 
flow path due to the heat exchange with the hot blade. The new 
coolant flow temperature computations are performed following 
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the computations of new blade temperatures for each iteration. 

The heat balance equations which are used to determine the 
coolant temperature rise over an incremental length of the 
passage are derived in Appendix B. The local blade temperature 
and the local coolant temperature are interdependent. For 
initial set of iterations required to achieve the desired 
accuracy, the temperature of the film coolant T^ is assTomed 
to be the same as the temperature of the coolant T^ in the 
impingement tube. The actual film coolant temperature T^ which 
is higher than T^ and which is determined using heat balance 
equations, and the physical properties of the coolant at this 
temperature are used in subroutine TEMAD to obtain new 
adiabatic wall temperatures which are closer to actual adiabatic 
wall temperatures. The new values of the film coolant temperature 
and adiabatic wall temperatures are then used in the boundary 
conditions of the next iteration. This process is repeated 
until the desired accuracy is achieved. 

The accuracy of the convergence process is specified by a 
number e. The iteration process is continued until the sum 
of the squares of differences in temperatures between two 
successive iterations is less than e. Mathematically v/e can 
write, 

(T - T , )^ < e (14) 

m ^ 1 , H n vL 

all poa_nts 

til 

where T^ is the temperature after n iteration, 

til * 

T^_^ is the temperature after (.n-1) iteration,- 
and e is the prescribed error limit - 

One of the problems of Gauss-Seidel method of solving 
equations is that it is relatively slov; to converge to the 
solution- The method of successive over-relaxation has been 

used to accelerate the convergence process. In this technique, 

« th 

T^ the temperature to be used for (n+1) iteration is obtained 

by using the following relation: 
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CO 


(15) 


n 


+ (I--) Vl 


where '^n-1 temperature after (n-1) iteration, 

is the temperature after n iteration, 
and 0 ) is the overrelaxation factor. , 

I 

For 0 ) = 0 the new value of the temperature T^, would be 

identical to the old value ^^-1' hence no progress would 

be made. For a = 1 the new value of the temperature T , 

would be the same as T calculated in the Gauss-Seidel 

n 

Procedure. Values of e greater than unity would represent an 

overrelaxation or an extrapolation beyond the Gauss-Seidel 

value for temperature. It is this extrapolation process that 

can accelerate the convergence of the iteration process. 

Since the iteration process will not converge for values of 

tu greater than 2, the value of o will be between 1 and 2. 

After the convergence of 'the solution is achieved, 

temperatures at all the grid points are printed. The program 

then calls subroutine TEMP which determ,ines the maximum and 

rainim\am temperatures in the blade along any I-line and finds 

the locations of all the temperatures lying betv;een these two 

temperatures. Subroutine ISOTH determines isothermal line 

locations in the. blade. The desired number of isotherms and 

their temperatures are specified by a DATA statement in a 

BLOCK DATA subprogram. Finally, the cooling effectiveness n. 

c 

is evaluated. This parameter facilitates the comparison of 
different cooling configurations and is defined as follows: 


T - T 

„ ^ _g s 

^c T - T 

g c 

where ri is the cooling effectiveness, 
o 

Tg is the local temperature of the hot gas, 

Tg is the local blade temperature obtained using 
particular type of cooling, 

and T is the initial temperature of the coolant. 


(16 ) 
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In the following sections the description of the input variables 
as well as the instructions for preparing the input data are 
given. 

Input Variables 

There are ten sets of input data. They are fed in the 
program in the following order. The variables appearing in 
these sets will be described in section (5.3). 

1. Physical parameters of the grid system: 

DX, DY, NXO, NXI 

NX, NY., Nl, N2, N3, N4 , N5, N6, N7 , NP 

2. Properties of the blade: 

LB , XK, CP (3P8.6) 

3. Properties of the coolant: 

PPNO, ALPHA, TCOLN, TEMPO 
DENSC, CPC, VI see, FLOC 

4. Iteration checks for required convergence: 

OME, SUMM, NTE, IMAX (2F8.3, 2 15) 

5. Blade outer and inner boundary point coordinates and their 
corresponding sequence numbers : 


(Y1(I) , 

1=1/ 

NXI) 

(12F6.3) 

(IP(I) , 

1=1, 

NXI) 

(2014) 

(Y2(I) , 

HI 

II 

I-* 

> 

NXI) 

{.12P6.3) 

(IP (I) / 

1=1, 

NXI) 

(2014) 

(XI (J) , 

J=l, 

NYl) 

(12F6.3) 

(IP(J) , 

J=l, 

.NYl) 

(2014) 

(X2 ( J) , 

J=l, 

NYl) 

(12F6.3) 

(IP(J) , 

J=l, 

NYl) 

(2014) 

(Y3(I) , 

I=N11, N21) 

(12F6.3) 

(IP (I) , 

I=N11, N21) 

(2014) 

(Y4 (I) , 

I=N11, N21) 

(12F6.3) 


C4F8.2) 

C2P13.10, E16.10, F8.4) 


(2F5.3, 215) 
(1015) 
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(IP CD, I=N11, W21) 

(2014) 

(X3(J), J=N31, N41) 

C12F6 .3) 

(IP(J), J=N31, N41) 

(2014) 

(X4(J), J=N31, N41) 

(12F6.3) 

CIP(J), J=N31, N41}' 

(2014) 

Gas temperatures and initial 

coolant temperatures : 

(TG(I,1), 1=1, NXl) 

(10F8.2) 

{TG(I,2) , 1=1, NXl) 

(10F8 .2) 

(TG(I,3), I=N11, N21) 

(10F8.2) 

(TG(I,4), I=N11, N21) 

(10F8.2) 

CTGX(1,J), J=2, NYl) 

(10F8.2) 

(TGX(2,J) , J=2, NYl) 

(10F8.2) 

(TGX(3,J), J=N31, N41) 

(10F8.2) 

CTGX(4,J), J=N31, N41) 

(10F8.2) 

Heat transfer coefficients on 
boundary surfaces : , 

gas side and coolant 

(H(I,1) , 1=1, NXl) 

(10F8.4) 

(H(I,2), 1=1, NXl) 

(10F8.4) 

(H(I,3) , I=N11, N21) 

(10F8.4) 

(H(I,4), I=N11, N21) 

(10F8.4) 

(HX(1,J), J=2, NYl) 

(10F8.4) 

(HX(2,J), J=2, NYl) 

(10F8.4) 

(HX(3,J), J=N31, N41) 

(10F8.4) 

(HX(4,J), J=N31, N41) 

(10F8.4) 


Velocities at points where adiabatic wall temperatures are 
to be determined : 


{VELG(l,a), I=N51, NXl) 
(VELG(b,J), 1 = 2 , N6) 


(10F8.2) 
(10F8. 2) 



9. Coolant mass flow: 

(MASC(I,3), I=N11, N51) 
{MASC(I,4), I=N11, N51) 
(DELM(I,3), I=N11, N5) 
(DELM(I,4), I=N11, N5) 

10. Surface distance between 
side boundary surface: 

(DELS Cl/ 3), I=N11, N5) 
(DELS(I,4), I=N11, N5) 


tlOF8.4) 

(10P8.4) 

(10F8.4) 

(10F8.4) 

coordinate points on coolant 

(10F8.4) 

(10F8.4) 


Explanation of the Input Variables 

There are two types of input variables, geometric and 
nongeometric. The geometric variables are shown in Figure 16. 
All input variables and other key variables used in the 
program are explained in this section. 


DX 

DY 

NXO 

WXI 

NX 


NY 


I 

J 

Nl, N2 


Grid spacing in the x direction (inches) . 

Grid spacing in the y direction (inches) . 

Total number of mesh line intersections with the 
outer boundary. 

Total number of mesh line intersections with 
the inner boundary, 
s 

Total ntimber of mesh lines in the x direction, 
if NX is greater than 30, the DIMENSION statements 
in the program need to be changed. 

Total number of mesh lines in the y direction, if 
NY is greater than 65, the DIMENSION statements 
in the program need to be changed. 

Mesh line perpendicular to the x axis. 

Mesh line perpendicular to the y axis. 

I mesh lines enclosing the blade inner boundary 
v/ithout intersecting it. 
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N3, N4 
N5, N2 
N3, N6 
Nl, N7 
NP 


LB 

XK 

CP 

PRNO 

ALPHA 

TCOLN 

TEMPC 
DEN SC 
CPC 

VISCC 

FLOC 

OME 

SUMM 


J mesh lines enclosing the blade inner boundary 
.without intersecting it. 

I mesh lines enclosing the coolant discharge slot 
on the blade outer surface, without intersecting it 

J mesh lines enclosing the coolant discharge slot 
on the blade outer surface, without intersecting it 

I mesh lines enclosing the impingement tube- 
without intersecting it. 

Variable to identify the coolant discharge surface. 
In case of coolant discharge on outer lower 
surface, put NP = 0 and for coolant discharge 
on outer upper surface, put NP = 1. 

Blade height (inches) . 

Thermal conductivity of the blade material 
(Btu/hr in °R) , 

Specific heat of coolant at constant pressure 
(Btu/lbm °R) . 

Prandtl number for coolant. 

Coolant discharge angle (degrees) . 

Temperature of the coolant in the impingement 
tube ( °R) . 

Temperature of the coolant for film cooling (®R) . 

3 

Density of the film coolant (Ibm/ft ) . 

Specific heat of film coolant at constant 
pressure (Btu/lbm °R) . 

Viscosity of the film coolant (Ibm/ft-sec) . 

Mass flow of the coolant through the slot per 
unit length (Ibm/ft-sec) . 

Over-relaxation factor. 

Maximum allowable error (e in equation 14) . 
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NTE 


Number of iterations after which a print out 
of the grid temperatures is required. 

IMAX Maximum number of iterations. After IMAX 

iterations the program terminates even if the 
required convergence is not achieved. 

Y1(I) Y coordinates of points on the outer lower surface 

intersected by I mesh lines (inches ) , (see 
Figure 16) . 

IP (I) Sequence numbers corresponding to the above 

points (see Figure 16) . 

NXl Equals to (NX ~ 1) . 

Y2 (I) , IP (I) Y coordinates of points on the outer upper surface 

intersected by I mesh lines Cinches) , and 
corresponding sequence numbers Csee Figure 16) . 

XI (J), IP(J) X coordinates of points corresponding to the 

first intersection of J mesh lines with the outer 
surface, i.e., near outer surface (inches), and 
corresponding sequence numbers (see Figure 16) . 

NYl Equals to (NY - 1) . 

X2(J), IP(J) X coordinates of points corresponding to the 

second intersection of J lines with the outer 
surface, i.e., farther outer surface (inches) and 
corresponding sequence numbers (see Figure 16) . 

Y3(I), IP (I) Y coordinates of points on the inner lower surface 

intersected by I mesh lines (inches) and corresponding 
sequence numbers. ' 

Nil Equals to (N1 + 1) . 

N21 Equals to (N2 - 1) . 

Y4(I), IP (I) Y coordinates of points on the inner upper surface 

intersected by I mesh lines (inches) and corresponding 
sequence numbers . 
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X3(J), IP(J) 


N31 

N41 

X4(J), IP{J) 


TGtl,a) 
TGX (b , J) 
H(l,a) 

HX (b, J) 


VELG(I,a) 

VELG(b,J) 

N51 

MASC(I,a) 


DELM(I,a) 


DELS (I, a) 


X coordinates of points corresponding to the first 
intersection of J mesh lines with the inner sur- 
face (inches), and corresponding sequence numbers. 

Equals to (N3 + 1) . 

Equals to (N4 - 1) . 

X coordinates of points corresponding to the 
second intersection of J mesh lines with the 
inner surface (inches ) , and corresponding 
sequence numbers. 

Gas temperature at I mesh line intersection point 
on the blade surface Ca = 1,2, 3, 4), C°R) . 

Gas temperature at J mesh line intersection point 
on the blade surface (b = 1,2, 3, 4), C®R) • 

Heat transfer coefficient at mesh line inter- 

2 

section point on the blade surface CBtu/hr in ’’R) . 

Heat transfer coefficient at J mesh line inter- 
section point on blade surface (Btu/hr in '’R) . 

Velocity at I mesh line intersection point on the 
blade surface tft/sec) . 

Velocity at J mesh line intersection point on the 
blade surface (ft/sec) . 

Equals to (N5 + 1) . 

(Total mass of coolant entering the fluid element 
at I mesh line, per unit time (see Figures 17 
and 18) , (Ib/hr) . For further description refer 
to Appendix B . 

Coolant mass added to the fluid element between 
I and I+l mesh lines from the impingement tube, 
per unit time (see Figures 17 and 18) , (Ib/hr) . 

Surface distance between -consecutive points on 
inner surface intersected by I mesh lines (see 
Figures 17 and 18) , (inches) . 
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a=l 

a=2 

a=3 

a=4 

b=l 

b=2 

b=3 

b=4 

TCI,J) 

TB (I, a) 

TBX(b, J) 

NSCLT 

TSCLT 


NNY 


Corresponds to outer lower surface-. 

Corresponds to outer upper surface. 

Corresponds to inner lower surface. 

Corresponds to inner upper surface. 

Corresponds to near outer surface. 

Corresponds to farther outer surface. 

Corresponds to near inner surface. 

Corresponds to farther inner surface. 

Temperature at a grid point/ C°R) . 

Blade temperature at I mesh line intersection 
point on the blade surface, (®R) . 

Blade temperature at J mesh line intersection 
point on the blade surface, { °R) . 

Number of isotherms for which locations are to 
be determined. 

Temperatures of the isotherms (NSCLT temperatures 
have been specified with the DATA statement in 
the BLOCK DATA subprogram. These can be changed 
if isothermal line locations for different set 
of temperatures are desired) . 

Dummy variable which counts the number of points 
determined on an isotherm using ISOTH subroutine. 
Initial values of this variable have been set 
to zero in the DATA statement. 


Preparation of the Input Data 

The input to the computer program to determine the temperature 
distribution in the guide vane blade consists of the sets of 
variables indicated in Sections (5.2) and (5.3). In the following, 
more information regarding the preparation of the input data 
will be given. 
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A rectangular grid network is used for the solution of the 
basic differential equation with specified boundary conditions. 

A typical mesh pattern is shown in Figure 16. The cross 
section of the blade is divided into grids of equal increments 
in both X and y directions. The mesh spacing in the x direction 
need not be the same as the. one in the y direction. The mesh 
lines perpendicular to the x direction are referred to as the 
I lines while the mesh lines perpendicular to the y direction 
are referred to as the J lines. The grid network should be 
drawn such that there are not multiple intersection points, 
i.e., each grid line intersects the outer boundary at most twice 
and the inner boundary also at most twice. It is convenient to 
draw the grid network with 1=1 line tangential to the blade. 

In order to read the coordinates of points on the blade 
boundary, each point is associated with a number for purpose 
of its identification. The numbering can start from any nodal 
point on the blade boundary but it is done in the following 
manner for convenience. For the blade outer surface the numbering 
starts from the node point where 1=1 line meets the outer 
blade surface and for the blade inner surface it starts from the 
node point near the leading edge of the blade inner surface 
(see Figure 16) . The (1=1, J=l) point is treated as origin and 
the coordinates of points on the boundary are given with respect 
to this point. The boundary coordinates and their corresponding 
sequence numbers are now defined. 

The next two sets of input data consist of gas temperatures, 
convective heat transfer coefficients and velocities which are 
calculated in Chapter 4. Note that the initial temperature of 
the coolant is assumed to be constant and equal to that in the 
impingement tube, T^. 

Finally, the coolant mass flow distribution may be obtained 
as follows. The total coolant mass flow entering the impingement 
tube per unit time and the impingement tube configuration are 
selected. As mentioned in Chapter 3, the impingement tube may 
have slots or rows of holes for impinging the coolant against 
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the blade inner surface. Knowing the amount of coolant injected 
through each slot or hole, it is possible to calculate the 
mass of the coolant crossing any I mesh line. Figure 17 shows 
typical coolant flow path. The typical control voltimes used 
to calculate the increment in coolant temperature are shov/n 
in Figures 18 and .19 (refer to Appendix B for further description 
of the control volumes) . It can be noted that 

m{i) + Am(i) = m(i+l) (17) 

where m(i) = Coolant mass entering the control volume 

per unit time at I mesh line, 

AmCi) = Coolant mass added between I and I+l mesh 
lines from the impingement tubs per unit 
time. 

m(i+l) = Coolant mass leaving the control volume 
at I+l mesh line per unit time. 

It is very unusual to prepare the input without any errors, 
the first time this program is run. Therefore, it is recommended 
that the first attempt for running the program should allow 
maxim\am of ten iterations , The resulting output should be 
checked carefully. Of particular interest are the values of 
the slopes of tangent to the blade boundary at all the boundary 
points. Any error in the blade geometry input will usually 
result in wild values of some of these derivatives. All other 
preliminary output should be checked to see that it is reasonable. 


Output 

The first part of the output consists of a printout of the 
list of input parameters for the purpose of reference and checking. 
Such data checks are useful in correcting keypunch mistakes 
on the input cards. The program then prints out the boundary 
point number, boundary point coordinates and the coordinates 
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of two adjacent points. The coefficients a, b and c in the 

2 

equation y = a + bx + cx of the parabolic curve passing 
through these three points and the slope of the tangent to the 
curve at the boundary point are also printed. The error e for 
every iteration is then printed. This allows monitoring of 
the convergence process. The grid temperatures as well as 
blade surface temperatures and gas temperatures are printed 
after NTE iterations and also after the convergence of the 
solution is achieved. In the output, the gas temperatures at 
points downstream from the slot are replaced by the adiabatic 
wall temperatures. The next part of the output is a printout 
of the locations along an I line, of blade temperatures which 
have values between the maximim and minimum temperatures along 
that I line. The isothermal line locations corresponding to 
the temperatures specified in the DATA statement of the 
BLOCK DATA subprogram are then printed. Finally, the values 
of cooling effectiveness at all surface points are printed. 

In Appendix A, a sample of the program output is presented. 
This output corresponds to the B3 cooling configuration of 
Figure 29. 


6. RESULTS AND DISCUSSION? 

A computer program has been developed to predict the 
temperature distribution of the cooled guide vane blades. The 
results of the investigation are presented as plots of isothermal 
lines in the blade cross-section, and the cooling effectivness 
variation over the blade surfaces for the various cooling 
arrangements of Table 1. The geometry of the turbine nozzle 
blades used in this study is that of Reference [5 ] and under 
the same main flow operating conditions. The hot gas inlet 
stagnation temperature is 2760 and the gas turbine flow rate 
is 4.9 Ib/sec. The integral nozzle casting consists of 20 
hollow vanes with inserts as shown in Figure 4A. The coolant 
mass flow rate in the nozzle vane v;as taken as 3% of the main 
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flow rate per channel for all the cooling configurations that 
are investigated. 

The temperature distribution and the cooling effectiveness 
corresponding to the A1 cooling configuration are shown in 
Figures 20 and 21 respectively. In this cooling arrangement 
efficient cooling of the leading edge inner surface is achieved 
by the normal impingement of all the coolant through the leading 
edge slot in the impingement tube number 1 of Figure 5A. The 
coolant is divided into two streams around the tube with a 
lower velocity resulting when they rej'oin at the end of the tube. 

The coolant is finally discharged on the blade pressure 
surface. The blade temperatures are low near the leading edge 
inner surface and increase rapidly on both the pressure and 
suction surfaces away from the leading edge. This is expected 
because for a slot jet, the local heat transfer coefficient 
undergoes a rapid decay after its maximum value in the stagnation 
region. Higher metal temperatures are observed in the central 
portion of the blade section (surface distance AS^/S^ between 
20% and 40%) . This is a result of the reduced coolant velocities 
in that region as explained earlier. The metal temperatures 
are low near the point of coolant discharge and increase at the 
points downstream of it. This is expected since the value of 
the film cooling effectiveness, expressed in equation (11), 

is higher near the point of coolant discharge and decreases 
towards the trailing edge. 

Figure 22 shows the isothermal lines in the blade for the 
A2 cooling configuration. In this case, the vane has a closely 
conforming tube which provides impingement cooling to the leading 
edge and gives good convection cooling also, on account of high 
velocity flow along the side wall passages. The cooling 
effectiveness variation on the blade' surface for this cooling 
configuration is shown in Figure 23. The temperature distribution 
in the blade and the cooling effectiveness variation for the 
A3 cooling configuration are shown in Figures 24 and 25 respectively. 
The impingement tube used in this case has a grid of 0.012 inch 
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diameter holes and a spanwise slot in the leading edge. In such 
cooling arrangement/ higher internal heat transfer coefficients 
are achieved in the leading edge as well as in the central 
portion of the cooling passage. As a consequence, lower blade 
temperatures result in both these regions for this particular 
cooling arrangement. 

For the purpose of comparison, the cooling effectiveness 
variations for Al, A2 and A3 cooling configurations are shown 
together in Figure 26. It is evident that the best heat transfer 
performance is obtained using the round jet array, configuration 
A’3. It is also clear that the temperatures are almost the same 
near the trailing edge region for the three different configura- 
tions. This could be expected as a result of the same coolant 
discharge position in all three cases. 

The temperature distributions obtained so far are for a 
blade having the same coolant passage with different impingement 
tubes. In the following, the temperature distributions in the 
blade for various coolant discharge positions are presented. 
Theoretically, it is feasible to exhaust the cooling air back 
into the gas stream at any point on the vane surface. This is 
possible because the static pressure around the surface ox the 
vane is everywhere lower than the stagnation pressure of the 
coolant, except at leading edge. However, when a tube insert 
is used as part of the cooling system, the coolant stagnation 
pressure is reduced and then it is only possible to discharge 
the cooling air back into the gas stream at points on the vane 
surfaces, which are close to trailing edge. The choice in 
practice is a compromise between the requirements of cooling, 
mechanical strength and good aerodynamic performance. The 
present analysis does not consider the aerodynamic effects of 
cooling air discharge on the turbine performance. The penalties 
in terms of turbine efficiency are minimized however for cooling 
air exhaust on the pressure surface near trailing edge. 

The blade metal temperature distribution obtained for the 
B2 cooling configuration is shown in Figure 27. For this 
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configuration the cooling air path near discharge location was 
chosen keeping in mind that the minimum thickness for casting 
is 0.02 to 0.025 inch. The cooling effectiveness on the blade 
surface corresponding to the above configuration is shown in 
Figure 28. The isothermal lines and the cooling effectiveness 
corresponding to B3 and C3 configurations are shown in Figures 
29 through 32. The effect of changing the coolant discharge 
position is demonstrated in Figure 33 which combines the results 
of A3, B3 and C3 cooling configurations. It can be seen that 
the B3 cooling configuration shov7S a better heat transfer 
performance compared to the .^3 and C3 configurations. 

A closer examination of Figure 30 reveals that, generally 
the cooling effectiveness is higher on the pressure side as 
compared to the suction side when the array configuration 
impingement tube of configuration B3 is used. Another cooling 
configuration B4 was studied in order to improve the suction 
surface cooling. In this arrangement, the tube has a closely 
spaced grid of 0.012 inch diameter holes which provides impinge- 
ment cooling on the suction surface, while the pressure surface 
is provided with convection cooling only. Figures 34 and 35 
show the temperature and cooling effectiveness distributions 
corresponding to this configuration. The cooling effectiveness 
variations for B3 and B4 configurations are plotted in Figure 36 
to provide a means of comparison. It is evident that more 
uniform temperature distribution is obtained by employing B4 
configuration and that this configuration demonstrates best 
heat transfer performance in the leading edge area, central 
portion of the blade as well as in the trailing edge area. 

The computer time naturally depends on the particular blade 
geometry under consideration, as well as on the initial guesses 
of metal and coolant temperatures, but did not exceed 35 
seconds in any of the cases reported here. 
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7. CONCLUSION 


A useful numerical technique has been developed to predict 
the temperature field in the cooled guide vane blade of a radial 
inflow turbine. The computer program which has been developed 
can be used to predict the temperature distribution in the con- 
ventional nozzle guide vane with impingement cooling at the 
leading edge and film cooling near the trailing edge. The 
impingement tube can have any geometry but the coolant can be 
discharged through a single slot. This computer program can 
also be used to study the effects of changes in the coolant to 
primary mass flov7 ratio and initial temperature of the coolant 
on the temperature distribution in the blade. In the present 
investigation, the temperature distribution in the cooled guide 
vane blade v/as obtained using four different impingement 
tubes and three different discharge locations. The results 
show that simple convection cooling configurations are not 
particularly suitable for turbine vanes under the combination 
of high gas temperatures and pressures. Of the various cooling 
configurations investigated the case of impingement cooling of 
the leading edge, cooling of suction side by round jet array 
and cooling discharge near the trailing edge was found to give 
the best temperature distribution. According to the numerical 
results obtained the areas which need special attention if ' 
high local metal temperatures are to be avoided are: 

CD Leading edge where external heat transfer coefficients 
are high. 

(2) The suction surface where external heat transfer 

coefficients are high after boundary layer transition 
from laminar to turbulent. 

C3) Trailing edge where difficulties exist in getting 
internal passages with adequate surface area. 
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LIST OF SYMBOLS 


Symbol 

A,B,C 

S 

d 

DN 


DX 


DY 

h 

I 

J 

k 


LB 

• 

m 

NU 

n 

Pr 

Re 


s 

T 

T(i,J) 


T. 

X 




n 


Description 

Letters to identify coolant discharge position 
(see table 1) . 

Specific heat at constant pressure. 

Distance of a point from the point of coolant 
discharge (see figure 11) . 

Normal distance (see figure- 3) . 

Grid spacing in the x direction (see figure 3) . 
Grid spacing in the v direction (see figure 3) . 
Heat transfer coefficient. 

Mesh line perpendicular to the x axis (see 
figure 16) . 

Mesh line perpendicular to the y axis (see 
figure 16) . 

Thermal conductivity. 

Length of the blade (equation b.l). 

Mass flow rate (see figure 18) . 

Nusselt number. 

Normal direction (see equation 2) . 

Prandtl niimber. 

Reynolds number. 

Distance along blade surface (see figure 8) . 
Temperature . 

Temperature at grid point (i/ J ) i (see figure 2) . 
Film coolant temperature. 

Temperature of the coolant in the impingement tube. 
Temperature after n ' iteration (equation 14). 
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p 

V 

a 

£ 




1,2, 3,4 


Subscripts 

aw 

c 

e 

g 


Description 

Temperature on the blade surface (equation 16) , 
Temperature of the gas (equation b.l). 
Temperature of the blade (equation b.l). 
Velocity of the stream. 

Slot width. 

Over-relaxation factor (equation 15') . 

Coordinate directions (see figure 1) . 

Cooling effectiveness (equation 16 ) . 

Film cooling effectiveness (equation 10 ) . 

Additional mass flow rate from the impingement 
tube (see figure 18) . 

Density. 

Viscosity. 

Coolant discharge angle (equation 11 ) . 

Error limit (equation 14 ) . 

Fractions (see figure 2) . 

Numbers to identify the impingement tube (see 
table 1) . 

Adiabatic wall. 

Coolant. 

Reference stream. 

Hot gas . 
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FIGURE 2, TYPICAL GRID SYSTEM WITH UNEQUAL SPACIN6S. 
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FIGURE 3. GRID POINTS NEAR BLADE SURFACE. 
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Width of Slot Jet = W =0.02 in. 

s 

Length of Slot Jet = = 0.339 in. 

Height of Impingement Tube = H = 0.18 in 
3% Cooling Air Through Leading Edge Slot 


Impingement Tube No. 1 (Slot Jet) 



Width of Slot Jet = 0.02 in. 

Length of Slot Jet = 0.339 in. 

Height of Impingement Tube = 0.5. in. 

3% Cooling Air Through Leading Edge Slot 


Impingement Tube No. 2 (Slot Jet) 


FIGURE 5A. SCHEMATIC OF THE AIR JET IMPINGEMENT CONFIGURATIONS, 
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FIGURE 5 



Width of Slot Jet = 0.01 in. 

Length of Slot Jet = 0.339 in. 

Height of Impingement Tube = 0.5 in. 

Round Jet Diameter = 0.012 in. 

Round Jet Spacing = 0.10 in. 

1% Cooling Air Through Leading Edge Slot 
1% Cooling Air Through Grid on Suction Side 
1% Cooling Air Through Grid on Pressure Sid 


Impingement Tube No. 3 (Slot Jet and Array of Round Jets) 



Width of Slot Jet = 0.01 in. 

Length of Slot Jet = 0.339 in. 

Height of Impingement Tube = 0.5 in. 

Round Jet Diameter = 0.012 in. 

Round Jet Spacing = 0.05 in. 

1.5% Cooling Air Through Leading Edge Slot 
of which 1% flows tov/ards pressure side 

1.5% Cooling Air Through Grid on Suction Si 


Impingement Tube No. 4 CSlot Jet and Array of Round Jets) 


b‘. schematic of the air jet impingement configurations. 
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(A) Coolant Discharge on 
Lower Surface 



(B) Coolant Discharge on (C) Coolanr Discharge 

Lower Surface Near on Upper Surface 

Trailing Edge. 


FIGURE 6, SCHEMATIC OF THE.SLOT GEOMETRIES FOP. COOLANT DISCHARGE, 




FIGURE 7. REGION 0 
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heat transfer coefficient 
with film cooling 

heat transfer coefficient 
without film cooling 
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0.0 0.4 0.8 1.2 1.6 

Surface Distance, (inches) 

FIGURE 12. HEAT TRANSFER COEFFICIENT DISTRIBUTION ON INNER BOUNDARY WITH 
3% COOLANT (CONFIGURATION Al). 




Surface Distance, (inches) 


FIGURE 13, HEAT TRANSFER COEFFICIENT DISTRIBUTION ON INNER BOUNDARY WITH 
ll COOLANT (CONFIGURATIONS A2, B2), 




FIGURE 14, HEAT TRANSFER COEFFICIENT DISTRIBUTION ON INNER BOUNDARY WITH 
3% COOLANT (CONFIGURATIONS A3, B3, C3) , 
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FIGURE 15. HEAT TRANSFER COEFFICIENT DISTRIBUTION ON INNER BOUNDARY WITH 
It COOLANT (CONFIGURATION BA)', 
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FIGURE 16 . BLADE GRID NET WORK. 
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FIGURE 18 , SKETCH ILLUSTRATING THE CONTROL VOLUMES FOR COOLANT MASS BALANCE. 
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FIGURE 19', SKETCH ILLUSTRATING THE CONTROL VOLUMES FOR DERIVING HEAT 
'■ BALANCE EQUATIONS. 
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FIGUE 20, NOZZLE VANE TEMPERATURE DISTRIBUTION WITH 11 COOLANT 
(COfIFIGURATIOM Al). 
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FIGURE 21. EFFECTIVENESS OF COOLING ON BLADE SURFACES WITH 3% COOLANT (CONFIGURATION Al). 
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FIGURE 23, EFFECTIVENESS OF COOLING ON BLADE SURFACES WITH 1>1 COOLANT 
(CONFIGURATION A2) , 
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FIGURE 24. NOZZLE VANE TEMPERATURE DISTRIBUTION WITH 3 % COOLA.NT 
(CONFIGURATION A3) 




FIGURE 25, EFFECTIVENESS OF COOLING ON 
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FIGURE 26. EFFECTIVENESS OF COOLING ON BLADE SURFACES WITH 3% COOLANT FOR 
DIFFERENT IMPINGING CONFIGURATIONS. 





FIGURE 27. NOZZLE VANE TEfIPERATURE DISTRIBUTION NITH 3% COOLANT 
(CONFIGURATION B2). 
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FIGURE 29. NOZZLE VANE TEMPERATURE DISTRIBUTION WITH 35? COOLANT 
(CONFIGURATION B3) . 
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FIGURE 30, EFFECTIVENESS OF COOLING ON BLADE SURFACES WITH 3% COOLANT (CONFIGURATION B3) 



FIGURE 31. NOZZLE VANE TEMPERATURE DISTRIBUTION WITH COOLANT 
(CONFIGURATION C3), 
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FIGURE 32. EFFECTIVENESS OF COOLING ON BLADE SURFACES WITH 3% COOLANT (CONFIGURATION C3). 
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FIGURE 33. EFFECTIVENESS OF' COOLING ON Bi 
DIFFERENT SLOT POSITIONS. 
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FIGURE 34, NOZZLE VANE TEHPERATURE DISTRIBUTION \lUW-3% COOLANT 
(CONFIGURATION B4) . 
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FIGURE 36. EFFECTIVENESS OF COOLING ON BLADE SURFACES WITH 3% COOLANT 
(CONFIGURATIONS A3', B3, B4) , 
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PROGRAM LISTING AND SAMPLE OUTPUT 
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r'i o o o o <'^ o r> o o o o o o o 

rvi KJ NJ 


ORlGINAIi PAGE ® 
OF POOR QUAIxTHI 


FCRTP>Cif\ CCf'FLTEP FPGGRAf' TC CET4TN TEMPERATURE DISTRIBUTION 
IN A SLADE USING ‘ THE FINITE DIFFERENCE VEThCD 


*=!= NAIN FRCGRAM NC.l ** 

FCF • CCCLANT CISCFAFGS CN LC’aER CUTEF SURFACE USE 
MAIN PROGRAM NC.l ALCNG VylTH ALL - THE SUSRCUTINES 


Cl MENS I CN TGTEMOG) ,TGXTEM(65) 

DIMENSION XISC( 30,3C ) »VISC(SC,3G3 ,TSCLT{3G) ,NNY{3Q) 

C IM ENSIGN CELS (20,4 ) ,MASC (3Qf ,C { 30r4) ,CELM(20, A) 

DIMENSION T13Q,e.E) ,TE( 30,4) ,T8X( 4, 6£) ,TG( 30,4) ,TGX( 4,65} , 

1H(3G,4J ,hX (4,65 ) ,V1 (30 ),Y2 (30) ,Y3 (3C ), Y4 (30 3, XI (65) ,X2 (65) , 
2X3(65) ,X4( 65 3 ,BX( 200) ,BY(2CQ) ,B> I ( 200) ,8YI (200) , IP(20C) 

CIMENSICN VELG(30,65) o 

DIMENSION EFF(3C,65) 

REAL LB,MASC 
INTEGER=>4 TSCLT 

CCiMMON T ,TB ,TBX ,TG ,TGX ,H , H X , VI , Y2 , Y3 , Y4 , X 1 , X2 , X 3 , X4 , B X , B Y , 8 X I , 8 Y I 
IrlF 

C CM^'CN CX , C Y , NX , N V , NX C ,N X I , I T EP , M , N 2 , N3 , N4 , NP , N 5 , N6 , N7 
CCMMCN/CAT/ NXl,NYl,Nll,N21,N31,N41,N5l,N52,N61,N62,N71,NSR, 
1CELM,?.'ASC,CELS,XK,CM.E,SUMM , IVAX,NTE,7CCLN,L3,CP 
CCMMCN/TEM/ VELG , PR NC , T EM FC ,CFC ,V ISCC ,ALFhA , FLCC ,CENSC 
CCVMCN/ISG/ XISQ,YTSC,TSCLT,NSCLT,NN\ 

CALL CATA 
45 CCNTINUE 

TEMFCPAPY STGPAGE FCP GAS TEMPS DCWNSTREAM GF SLOT 

DC 250 I=N52,NX1 

TGTEM(I)=TG(I ,l ) 

50 CCNTINUE 

DC 255 J=2,N6 

TGXTEM(J)=TGX( 1 ,J) 

55 CCNTINUE 

CALCULATES ACIABATIC UALL TEMPERATURES CG^NSTREAM GF SLCT 

DC 260 I=N52,NX1 

CALL TEMAC( I ,1,TAC) 

260 TG(I,1)=TAC 

CC 265 J=2,N6 

CALL TEMAD( 1,J,TAD) 

265 TGX(1,J)=TAC 

DC 266 I=N52,N21 

TG( I ,4)=TG( 1,1) 

2ft6 continue 
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o n o o I't r > non 


DC 267 J=N31,N6 

TGX(A, J )=TGX( It J) 
267 CCNTIKUE 
2C0 CCNTiNUE 


STAFTING OF GALSS-SS ICELL NETI-CC CF ITERATION ON I LINES 


SL‘N=0 . 

CC 510 1 = 2, NXl 

IL=IFI X( Y1 ( I ) /DV+Q.0001)+2 

Ih=IF!X(Y2{ I J/DY+C.0G01)+1 

IF (ABSI Y2( n-FLCAT{ TF-L J^CY) .LToO. 00001) IF= IF-1 
CLTEP LCi^cF eCUNCARY 


CAPCS WITH NECA CCVVENT TO SPECIFY SLOPES AT SPECIAL POINTS 

E2=FLCAT( I-1)*CX 
C2=Yi(I ) 

CALL CUR( E2 ,C2 ,e ,C tC) 

XNN=E+2.»C*82 
IFd.EC.NXn XNN=0,01 
Y7 = FLCAT(IL->1)^CY 
X= (C2-YT)=i‘XMN + E2 
IF(I.EG.6) XPN=-1.0 

IFU .GT,N5l.ANC.I ,LT ,N2) GC TC 352 
IF{ I.GE.N2 i GC TC 305 
IF( X,GT.X3( lU ) GC TC 330 
305 Ir(ABS{XMN)«LE. 0.015) GC TC 335 
IF(X.LE.E2JGC TC 320 
XR=E2+CX 

IF{ I .GE .N2) GC TC 310 

IF nX3(IL)-XR) .GT.C.OOOl) GC TC 310 

TTT = T{ I, IL ) + {T8X(2,IL)-TU, ID )=^(X-62)/( X3nU-B2) 

GC TC 325 

310 IF( (X2(lLi-XR).GT.G.C0Cl) GC TG 315 

TTT = T( !,IL )+(TBXl2 , IL )-TU,lL ) }^( X-82 )/{X2( IU-B2) 

GC TC 325 
315 CCNTINUE 

IF{ {X-XP ) .GT .0.0001 ) GC TC 33? 

TTT = T( ItIL)+(T(I + ltIL)-Tn,IL))^(X-E2)/CX 
GC TC 325 

320 TTT=T{I, IL )-{Tn,U)-TU-l,IL))=^(e2-X)/CX 
32 5 XNL=( X-E2 )=* IX~B2} + {C2-Y7 )«{C2-YT ) 

XNL = SCRT(XNL ) 

GC TG 34Q 

330 TTT=iaX(3,IL) 

X.NL = {X3(IL) -B2) =S'tX3( IL)-E2)+(C2-YT )^(C2-YT) 

XNL=SCRT(XNL) 
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iLl lU 10 


GC TC 34C 
335 TTT=T(ItIL) 

>NL=YT-C2 
GC 7C 340 
331 TTT=T(I+1,IL) 

XM = CX^CX+ (C2-YT)=^(C2-YT) 

XM = SGPT(XM) 

340 ^L=h(I,l)/XK 

TN5U=(XM*/iL*TG ( 1,1 ) + TTT )/ { 1,4XNL4;SL ) 
350 CCNTIKUE 

{l.-C/^E)=?TE ( I ,1 )4C/^E=^TN£W 
SL.'' = SL:«+{TE( I,1)TT^EV^)4(TB(I,U-T^EV.) 
reel, 1 )=TN£W 
352 CCKTINLE 

I F( I .GT.M.AND. I .LT,^2) GC TC 260 


INTERIOR POINTS STARTING FRC.v LQ^ER END FOR I LINES WHICH CC MGT CUT 
THE INNER BOL.NCAPY 

CC 355 J=IL,IH 
CALL NEShd , J,TNEW) 

TNEW=(1,-CHE)=J=T( I , J )4CNE4TNEW 

SLjM = SUN+{TNEW-T{ I , J) )4(TNEW-T(I ,J) } 

55 TU,J) = TNEW 
GC TC 480 


FCR I LINES INTERSECTING INNER BCUNCARY-ALL INTERIOR POINTS FROM 
LCWER SURFACE CF CLTER BCUNCAPV TC LCWER SURFACE CF INNER 8CUNCARY 

6C- CCNTINLE 

ILT=IFIX(V3(n/CV + C.0001 ) + X 

IP(ABS(Y3{ I )-FLCAT{ ILI-134CY) .LT.C.OQOOl) ILI=ILI-1 

IHI=IFIX(Y4{ I }/CY+0«0001 )+2 

IF(I,GT.N51.AND,I .LT,N2) GC TC 375 

IF ( IL .GT .!LI ) GC TC 375 

CG 365 J=IL,ILI 

call *^ESH( I , J ,TNEW ) 

TNEW=(l.-CME)=^T{I, J)+CNE*TNEW 
SUM^SUM+(TN£W-T( I,J) )’S=(TNEW-T( I t J) ) 

65 t(i,j) = tN£Vv 

INNER BCUNCARY - LOWER SURFACE 
375 a2=FLCAT(I-l)¥CX 
C2=Y3( I) 

CALL CUR(S2,C2,B,C,1 ) 

XMN = 8 + 2 .*C^62 
IFd.EC.Nli) XHN = -1.0 
YT = FLCAT(ILI-1)>!'CY 
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OP POOR QUA.Uii 


X={C2-VT)=«'XiMN + B2 

IF(I,EG.7) XI^N=->1,0 

IF( I.GT.N51.AN-C.I.LT.N2) GC TC A12 

IF (X ,LT,X1 ( ILH ) GC TC 380 

IF { (C2-YT) .G£.(C2-Yin > ) ) GG TC 380 

GC TC 382 

3S0 T7T=TBX{ 1 ,ILI ) 

XKL=(C2~YT)*(C2-YT)-}-{Xl( ! L I )-E 2 J O 1 { I L I ) -8 2 ) 

XM = SC8T( XM) 

GC TC ^00 

382 IF { flBS (XN<N ).LE.0.015 3 GC TC 395 
IF(X,LE.B2) GC TC 283 

TTT = T( I , III ) + (T{ I + l, ILD-TI I, IL n (X-62 3 /CX 
GC TC 390 

383 CCNTINU6 
XL=B2-CX 

IF( { XL-Xl (ILI ) ) .GT^G.QOOl) GC TC 385 

TTT = Tn,IUI)-{TtI, ILn-TBX{l,ILI) 3 ^ { B2- X ) / { 82-X 1 ( I L I ) ) 

G G T C 390 

365 TTT=T( I, ILI3-{TCI,!LIJ -1(1-1,111} )*(62-X)/CX 
39<3 XM = {X-B2}*(X-B2)4{C2~YT )*(C2-YT3 
>M=SCPT( Xi\U 
GC TC 400 
395 TTT=T(I,IL1) 

XNL=C2-YT 
400 AL=F(I,3)/XK 

TNEV» = (XM=^AL^TG( I ,3 3+TTT) / { 1» + XM4;5L) 

AIG CCNTINUE 

TN£ls=n.-CF’E3X^TE(I,3 E=^TNEW 

SL'8 = SUN-s-nB{ i ,3)-TN£U)=!'(TB(I ,33-T^EV^) 

T£ ( I ,3) = TNEVi 
412 CCNTINLE 
C 

C CALCULATES NEv'f CCCLANT TEMP USING HEAT BALANCE EGUATIGNS 
C 

IF.U-GT.N5 .A^D.I ,LT,N23 GC TO 43C 
C(T ,3)=H( I,3)=*(TE{ 1,3 )-TG( 1,3 3 )X‘CELS( 1,3 )=>LB 
C( I ,4 3=H( I ,4}*(TE( I ,4 3-TC-{ 1 ,4 3 3^CtLS( I ,43>»Le 
IF{ I .6C.N71 3 GC TC 4 15 
IF(I.LT.N7L3 GC TC '-20 
IF(I.GT,,N71) GC TC 425 

415 TG^ + 1,3} = (^'ASC{I,3 3=^CP=^TG(I,33 + ^'ASC( 1 , 4 ) *CP4TG ( 1 , 4 ) +C (1 , 3 3 +C ( I , 4 3 
1) /(8ASC{ I + l ,3)*CP3 
GC TG 430 

4 20 TG { I + l ,3 3 = (MASC (I ,3 )4CF^TG { 1 ,3 ) + CEL.V (I , 3 ) =^CP^TCOLM+ Q ( 1 , 3 3 3 / ( VASC { 
1I^1,3)=^CP3 
GC TG 430 

425 TGI I + 1,3)==(NASC(I ,3)=^=CPXsTG { I ,33+ C(I,33+ G ( 1 , 4 3 3 / { « ASC (I + 1 , 3 3 ^CP 3 
C LFF6F SLRFACE CF INNSP BCUNCArY 
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c 

^30 C2=V4(n 

C4LL CUR { E2 ,C2, g,C, 1 ) 

XyN = 3+2,=^'C*B2 

IF { I .EG. Nil) XyN=1.2 

Y7 = FLCAT( IhI-l)=}'CY 

IF (A6S( XMN ) .LE,0.ai5 ) GC TC 445 

X=(C2-YT)^xyN+B2 

XR=B2+CX 

XL = B2-C'X 

IF (X.LE.82 ) GC TC 433 
IF ( X.GT*X2( IHI ) ) GC TC 446 
IF { (X-XR) .GT. 0*000 1 ) GC TC 447 
IF( (X2( IFIJ-XPj.LE.O.OOOl) GC TC 434 
TTT = T { I, Ihl ) + (TU + l , IFI )-T(I ,IHI ) (X-S2) /CX 

GC TC 440 

4 24 TTT=T(I ,IHI J + (TBX(2,IHn-T<I , I h IJ 5 ( X~e2 ) / CX 
GC TC 440 
425 CCMIKLE 

IF ( (X-XK IHI) ).LT. G.OOOl) GC TC 448 
IF ( (XL-Xl UFI ) ) .LT*0 .cool ) GC TC 427 
TTT = T(I ,IHI T( I , IHI )-T ( I'l , IFI) (B2-X) /CX 

GC TC 440 

437 TTT = T ( I , I F I }- (T ( 1 , I F I) -T BX 1 1 , I F I ) ) ( E 2-X ) / ( B 2-X 1 II H I )) 
44 C XNL=(X-B2) {X-E2) + (C2-YT J*1C2-V7) 

XM = SCRT(XM) 

GC TC 45Q . 

445 TTT = TU,IFI) 

XM = YT-C2 • 

GC TC 450 

446 IF{ {Y2(I+1 J-YT) .GT.Q.GOOl) GC TC 447 
T7T = TBX{2 f IHI ) 

XNL={b2-X2{ IFI) ) =» ( B2-X2 { IHI ) ) + CY^CY 
XM = SCR1 ( XM ) 

GC TC 450 

447 TTT = T ( I + l , IFI ) 

XM = DX=!=DX+ { C2-YT ) ( C2-YT ) 

XM = SGRT(XNL) 

GC TC 450 

448 TTT = T3X( 1 ,IHI ) 

XKL= (8 2-XU IFI ) )=«■( E2-X1( IF I ) I + CY^i'C Y 
XM = 5CRT (Xi\L ) 

450 AL=H{I,4)/>K 

TN6Vi= (Xi\'L=^AL*TG ( 1 , 4 ) -iTTT } / (1 *+ XM « AL ) 

TNEF, = ( 1 ,-CiVE) -7 8(1 ,4)+CFE*TNEW 
SUiY=SUY+(Te( I ,4)-TNEi^)=S={T8 (I .AJ-TREV' ) 

78 ( I ,4)=TNEV 

CALCULATES NEW CCCLANT TEiVP USING F.EiT BALANCE ECUATIONS 
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IFU.GT.N5 .LT.^2J GC TC 47C 

IF(I.EC.N71) GC TC 455 • 

IF(1,LT.N71) GC TC 460 

!F( UGT.N71 ) GC TC 463 

455 TG{ I + l ,4)==(NASC ( 1,3 )=^-CP^TG U ,3)+r^ASC { 1 1 4 )’!'CP=^TG{ 1 ,4)+ QU,3)+ Qd, 
14)}/(PASC(I+lt4)’>CP> 

GC TC 470 

460 TG( I + l ,4) = (MASC ( I ,4 ) ^CP^i^TG { I ,4) +031^(1 , 4 ) P=?TCCLN+ C { I t 4 ) ) / ( P ASC ( 
ll + l ,4 )^CP J 
GC TC 470 

465 TG( I + i,4) = (PASCU ,4)^CP«TG(I ,4)+ C{I»3}+ C (1 ,4n / ( PASC ( I + 1 ,4 ) ^CP 5 


lATSPIDR POINTS FRCP UPPER SURFACE CF INNER 8CUNDARY TG THE UPPER 
SURFACE CF cuter BCUNCARY 

470 DC 475 J = IhI ,IH 

CALL VESH( I, J,TNE'rt) 

TN£'rt= (1 .-CPE3=«=T n , J 3 +CPE=»TNEW 
SLP = SU«+{TNEls-T( d J) 3*(TNEU-TU ,J3 3 
475 T( I , J )=TNE'a 
4E0 CCNTINLE 

UPPER eCUNCARY CF CUTEP SURFACE 
C2=Y2( I) 

CALL CUR(B2,C2tBtCtG3 
XPN>=e^-2.-C=^62 
IF{I.SQ.NY13 XPN=-1.8 
IFU.EC.16) XPN=-2,8 
IFd.EC.193 XPN = ~2.e 
YT=FLGAT{ IH-1)*CY 
!F( AES(X,MN).LE.0.015) GC TC 500 
X={C2-YT3^XF^Ni + S2 
XR=E2+CX 
XL=82-CX 

IFd.EC.3} XPN = 1.0 

IF( I.LE.NL ,CP.I .GE.N2 3 GC TC 461 

IF( IC2-{Y4(I)+Y4{ 1+13 3/23 .GE.CY3 GC TC 481 

TTT = TS (I t4 ) 

XNL=C2*-Y4(I ) 

GC TG 505 
461 CCNTINLE 

IF(X.LE.32) GC TC 463 
IF(IF.G€.N4.0P. Ih.LE .N23 GC TG 464C 
!F{ (X-X3( IH3 3 ,GT. 0.00013 GC TC 492 
IF(X3I IH3 .LE.XR 3 GC TC 482 
4 64 0 TTT = T{ I, IF ) + (Td + l t IF3-Td, IH3 3*(X-E23 / CX 
GC TC 495 

4 82 TTT = T( I, if ) + {T3X(3 , if )-T{ I , I H 3 3 =* { X- E 2 3 / { X 3 { IF 3-B2 3 
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^83 


^85 




495 


49? 

498 

500 
5C 2 

503 

505 


510 


GC TG 495 
CCMINUE 


ORIGINAL’ PAG® ig 
OP POOR QUAUTJJ 


IF{ Ih.GE,N4.CR. Ih.L£.N3) GC TC 435 

:F( X.LE.X4 { Ih ) ) GG TC 49S 

IF( {XL-X4(II- ) ).GT, 0,0001) GC TC 490 

TTT=T ( I , IH3-( T( I,II- )-TBX (4,IH >.) { B2-X) / ( B2-X4 ( IH) ) 

GC TC 495 

CCMINCe 

IF{X.LS.XlUh) ) GG TQ 497 
IF( {XL~X1(IH)),GT. 0,0001) GC TC 490 
TTT=T{! , Ih)-(T( I , IH }-TBX( 1, IH) ) I B 2-X ) / ( B 2-X 1 ( I H ) } 
GC TC 495 


TTT=T( I , IH)-(T{ I Tlh)-T( I-l,Ih) )^ 1B2-X) /CX 
GC TC 495 
TTT = T3X(3 ,Ih) 

XM= (E'2-X3 ( IH ) (82-X3( IF) I + CY’iCY 

XRL=SCRT{XM) 

GC TC 505 

XM=(X-B2)=i'(X-e2)+{C2-YT )=^‘{C2-Y7) 


XM = SCRT(XM) 
GC TG 505 
TTT^TBX (1 , IF) 
GC TG 503 
TTT=TBX (4, IH) 
GC TC 502 
TTT = T( I , IH ) 


XNL=C2-YT 


GC TC 505 

XM= (8 2-XK IH ) ) {B2-XU IH) )+(C2-YT)^{C2-YT) 
XM = SCRT(XNL) 


AL = H( I ,2) /XK 

TNEW= (XML=^AL^TG( 1,2 ) +TTT ) / ( 1 , + XN L ) 

TNe>A={l,-CH.E)’fTE{ I ,2 ) + CNS»^TNEW 
CLf'=SUM+{TE{ I ,2)-TKEliJ=!'{TB( I ,2)-T^EV^) 
TE( 1,2 ) = TNEis 
CCNTINCE 


ITER. FCR J LINES CMY BCUNCARY POINTS INTERSECTED BY J LINES 

CC 650 J=2,NY1 
IN =IFIX(Xi ( J)/CX+0.0001)+2 
lA =IFIX(X2( J )/CXi^0.CCi01 ) + l 
IF(AES(X2 (J)-FICAT UA-D^CX) ,LT, 0.00001) IA=IA-1 

N-EAR SURFACE CF CUTER ECUNCARY 


B2=X1( J) 

C2 = FLCATIJ-1 )X=CY 
CALL CU’R(R2 ,C2,B,C ,C) 
X,VN = 3 + 2 .=i=C=*a 2 
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XT = FLC.iT( IN-1 )=»CX 
IF(J.EC.25) XNN=-1.0 

IF( J,G£.N31-AN0.J.LE.N6i) GG TC 5A7 
IF (ASS(XiVN-) .GE,5- )GC TC 530 
XL=XT-DX 
X=-0Y^XMN-+8 2 
IF(X.LEo82I GC TC 520 
IF{ J.LT.N6 .CR .J.GE.N^n GC TC 512 
IF(X.GT.X3{ J + 1) 1 GC TO 542 
IF((X3{J+1)-XT),GT,G.0Q01) GC TC 515 

TTT=T( IN-1, J+1)+(TEX( 3,J+1)-T( IN-1,J+13 )*( X-XL) /{ X3( J+1)-XL) 
GC TC 525 
512 CCNTINUE 

IF(X .GT.X2IJ + 1 3 3 GC TC 522 
IF( (X2{J+1 3-XT) -GT. 0.00013 GCTC 515 
TTT = TUN-1,J + 1 3+(TEX{2, J + l)-T( IN-i,J+l) 3 ? ( X-XL3 / ( X 2 ( J + 1 3“XL 3 
GC TC 525 
515 CCiNTINUE 

IF( (XL-X1{ J-H ) 3 .GT.0.0001 3 GG TO 5ia 
IFIX.GT.XT) GO TC 533 

TTT=T( IN, J+1) + <TSX{ 1 , J4i )-T< IN, J + 1 3 3 ^ ( XT- X ) / ( XT-X 1 ( J4 1 } 3 
GC TC 525 
51£ CCNTINLE 

TTT = T{ IN,J413 + {T{ IN-l,J + i3-TnN,J+13 >*{XT-X) /DX 
GC TC 525 
52G CCNTINUE 

X = CV*XNN-i-E2 

IF{XT.GE.X3(J-13) GC TC 524 
IFCX.GT.XT) GO TO 523 

TTT = T{ IN,J-13 + (T( IN-1 ,J-13-T(IN,J-1 3 )^{XT-X3/CX 
GC TC 525 

523 TTT=T{IN,J-13 
GC TC 534 

524 TTT=TBX (3,J-1 3 

XNL=(B2-X3( J-13 3^{E2-X3IJ-1) )+CV^CV 
XM = SGRT(XNL3 
GC TC 535 

525 XNL=(B2-X3=^=(e2-X)+CYX=CY 
XM=SCRT(XNL ) 

GC TC 535 
530 TTT=T(IN,J) 

XNL=XT-B2 
GC TC 535 
532 TTT=TB(IN,23 

XNL=I82-XT 3=^(B2-XT )-K C2-Y2 ( i N 3 > ^ { C 2- Y 2 (I N 3 3 
XNL=SCRT(XNL) 

GC TO 535 

532 TTT=T(IN,J+13 
524 XNL=(32-XT)=^'(B2-XT} + CY*DY 

XM = SCRT(XM. 3 
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535 JiL=l-XC 1 , J J /XK 

TKEW = (TTT + AL=^=XNL’^TGX( 1, J) )/{l. + XM^4L) 
GC TC 5A5 

5^2 TN£V,= (TB( I N , 1) TB ( I K- 1 , 1 ) ) /2 
5A5 CCNTINUE 

TNEl^={ l.-Ci'^Ei^iaXC 1, J) + C^e=^TNEVi 

SUf^ = SU/^ + { TNEW-TcX ( 1 , J I J^CTNEW-TEXl 1,J) ) 

TEX{1,J) = TNEV^ 

5^7 CCNTINLE 

IF(J,GT.N3,ANC.J.LT.N4) GC tq 55 c 
GC TC 625 


FCP J LINES INTERSECTING INNER ECUNCAP V-BGUNCiRY PnTf^TS OM THE 
NEARER SLRFACE CF INNER SCLCARV 


550 CCNTINUE 

INI=IFIX( X3 { J )/CX+C.0001)+l 

IF(ABS(X3 ( J3~FLCAT( IM-1 )=*CX ) .LT.O»OQOOU IN I = IN I- 1 
IAI=IFIX( X4(J) / CX + C-QOOD+2 
B2 = X3( J } 

CALL CIJRC B2 ,C2,E ,C ,1) 

XNN = E+ 2 .>?C*e 2 
IFU.EC.N31) XNN=-1.7 
XT=FLOAT( INI-1)*CX 
XR^XT+CX 
XL=XT-CX 

IF(J.£S. 2 E) XNN= 1 . 2 C 
IF(AeSlXNN) .GE. 5 . )GC TC 565 
IF { J.GE.N 3 .ANC. J. LE-N 61 ) GC TC 587 
IF { {C 2 -IV 1 ( IND+YU INUin/ 2 .) ,LE,CTf J GOTO 570 
X = CV’i'XiVN' + E 2 
IF{X.GT, 62 ) GC TC 555 
IFIX.LE.XK J- 1 ) ) GC TC 570 
IF( (Xl( J-D-XT) .GT.O.OOOi) GC TC 552 
IF { (XT-X) .GT.O.QOOl ) GC TC 562 
TTT=T(INI, J-l ) + {T{ INI + 1 ,J- 1 )-T (INI , J-IJ ) =!= (X-XT) /OX 
GC TC 560 

5 52 TTT = TBX ( 1 , ) + (T ( IN + 1 , J~ 1 ) -T EX ( 1 , J - 1 ) ) * { X-X I ( J- 1 )) / ( XR- X 1 ( J- 1 ) ) 
GC TO 560 
555 X=-CV*XNN+e 2 

IF{ X.LE.XK J + 1) ) GC TC 572 

TTT = T< IiM ,J + 1) + (T( INI, J + D-T (INI+lt J+1) )«( X-XT)/DX 
GC TC 560 

562 IF( { XL-Xl(J-lJ ) .GT.O.OOOI) GC TC 563 

TTT=Tax( 1,J-1)4(TUM ,J-1)-TEX{ i , J- 1)) ( X T-X ) / (XT-X U J-1) ) 

GC TC 560 

56 3 TTT=TnNI-l,J-l)-MT( INI , J-D-T (INI-1 ,J-1) )*(XT-X)/CX 
560 XNL=(22-X)>^ (E2-X }4CV’>CY 
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XM = SCRT( XM) 

GC TC 575 

565 TTT = T(IM,J) 

XKL=B2-XT 
GC TQ 575 
57G CCMINLE 

IFUNI .EQ.N51) GC TC 577 

IF { (XL-Xl ( jn ,GT. 0.0001) GC TC 574 

TT7=(T5(IN! ,1)+TBX(1 ,J) 3/2. 

XNL = B2-(XT+X1 (J) )/2, 

GC TC 575 
572 TTT=TB{INI,2) 

XNL^ (XT-62)=i'{XT-P2 3+{Y2( INI )-C2)^('f2( IN’I )-C2) 
XM = SCRT(XM) 

GC TC 575 

574 TT7 = TB{INia) 

XM = (32-XT3«(B2-XT3-MC2-yi(IM3 3 * I C2-Y 1 { I M ) J 
XM=SCRT(XM 3 

575 AL=HX(3f J)/XK 

TNcV^='(TTT + AL^XNL4TGX{3,J) )/{ l. + XM^AU 
GC TC 585 

577 TNEV( = T6{ IM , i) 

585 CCNTINUE 

TNEV(=( l.-CME)*TBX (3 , J)+C^E’^T^EVt 
SUM = SUN + { TNSV^-TEXI 2 , J) ) =* { 7 iN E W-7 E X ( 5 , J ) ) 
TEX{3,J3 = T'NEW 
567 CCNTINUE 

CC 590 I = M,N2 
XR = FLDAT( I-l )^CX 
IF{ (XR-X3(J)),GE.C. 00013 133=1-1 

!F( (4ES{XR-X3{ J 3 3 .GE. (0.9595=^CX 3 3 3 1 3 3= 133-1 
IF{ (XR-X3 (J3 3 -GE .0.0001 ) GC TC 553 
55G CCNTINUE 

593 IF( I33.lt. N il) I33 = Mi 
TGX(3»J)=TG{ I33 t3) 

FARTHER SURFACE CF INNER BCUNCARY 

ez=x4( j) 

CALL CUR ( 62tC2,e ,C ,1 3 
XNN=e+2,^C«E2 
IF-(J.Ee.N31) XNN=-1.7 

XT = FLCAT( I AI-13=^CX 
XL=XT-DX 

IF (A3S(XNN) .GE .5. 3GC TC 6G5 
X = -C-Y^XMN + EZ 

IF( (X-XT) .GT.O.OOQl) GC TC 595 
IF(X.LT.X2( J + 13 ) GC TC 557 
595 CCNTINLE 

IF ( (Y2( IAI3-(C2+CY) 3 .GT. 0,00013 GC TC 596 
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TTT=TBX(2,J+1) 

XNL=( B2-X2 { J+1 ) i’f (E2-X2 (J + 1) ) + CY*CY 
XM = SCRT(XNL) 

GC TC 610 

596 1TT=T( lA! ,J + 1) 

XM=(62-XT ) + (e2-XT)-5-CY*CY 
XM = £GPT{ XM ) 

GC TC 610 

597 CCNTINUE 

IF(X2( J^l) .LE.XT) GC TC 598 
TTT = T(IAItJ + l)-KT{ I A I - 1 , J + 1 ) -T ( 1 4 I , J+ 1 ) ) =*> ( XT- X ) /D X 
GC TC 600 

5 98 TTT=TBX(2 , J+1J + (T( IAI-1, J+1)-TEX{2,J+I)l^a2( J + i)-X)/ 1X2 ( J + 1)-XL) 

600 XM= (e2-X { E2-X )4CV^CY 

XM = SCRT{ XNU 
GC TG 61C 
605 TTT = TUAI,J) 

XM = XT-62 

610 AL = KX(A, J )/XK 

TNc^»=(TTT 4AL*XNL^TGX (4,J))/tl.+XM=»AU 

6 12 TN-EW= ( 1,-CME )X=TEX{A, J)+G^'E^T^Ei^ 

SLN = SUf^+{TN£h-TEX14, J) ( TNEW-T EX ( A , J ) ) 

TEX (A, J 3 = TNE!^ 

CC 615 1=133, N2 
XP=FLCAT( I-13*0X 
IF ( (XP-X4(J ) ) .GE.O.OOOn IA=I 
IF (AES(XR-XA( J ) ) ,LE,0.0001J IA=IA+1 
IF 1 (XP-X4( J J ) .GE^C.QOOl ) GC TD 620 
615 CCMINtE 

620 IF1I4.GT.N21) I4=A21 

TGX(4,J)=TG1I4,4) 

C 

C FARTHER BOUNDARY FCUTS CN CUTER SURFACE 
C 

625 32=X2(J) 

CALL CUR(E2,C2, e,C,G) 

XNN = E42.x<C=^e2 
IF1J.EC.15) XNN=-2.0 

IFIJ.EG.IT) XNN=-2.C 
IFU.EG.2 ) XNR = 6.0 
■ XT = FLCAT1 IA-1)X<CX 
IF(AES(XNN).GE.5,) GC TC 645 
XP=XT+CX 
XL=XT-CX 
X=CY*XNN+E2 

IF(J.LE.N31J GC TC 635 
IF! lA .GT.N21 ) GC TC 62? 

IF! {C2-Y4 1 lA) ) .GT.CY) GC TC 627 
T7T = TB{ lA ,4 ) 

XM=(C2”Y4(IA) (C2-Y4 1 lA J )+ {B2-XT )^ { S2-XT ) 
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)!M = SCRT( XNL) 

GC TC 6^6 
627 CGNTINUG 

IF{ (XT-X4(J-n ) .GT. 0.00001) GO TC 62G 
IF { (X4(J-1)-X) .GT .0.0001 ) GC TC 623 
IF{ (XP-X2(J-1) ) .IT. 0.0001) GC TC 62S 

TTT=Tex(4,J-l)-{TEX(A,J-l)-TBX(2,J-l) )^( X-X4I J-1) )/(X2( J-l)-X4 
1 (J-13) 

GC TO 642 

62 8 TTT = TBX(4,J-U-{TeX(4 , J-l)-T( 135 + 1, J- 1 ) ) * ( X-X4 { J-1 ) ) /( XR-X4 { J-1) ) 
GC TC 642 
62S CCNTINUE 

IF(X.GT.XT) GC 7C 6AG 

630 CONTINUE 

IF( (X-X4( J-1) ) .GT.O.OQOU GC TC 632 

631 TTT=T ( lA, J-1 ) 

XNL = {e2-XT)=?(B2-XT J + (CV=»DY) 

XNL=SCRT{ XNL) 

GC TC 646 
622 CGNTTNUS 

IF { (XL-X4{J-U ) .GT.G.GGOl ) GC TC 634 

633 TTT:=TEX(4, J-1 ) 

XNL=(S2-X4 (J-1) )*{B2-X4 ( J-1 ) ) + CY=i=CV 
XNL=SCRT (XNL) 

GC IC 646 

634 TTT = T( I A- 1 , J-1 ) + (T ( I A-1 , J-i ) ^T ( I A , J-1 ) )^(X-XL)/CX 
IF { (XL-X) .GT.C.OGCI } GC TC 631 

GC TC 642 

635 CCNTINUE 

IF( IC2-YI IIAJ) .GT.CY ) GG TC 636 
TTT = TB( lA ,1) 

XNL=(C2-Y1{ lA) )^<C2-Y1( lA ) ) + ( B 2-XT ) ( 8 2-XT ) 

XM = SCRT( XNL) 

GC TO 646 

636 CCNTINLE 

IFU XT-XK J-in .GT.C.OOCOl) GC TC 639 
IF ( (XI ( J-U-X ) .GT.C.QOOl ) GC TC 644 
IF{(XR-X2( J-1) 3 .LT, 0.00013 GC TC 638 

TTT = TBX(1 ,J-1)-(T6X( i , J- 1 ) -T B X ( 2 , J- 1 3 ) ^,{ X-X 1 ( J- 1 } ) / ( X 2 ( J- 1 )-X 1 
1 (J-D5 
GC TG 642 

638 TTT=TBX(l,J-i)- (T6X( 1,J-1)- T(IA+1,J-1))^(X-X1(J-1)J/(XR-XI(J-1)3 
GC IC 642 
635 CCNTINUE 

IFIX.GT.XT} GC TC 640 
IF( (X-XK J-1) ) .GT.C.QQQl) GC TC 643 
GC TC 631 
64C CCNTINUE 

IF( [XP-X2(J-1)).LT.0,0C013 GC TC 641 

TTT=T(IA, J-1) + (TSX(2,J-1)-T( lA, J-1) ) * ( X-XT 3 / ( X2 ( J-1)- XT 3 
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GC TC 642 

641 TTT=T{ lA ,J-L) + (TnA + l,J-l)-7n^ ,J-1) )=^={X-XT) /CX 

642 X^•L=(E2-X)’^'(E2-X)-^CX’^CY 
XM=SCRT{ XNU 

GC TC 646 

643 CCNTINUE 

IF( (XL-XK J-1) ) .GT.C.OQQU GC TC 634 

644 TTT=TBX( 1 ,J~1) 

XN'L=(B2-X1( J-L) )=!'{E2-X1( J-1) ) + CV^CY 
XNL = SCRT(XNL ) 

GC TC 646 

645 TTT=T{IA,J} ' 

XM=B2-XT 

646 AL = HX(2» J )/XK 

TNE^ = ( TTT+AL=S'XM=^'TGX(2,J) )/(i- + XM^AU 

647 CCNTINUE 

TNE^ = {l.-C,ME)*TeX(2, J)+CN8 =^TNE'a 
SUN=SUM+(TBX< 2, J)-TNEV»)=l'(TEX (2, J)-TN£W) 

TBX{2,J)=TNEU 
650 CCNTIME 

CG 655 I=1,NX1 
CC 655 J=1,NY1 
X = FLGAT{ I-1)=^'CX 
TE{l,l)=TeX(ltJ) 

Y=FLCAT( J-1)«CY 

IF{ABS(X-XUJ)),L7 .0 .00001. AND. ABS ( Y-Y 1 ( I > ) , LT. G.OOOO 1 } T(I , J ) = T8 < 

11,1) 

IF(A8S{X-XU J) ) .L7.C.0OOQ1.AND. ABS (Y-Y2( I) ) . LT .0 ,0000 1 ) TH,J )=TB{ 

11,2) 

IF {ABS(X-X2( J) ) -LT.Q. 00001. ANC. ABS (Y-Y2( !J ) , LT . 0 .0000 1 ) T{I,J)=TB{ 

11,2) 

655 CCNTINUE 

CC 657 T = M1,N21 
CC 657 J=N31,N41 
X=FLCAT{ I-1)^CX 
Y = FLGAT( J-1)=<'CY 

IF (ABS(X-X3 ( J) ) - LT .0 .00001 .AN C , ABS (Y-Y3(I ) ) . LT . 0. GOCO 1 ) Id , J)=TB ( 

II, 2) 

IF (ABS (X-X4( J n .LT.Q,0QG<U.AND.ABS(Y-Y4{ I ) ) . LT . 0 -0000 1 ) T( ! , J ) = TB ( 
11,4) 

IF ( ABS {X-X3( J) ) -LT.O -COOOl -AND . A BS ( Y-Y4 ( I ) ) .LT. 0,00001) T(I , J ) =TB ( 
11,4) 

657 CCNTINUE 


CCNFLETIGN CF AN ITEFAT ICN-CFECKS FCR PECUIREC CONVERGENCE 


IF { ITER.EC.l ) WPITE{6,660) 
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66C FCR.v^T(lHl) 

WRITE(6,665 ) ITER,SU^^ 

£65 FCRNAT{5X , ’ ITER^TICN=» ,I3,3X, ’EPPCF= » ,£12. 5 ) 

IF( ITER.GE, If-<AX > C-C TC 2CC0' 

IF{SLN.LE.SUHV) GC TC 670 

IF( {!TER/NTt*K7E) .EG.ITER) GC TC 67C 

IT£R=ITER+1 

GC TC 300 

67C CCNTINUE 

WRITE (6,7C0 ) 

70C F0R^'AT(// ,10X,» T-RATPIX, GIVES ALL UTERICR PC INI S . I MT I AL SETTING 

1 CF T IS 1600-0* /lOXt >CNLY THCSE INSIDE PLACE ARE CHANGEC AND WILL 

2 8E DIFFERENT FRCP 1600.00') 

WRITE(6,7Q53 

705 FGRPAT (/ ,2X, • J = ' , ^X , • 1 ’ , 7X , ' 2 ' , 7X , ’ 3 ' , 7 X , • 4 ' ,7X , ' 5 * , 7X , ' 6 ' , 7 X , '7 
1' ,7X 8' ,7X, 'S' ,6X 10' ,6X, • 11' ,6X 12' ,6X , ' 13' ,6X, • 14' ,6X, * 15'/ 

22X, ' I ' ) 

IF(NY.GT.15) WRIT£{£,7IQ) 

710 FCBMATtaX,' 16',6X,'17',6X,»18',6X,'19’,6X,'20* ,6X,'21',6X,'22',6X 
I, ' 23 * ,6X, *24 »,6X , ' 25 ' ) 

CC 715. 1 = 1, NX 

715 WPITE(6,72C-} I, [ T (I , J ) , J= 1 , N Y ) 

72 0 FCRP.AT(/,2X,I2, iX,15FE-2,2(/5X,15F8.2)3 
WRITE(6,725) 

725 FCRMATI ' 1 ' ,y/15X, 'GAS TEPFEPATURE CISTRIBUTION ON THE BLADE INNER 
lANC CUTER SURFACE' ,// ) 

WRITE(6,745) NXl 
WRITE(6,720) ( TG { I , 1 } , I = 1 ,N X 1 ) 

WPIT£{6,750) NXl 
WPITEI6,720) { TG { 1 ,2 ) , I=i ,NX1) 

WRITE (6,755 ) M1,N21 
WRITE (6, a OG) ( TG( I ,3) ,I=M1 ,N21) 

WFITEC£,760 ) M1,N21 
WRITE(6,8CQ) ( TG { I ,4 ) , I =M 1 , N2 1 ) 

WRITE (6,770 NYi 
WRITE(6,8G0) (TGX(1,J) ,J = 2,NY1) 

WRITE (6 ,7SG ) NYI 

WPITE(6,800) (TGX (2 , J) ,J = 2 ,NY1) 

WRIT£(6,785) N31,N41 

WRITE (6,800) (TGX(3, J),J=N31,N41) 

WRIT£(6,7S0) N31,N41 

WRITE (6, 800) (TGX(4, J) ,J=N31,N41) 

WP ITE(6 ,7S5 ) 

V<jPITE (6,740 ) 

740 FOPP.ATI ' 1 » ,// ,5X , 'TEPFERATUPE CISTRI8UTICN ON THE BCUNDARY • , ////) 
WRIT£(6,7^5) NXl 

745 FORMAT!//, lOX, ' EGUNCARY FCIN'TS CN T-LTNES FOR LOWER SLRFACE CF 
1 CUTER eCUNCPV' /,12X, ' (START FRCP 1= 2 TC I = ',I3,' )') 

WR!T£(6,720) { T E ( I , 1 ) , I = 1 , NX 1 ) 

WR!T£(6,7£Q) NXi 
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750 FCPVAT(//,10X, » 8CU^CABY FCIKTS CN I LIN6S FOR UPPER SURFACE OF 
1 CUTER BOUMCRV* / ,12)« » M START FRCf' 1= 2 TC I = *,I3*' )M 
VtPITE{6,72G) (Te(I,2),!-lrNXl) 

VFITE (6,755 ) M1,A21 

755 FOPMAT{//,iOX, 'BCINCAFY PCUTS CN I LINES FCR LOWER SURFACE CF 
1 INNER BCU.NCARY»/,12X, ' (START FRQP TG I = >«) 

'aRI TE(6,8G0) (Te(I,3) ,I^M 1 ,N 21 ) 

WRrTE(6,76C) K11,N21 

76G FCRMAT{ //,10X, ' eCLNCAPY FCINTS CN I LINES FCR UPPER SURFACE OF 
1 INNER a0UNCARY*/,a2X, ‘(START FRCP I=»,I3,» TC )’) 

WRITE{6,80Q) (TB ( I,A ) ,I = Mi,N21J 
WRITE (6,770) NYl 

770 FCPNAT{//,lQX,'eCL-NCARY POINTS ON J LINES FOR NEARER SURFACE OF 
1 CUTER 80UNCARY ‘/,12X, '( START FRCN J= 2 TC J=',I3,‘ )’) 

WRITE (6,800) (TBX( 1 , J ) , J = 2 ,NY 1 ) 

WRITE(6,780) NYi 

7S0 FORHAT(//,10X, ‘6CUNCAFY POINTS CN J LINES FOR FARTHER SURFACE OF 
1 CLTER EGUNCARY' / ,12X, • (START FRCP J= 2 TO J=*,I3,» )») 
WPIT£(6,800)(TBX(2,J) ,J = 2,NY1) 

WFITE(6,785) N31,NA1 

785 FORMAT!//, lOX, ’BGUNCARY FCINTS CN J LINES FCR NEARER SURFACE CF 
lINNER eCUNCARY V12X, '(START FRCM J=‘,I3,» TC J=‘,I3,' )’) 

WRITE (6, 800) (TBX(3,J) ,J = N31,N41) 

WPITE(6,79C) N21,NA1 

790 FORMAT(//,10X,*BGUNCARY FCINTS GN J LINES FOP FARTHER SURFACE OF 
lINNEP BCUNCARY. »/l2X ,' (START FPCN J=»,r3,' TO J=»,I3,« )») 

WRITE (6,8D0) (T8X ( A , J ) , J = N31 t NA 1 3 

WRITE(6,795) 

795 FCRMAT(lHl) 

800 FGPMAT{3(/ ,5X .15F8.2 ) ) 

WRITE (6,810) 

810 FORMAT ( IHl/// ,I5X, 'TENFERATURE CISTRIBUTION INSICE THE BLADE AND 0 
XM ITS SURFACEV//2X, ' J=» ) 

DC 885 JL=2,NY1 
CQ 815 1 = 1, NX 
J=NY+1-J1 
XR=FLCAT(I-1)=^DX 

IF( (XR-XKJ) ) ,G£,C.Q0Q1 ) 11=1 
IF (ABS{XR-X1 (J) ) .LE.O.OOOl ) 11=1+1 
IF( (XR-XKJ) ).G£. 0.0001) GC TC 82C 
815 CGNTINUE 
820 DC 825 I=IL,NX 

XR=FLCAT(I-1)*DX 

IF( (XR-X2(J) ) .GE.C.OOOl) 12=1-1 
IF (ABS(XR-X2( J) ) .GE. (C. 9999^0 X) ) 12=12-1 
IF{ (XP-X2(J )} .G£. 0.0001 ) GC TC 830 
825 CCNTINLE 
830 CCNTINUE 
J=NY+1-J1 

IF( J.GE.NA.CR. J.LE.N3) GC TC 86G 
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835 CCNTINUS 

CC 8^0 I=K1»N2 
J=NV+1-J1 

XR. = FL0.AT( I-D’^CX 
IFC (XR-X3 ( J) ) -G6,O.OOOU 133=1-1 
IF( {AeS(XP-X3(J)).GE.(0,<;S99*CX)) ) 133 = 133-1 
IF( (XR-X3 ( J > ) ,GE,Q.Q001 ) GC TC 8A5 
£^C CCNTiNUe 
845 CQ 850 I=I33»N2 
XR=FLCAT( I-L)«CX 
IF ( (XR-X4 (J ) ) .GE .C-OCOl 1 14=1 

IF(A8S(XR-X4( J) ).LE. 0.000 11 14=14+1 
IF((XP-X4( J)1 .GE.C.aOOi) GC TC 855 
85C CCNTINLE 
£55 CCNTINie 

J= NY+l-Jl 

’rtP.ITE{6,87 5) J,TBX(ltJ3 ,I = U,I33),TeX(3,J) 

WRITE{6,880) Tex{4,J),(T(I,J},I=I4,T2),TEX(2TJ) 

GC TC 865 

£60 WRIT£(6,870> J ^ TBX { 1 » J J , ( T { I , J 1 , 1 = 1 1 , 1 2 3 ♦ T8X { 2 ^ J ) 

865 CCNTTNUE 

£7C- FCPf^AT (/2X,I2,4X^15F8.2) 

£75 FGFI^AT{/2X, I2,4X»15F8.2) 

8 80 FGFMAT</16X,‘ »3X,12Fa.21 

885 CONTINUE 

IF(SUM.LE.SUiVN) GC TC 890 
ITER=ITER+1 
GC TO 300 
890 CONTINUE 

IF{ABS(TENPC-TG(K51,3)}.LE.5) GC TC 910 


^'£i^ VALUE OF FUN CCCIANT TENF IS CET4INEC 
TFEN GOES BACK TG GET FINAL SET CF TENFS 

TENFC = TG{iN51 ,3) 

CENSC=i98.57/TG(N51,3) 

CFC=(9.87-36,1/ {TG(N51 ,3 ) . 5 1 -2389 ./TG ( N 5 1 , 3 ) + 9 . C6E C5 / (TG { I , J > 

12. ) )/28.97 

VISCC=7,77E-07* (TG <N51 ,3 )3=^*1.5/{TG(N5l7 3 3 + 198,3 
REASSIGNING GAS TENFS CCV.NSTREAN CF SLCT TO CRIG. VALUES 


CC 900 I=N52»NX1 

TG(I ,1)=TGTEN(I 3 
900 CONTINUE 

CC 905 J=2tN6 
TGXLl , J3=TGXTEN ( J) 
905 CONTINUE 
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GC TC 245 
910 CCNTIKUE 


CSTAINS ISCThEP/^;5L LU6 -LCC^TIGNS 
940 V.PITE {6,945 ) 

945 PGRi^AK IHl ,38X, • ISCTHEPNAL LIAE LCCATICNS ' / ,9X , ' I * tA.X, » J ’ , 3X , » T • , 
14X,'T - l»,6X,»T - 2 SE)!, 'FRAC* ,/) 

CALL TE.^'F 
990 ■ CCNTIMUe 

WPITE{6,S92) 

992 FCRf-'ATdX , IHl) 

CC 999 NN=l,NSCLT 
V^RITE (6,995 ) TSCLT(^:^) 

995 FCPiVAT (IX, 'THE ISCTHERH AT ',I5,/1 
ViFITE (6,997 ) 

997 FCPHAT(5X, ’ X-COGRC « , IGX , ‘ 7-CCCPC »,/) 

,NNZ = NNY(NN ) 

WRITE (6,99 8 ) ( ( X I S C { NN ,NM ) , Y I S C ( NA , NNL ) ) , fvNL- 1 , NNZ ) 

998 FCFyAT(5X,F8-4 ,1CX,F8.4} 

999 CCMINUE 

CC lOOO I=^52,^X1 
TG( I , 1 )=TGTE.M ( I } 

IGGG CCNTINUE 

CG 1010 J=2,N6 
TGX (1,J)=TGXTE^(J) 

ICIC CCMUv'UE 

CG 680 1=2, NXl 

EFF( I ,1 ) = (TG( I, D-TB ( I, 1 ) i/(TGU, 1 )-TCCLN ) 

EFF( I ,2)=(TG{I ,23-T8{ 1,2) )/ (TG( I ,2)-TCCL^) 

680 CCNTINUE 

CC 685 J=2,NY1 

EFF (1, J) = (TGX( 1 ,J) -TSX (1,J) )/(TGX ( 1 ,J)-TCCLNJ 
EFF (2, J )= (TGX (2,J )-T 6X ( 2 , J ) ) / ( TG X ( 2 , J ) -TCGLN ) 

685 CCNTINUE 

WRITE{6,687) 

687 FCRMAT( 'I',//, IflX, 'EFFICIENCY CF CCCLI NG ',//// ) 

WPI'TE(6 ,745) NXl 
WRITF(6,690)(EFF(I,1),I=1,NX1) 

6 9C' FCRNAT(/,2X,, I2,1X,15F8.4,2(/5X,15F8.4) ) 

WRITE(6,75o')NX1 

WPI 7E (6,690 ) ( FFF(I,2),I=l,NXl) 

WPITE(6,770) NYl 
WPITE{6,692)(£FF(i,J),J=2,NYl) 

692 FCPNATO(/,5X,15F8.4 ) ) 

WPITr(6,780) NYl 

WRITE (6,6 92 ) ( EFF ( 2 , J ) , J = 2 , NY J ) 

20QG CCNTINUE 
STCF 


SNC 
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C FCRTRySN CCNPUTER PROGR>ilw TO OBTAIN TEWPERATL-RE DISTRIBUTION 

C A BLADE USING ThE FINITE DIFFERENCE VETHCC 

C 
C 

C MAIN FRCGRAf^ NC.2 

C 

C FCR CCCLA^T DISCHARGE CN UPPER OUTER SURFACE USE 

C MAIN PROGRAM NC.2 ALONG Vi I TH ALL THE SUERCUTINES 

C 
C 

DIMENSION >CI£0(3C,3C) ,YIEC(30,30) .TS CLT { 30 ) . NKY {30 I 
C IMENS ION DELS! 3 0.4 ) , MASC { 30 , A ) . G ( 30 ♦ 4) .DELM{ 30 .4) 

DIMENSION T (30 ,o5).TE{30,4).T6X{4.65).TG{30,4J,T0X(4,6£). 
lH(3Q,4),HX{4,e£),YU30),'lr5{3C).Y3(33) ,Y4(30) ,X1 (65 > .X2 (65 ) . 

2X3 (65 ) »X4 (65 ) , ex (200 } , 6Y (200 ), ex I( 200 ) .8YI( 200 ) , IP( 200 ) 

DIMENSION TGTEM( 30 . TC-XTEM(65 ) 

DIMENSION VELG(3C,65) 

DIMENSION EFF(30,65> 

REAL LB, MA SC 
INTEGER^A TSCLT 

CCMMCN T .TE .TEX ,TC,TGX ,H .HX ,Y1 , Y2, Y3 , Y4 ,X1 , X2,X3,X4, BX.BY. 8X I.BYI 
X, IP 

CCMMCN CX,CY.NX,NY,NXC,NXl.ITER.M,N2,N3,N4,NP.N£,Ne,N7 
CCMMON/DATy NXl ,NY1 ,M 1 ,N21 , N3 1 , N4 1 . N51 , N52 , N6 I , N62, N71.NSR, 

1 C ELM, MAS C, CELS. XK ,OME ,SU MM , IMAX.NTE, TCOLN ,LS .CP 
CCMMCN/TEtVy VELG .PRNC.TEMFC.CPC.VISCC. ALPHA, FLO C, DENSG 
CCMMONVI SC/ XI SO , YI SO , TSCLT ,rsSCLT ,NNY 
C;LL CAT a 
245 CCNTINLE 
C 

C TEMFCFARY STORAGE FCR GAS TEMPS DOWNSTREAM OF SLOT 

C 

CC 256 I=N51,NX1 
TGTEMI I ) =TG { I ,2} 

256 CCNTINUE 

DC 258 J=2,N62 

TG>T£M{ J) =7GX(2, J) 

25E CCNTINUE 
C 

C CALCULATES ACIAEATIC WALL TEMPERATURES DOWNSTREAM CF SLOT 
C 

266 CC 267 I-N51 ,NX1 

CALL TEMADI I ,2 .TAD) 

267 TG(I,2)=TAC 

DC 268 J -2 ,N62 
CALL TEMAD ( 2 , J, TAD) 

268 TGX(2,J)=TAC 

DC 26= I = N51 ,N21 
T C- { I . 3 )=T 0 ( 1.2) 

26S CCNTINLE 
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POOR QUAufey 


3CC CCNTlNt'E 
C 

C ST;iFTIN€ of G/^USS-SE IDELL NEThCC CF ITERATION ON i'lINES 
C 

Stt^=0 t 

DC 510 1=2, ^X1 

IL= IF IX( Y 1( J )yDY+G. OOC 1 J +2 

IF= IFJX(Y2( I)/DY^C.000 l) + l 

IF<AB£(Y2(I)-FLCAT(IH-1 )=»CY).LT.O .00001 ) IF=Ih-l 
C 

C CL7EF LCitEP ECUNCAFY 
C 

C CFFCS WITF NECA COMMENT TO SPECIFY SLOPES AT SPECIAL PCIMS 

C 

C 

E2=FL0AT( I-l )5,OX 
C2=Y1 ( I ) 

XP=B2+DX 

CALL CUF (E2, C2 , E,C, 0 ) 

>NN=B+2,^C*e2 
IFII.EQ.NXn XMN=C,01 
YT=FLGAT ( IL-1 )=i:CY 
>=(C2-YT)* >MN+82 
IF ( I .EQ .6 ) XMN=-1,0 

IF(IL.LE.N3) GC TQ 305 
IFCI.GE.NG) GO TO 305 
IF (X . CT .X3 ( ID ) GO TO 330 
305 I F{ ABS( XNM .LE .0 .01 5) GC TC335 
1F(X .LE.E2 IGO TO 320 
>R=82+DX 

IF(IL.LE.N3) GG TO 31C 
IF { I .GE .N2 ) GC- TC 3 10 

IF ( ()<3( ID -XR) , GT.O. OCOl ) GC TC 310 

TTT=T ( I, IL )+{T3X{3, IL )-T{ I , ID )* ( X-B2}/( X3U D-B2) 

GC TC 325 

31C IF ( { X2( IL )- >R } .GT.O. OCCl J GC TC 315 

TTT=T{ I,IL) + (TBX{2,IL)-TU,IL))^(X-B2)/(X2{IL )-82) 

GC TC 325 
3 15 CONTINUE 

IF ( (X-XR) .GT .0 .0001 ) GO TO 337 
TT7=T(I,IL)+{T(J + 1,IL)-T(I,1L))^{X-E2>/DX 
GC TO 325 

3 20 TTT=T{I,IL)-{T{ I,IL)-T(I-l,IL)) ? <E2~X )/DX 
325 >NL=( X-82)^(X-B2)+{C2-YT)^{C2-YT) 

XNL=SCRT (XNL ) 

GC TC 340 

320 TTT=TEX( 3, ID 

XFL=<X3 ( ID-E2 )■* (X3 ( IL )-E2 ) 4 { C2- YT ) ^ { C 2- Y T ) 

>NL=£GRT( >NL ) 

GC TO 340 
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r, n 


325 


/ 


c4C 


3 £C 


17T=T( I ,IL) 




>M_=YT->C2 
GC TC 340 
m = T(I+l,IL) 

X^L=CX^DX■f C C2-YT ) *( C2- YT ) 

> ^L-£GRT (X^L ) ' 

4L = H( I , 1 ) /XK 

TNEW = (XhL’!f /SL^TG (1,1 ) + TTT ) / ( 1 .+XNi_=iAL ) 
T^EVl=( 1 . -C^El^TO (I,1)+C^E^T^EV^ 
SUA'=SUM-KTE ( I , 1 }-TNE A )» ( TB ( I , 1 ) - TNE W J 
T£ ( I ,1 } =TNEW 

IF( I.6T.M .AND.I oLT,N2) GC TC 360 


C 

C 

C 

C 

C 


IMEPICR POINTS STARTING FfiCN LG^EF END FOR I LINES WhICh DC NOT CUT 
TEE INNER ECUNCARY 


CC 2££ J=IL , IH 
CALL MESM I,J,TNEW) 

TNEW=(1 .-CVE)«T( I , J ) + CAE4TNEW 
£LM=£LM+( TNE1^-T(I,J))*{T^EIA-T(I,J)) 


355 T ( I, J )=TNEW 
GC TC 480 
C 
C 

C FCF I LINES intersecting INNER EGUNCARY-ALL INTERIOR POINTS FROM 

C LOViEP SERFAGE OF OLTER BCUNDAPY TC LCwER SURFACE CF INNER ECUNCARY 

C 

360 CCNTINLE 

IL r=IFIX( Y3( I )/0 Y+0. 0 0 01 )+l 

IF < AES (Y3 C I )-FLC AT ( IL I-I J^CY ) .LT « C . G GOO 1 ) ILI = ILI-1 
IHI = IFIX{ Y4(I)yDY + 0,OOCl )+2 
IF { IL .GT . IL 1 ) GC TO 375 
CC 365 J=IL,ILI 
CALL MESH{ I , J ,TNE1A) 

TNEW= (1 .-OME ) =?T( I, J ) +01/E =»TNEW 
£L.v=£LM+{TNEW-T( I,J))T^^T^£W-T(I,J)> 


365 T ( I , J )-TNEW 


INNER BCUNCARY - LCUER SURFACE 
C 

375 E2 = FLCAT ( I-l I’i'CX 
C£=Y3{ I ) 

CALL CUR (E2 ,C2,S,C, 1 ) 
X^^=B+2,^C-XE2 
IF(I.EQ.Ml) XMN=-1.C 
YT-FLCAT ( IL I-l )*CY 
X= { C2-YT)4=XNN + e2 
IFII.EQ.?) XMN=-1.0 
I F (X ,LT. XI ( IL I ) ) -GC TC 380 
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IF( (C2-YT) .GE.(C£-Y1( I )) ) GO TO 380 
GC TC 382 

3£C m = TBX( 1 ,I LI ) 

X^L=(C2-YT J=»(C2-YT)+( X H ILI >-B2)*C XKILI )-8 27 
X^L=£GRT(X^L) 

C-C TO 4CC 

382 1F{;>ESCXMN ) .LE .0 .015) GO TO 3^5 
IF<X.LE.B2) GC TC 383 

TTT=T( I,ILI) + (T( I + 1,ILIJ-T(I,ILI))=J:( X-32) /OX 
GC TC 390 

383 CCNTINLE 
XL=E2-CX 

IF ( (XL-Xl { ILI) ) -GT .0 ,000 1 ) GO TO 38 5 

T7T-T( I, ILI )-(T( I,ILI )-TEX( I.ILI) )^{e2-X)/(E2-Xl (ILJ)) 

CC TO 390 

385 m=T(I,ILI)-(T(I. ILI) -T(I-1,1L !))=*■{ E2-X)/CX 

390 XNU={ X-S2)=? ( X-B2) + { C 2- YT )* { C2- V T ) 

X hL=SCFT (XNL) 

GC TO 40 C 
395 TTT=T (I, IL I ) 

XbU=C2-YT 
4 C C AL=h{ I ,3 ) /XK 

TNEV>= {XNL + AL+TG( 1,3 )+TTT )/( i, + XNL^AE) 

41C T^E«={1.-CME)^'TB(I,3)+CN E^TNEW 

SUf^ = SUM+{ TB( I, 3) -TNE'iV) »( 16 < I ,3)-TNE7i) 

T£ { I ,3)-TNEW 
C 

C CALCULATES NEW COOLANT TEMP USING HEAT BALANCE ECUATICNS 
C 

IF(I.GT.NEF) GC TO 43C 
G( I, 3 )=h( 1,3) TB{ I , 3 )- TGI 1,3)) +OELSI I ,3) VLB 
G I I ,4)=h ( I ,4 )4 ( TEI 1,4 )-TG I 1,4) ) ^CELSI I, 4 )^LB 
IF(I.EQ,N7i) GO TO 415 
i F I J ,LT .N71 ) GO TO 420 
IF( I.6T.N71 ) GC TC 425 

415 T 6 (1+1,3 )=-(MASC( I , 3) *CP*TGl I ,3) + MASC ( I ,4) *CP*TG( I ,4)+GIJ,3) + G(I,4) 
1 )/ {[V ASCI I+i ,3 )^CF) 

GC TO 43 C 

420 TCI 1 + 1,3 )=(MASC( 1,3 )»CP*TG{ I, 3)+DELM{ I ,3)*CP*TC0LN+ Cl I ,3 ) ) / I MASC I 
1 I + 1 ,3 y^CP ) 

GC 70 430 

425 TG ( 1 + 1 ,3 (MAS Cl 1,3 ) *CF + rc I I, 3 )+ Q(I,3)+ Q I I , 4 ) ) /< MA SC ( I + 1 ,3 > ^ CP ) 

C 

C UPPER SURFACE OF INNER OCUNCARY 
C 

43C C2=Y4(I) 

CALL CUR (E2,C2,B ,C, 1) 

XMN=B+2#*C*E2 

IF( I.EQ.Nl I ) XMN = 1 ,2 

YT=FLCAT { IF I-i )*DY 
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IF ( I . 6T. NS. ;JNC. I .LT.N2 ) GC TO 452 

IF <466 ( XMN ) .L-E .0 .CIS ) GO TO 445 

>= {C2-Y7 )*X.VN + E2 

>R=B2+DX 

XL=E2-DX 

IFIX.LE.B2) GC TC 435 
IFI X.GT, X2( IHI ) ) GO TO AA6 
IF ( (X-XR) . GT .0 .0001 ) GO TO 447 
I F ( < X2( IHI ) -XR) . LE.O. COO 1 ) GC TC 434 
TTT = T{I. IhI) + (T( I+l,IhI>-T(I»IHI)I*( >-82> /DX 
GC TC 440 

4 34 m=T(I,IHI )+{TBX(2.IHI)-T(I.IhI)I*{ X-B2 )/CX 
GC TO 440 
435 CC^TI^CE 

IF( (X-XK IHi ) J .UT.O. OCCi ) GO TO 448 
IF UXL-Xl UFI ) ) .LT .0 .OCO 1 ) GO TO 427 
TTT = T{I,IHI)-(TU,IHU-T(I-1,IHI ) )4< (E2-X )/CX 
GC TO 440 

437 TTT = T Utihl )-{T( I,IhI )-T EX ( 1 , IFl ) )^ ( E2-X ) / ( 8 2- X 1 ( IH I ) ) 
44C >NL=< X-B2I^( X-B2)-t-(C2-Yl}:*;(C2~yT) 

X^L = SCfiT (XNL ) 

GC TC 450 

445 7TT = T( I , Ihl ) 

XNL=YT-C2 . 

GC TC 45C 

446 IF ( ( Y2 ( 1 + D-YT ) . GT.O . CCO 1 ) GC TC 447 
TTT = TEX(2, IFI) 

XNL=(B2-X2( IHI ) ) *(B2-X2C IHl) )+CY*CY 
XNL=SCRT (XNL ) 

GC TC 450 

447 TTT = T{ 1+ l» IHI ) 

XNL = CX=f CX+ { C2-YT )>?( C2-YT ) 

XNL=£GRT ( XNL) 

GC TO 45C 

445 1 TT-TEX (1 , 1 hi ) 

XNt_=(e2-Xl( IHI))v(B2-XlC IHI)) + CY’^=CY 
XNL=SGRT (XNL ) 

4SC 4L-H ( I , 4 )7XK 

TNEW- (XNL»AL^TG( I , 4) +TTT )y ( 1 . + XNL=t=AL} 

451 CONTINUE 

TNE1a=( 1 .-OME) *TB (I»4) + CNE=*TNEV4 
SLN=SUM+(TE{I. 4 )-TNEW )*( TB ( I ♦ 4 )-TNE'A ) 

TE ( I .4) =TNE(-i 

452 CONTINUE 
C 

C CALCULATES NE'A CCCLANT TENF USING HEAT EALANCE EGUATIGNS 
C 

IF (I.GT.NSR ) GC TC 470 
IF(I.EQ.N71) GO TO 45E 
IF ( I ,LT .N71 3 GC TO 46C 
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ORIGINAL PAGB IS 
OF POOR QUAUnj 


IF(^,GT.^7l) GC TC 465 

455 TG ( I + l ♦ 4 )=( MASC( 1 . 3 ) 7G ( I ,5)+ !^ASC< I t 4) *CF=?TG { I , 4) + C(I,33+ Q(l, 

14))/(NASC<I41,4J=?CP) 

GC TO 47 C 

460 TG{I + l,4)-( y.ASCC 1,4)=»CF^TG( 1.4 ) + DEUM ( I . 4) ¥CP*TCOLN+ Q( I ,4) ) / ( MASC{ 
1 I+l ,4 )^CF> 

GC TQ 470 

463 TG(I+1»4)=7NASC( I.4)?CP^7C( I.43+ Q(I,3)+ Q ( I . 4 ) ) / ( MA SC ( I + 1 , 4 ) ^ CP ) 

C 

JNTERJCfi PCI^TS FRCM UPPER SURFACE OF INNER BOUNDARY TO 7HE UPPER 
C 

470 CCMINUE 

I F { I « GT.. NS • 4NO • I ,UT ,N2 ) GO 70 460 
IF{IHI.G7. IhJ GC TO 480 
IF < YT .GE .72( I } ) GO TO 480 

471 CC 475 J=IFI,IF 
CALL ^ESH( I , TNE\s) 

TNEa= (1 .-CM E T, J ) +GME^TNEW 

SLN=SLM+ (Ti\e\A-T( I,JJ)^(TNEW-T(I,J)) 

475 T(I,J) = TNEV( 

C-C TC 480 

47E 7U ,IH)=7G< N52.3) 

460 continue 
C 

C LPFER BCUNCAfiY CF CL-TEF SURFACE 
C 

C2=Y2 ( I ) 

CALL CLR(B2.C2 .B.C ,0) 

X(^N=E+2 E2 

IFd.EQ.iej XiVN=~2<.8 
IF(I.EQ.IS) XMN=-2. a 
YT = FLCAT (IF-1 )^CY 
IF (ABS( XMN ) ,LE oO aOl 5 ) GO TC 500 
X= ( C2-YT )^XMN4B2 
X P = E2 + CX 
>L=8 2-D X 

IF{I.£Q.3) XMN=1.0 

IF( I . GT . N£ . AND. I . LT . N2 ) GC TO SOS 
IF( I .LE. N 1 . CR . I . GE. N2) GC 7C 481 
I F ( (C2- ( Y4 ( I ) +Y4 ( I+l ) )y2 ) .C-E «DY ) 60 TO 481 

777=TE( 1,4) 

XNL=C£-Y4{ I ) 

GC TC 505 

461 CCNTINLE 
1F(X.LE.E2) GO TO 485 

IF( IH.GE , N4 .CR. Ih.LE. N33 GC TC 4840 
IF( { X-X3{ Ih ) ) .GT. C. OCC 1) GG TO 492 
I F(X3 < Ih) .LE .XK } GO TC 4E2 
4E4C TTT=T(I,lh)+{T{I + l ,Ih)-T(I,IF))^(X-E2)/CX 
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483 
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4S0 

452 

4 9 5 

497 

498 

500 
£ C2 

5C3 

505 

505 

5 1 0 
C 

C 

c 


c 

c 

c 


OEIGINAi; PAGE Q 
OP POOR QUALITY 


GC TO 495 

T17 = T(I., IH)+(TEX<3,IH)-T(I,Ih))=¥(X-E2)/(X3{H-^-BZ) 

GO TO 4S£ 

CCM I NUE 

IF{ IH.GE.N4.CR. IH.LE.N3) C- C TC 485 

I F (X .LE «X^ ( Ih ) ) GO TO 498 

IF ( (XL-X4( IF ) ) ,GT ,0 .0001 > GC TC 490 

TTT=T{I,Ih)-<T(I,IH)-TEX(4,IH))« ( E2-X) X ( E2-X4 ( IF ) ) 

C- C T C 49 5 
CCNTINLE 

IF( X .LE. X 1 ( IH) ) GO TO 4 97 
IF { (XL-X 1 { Ihl ) .GT .0.0001 > GO TO 490 
TTT=T( I ,Ih)-( T( I , IHl-TEX (1 , Ih ) (B2-X )/ {E2-X 1 { Ih) ) 

GO TC 49£ 

Tn = T(I,IF)-(T(I,IH >~T (1-1,H-))7CE2-X)/CX 
GO TO 455 
TTT=TEX( 3, IF ) 

XNL=( B2-X3 { IHJ )* (E2-X3 ( I M )+CY^CY 
XNL=SQRT( XKL ) 

GC TG 505 

XNL=( X-8 2) 4 : ( X-82)+( C2- YT ( C2-YT ) 

XNL-SGRT (XNL ) 

GC TG 505 
Tn=TBX( l,Ih) 

C- C T G 50 3 
I 7T=TBX ( 4 . I K1 
GC TO 50Z 
TTT=T ( I , IF ) 

X^L=C2- YT 
GC TO 505 

>FL=(E2-X1 ( IF) )=?<e2-Xl (IF } ) + { CZ-YT ) t ( C2-YT > 

XNL=5QRT( XNL) 

^L=F( 1,2 )/XK 

^^EVl=(X^L^4L’^TG{ I ,2 )+TTT ) / (1 ,+XNL=!'4E ) 

CONTINUE 

TNE\^=(l.-C^iE)4=Te(I.2)4GME#TNEV^ 

£LN = £CW+<Te(I,2) -TNEVi ) ( T E ( I ,2 ) -T N E W ) 

T£{ I, Z )-TNEVi 
CONTINUE 

ITER. FOR J LINES ONLY BOUNDARY POINTS INTERSECTED BY J LINES 

DC €50 J=2,NY1 
IF =IF IX (X 1 ( J )/DX+0 .OGCl 3+2 
14 =IFIX<X2 (J) /OX + O.OOCl ) + l 
IF ( AES { X2 ( J )-FLQ4T( IA-1)*DX>«LT.C.00001) IA=IA-1 

NEAR SURFACE CF OUTER BCUNCAFY 


EZ=X 1 ( J ) 
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OPPM^ ® 

Wr quality 


C 2=FLGAT( J- 1 )^DY 
OiLL CUR {E2-,C2,a,C. 0 ) 

X NN=B + 2 . +C=«=E 2 
XT=FLOAT ( IN- 1 )’?=DX 
IF(J.EG.25) XMN=-i .0 
IF ( ABS{ Xh<N) oGE. 5. )GC TC £30 
XL=XT-CX 
X=-DY*X ^h + E2 
IF(X.LE.82) GO TO £20 

IF< J,LT.N3.CR .J.GE.N41) GO TO 512 

I F ( X. GT. X£( J+1 ) ) GO TC 542 

IF C (X3{J +1 )-XT) .GT «Q .000 1 ) GO TC £15 

7TT=T( IN-1 » J+1 )+ {Tax { 3 , J + n -T ( IN-1. J +1 ) ) < X-Xt_ )/ (X3( J+ 1 3-XL ) 

GC TO £25 
512 CONTINUE 

IF ( X, GT. X2 ( J+1 ) ) GO TC 532 
IF((X2CJ + 1 )-XT).GT.C.CCC 1) GC TC £15 
TTT=T{ I N-1 , J+1 ) + (TeX<2, J41 )-T ( IN-1, J+ 1 ) > i CX-XL )X( X2( J-f 1 )-XL ) 
GC TO £25 
515 CONTINUE 

IF ( C XL-Xl ( J+1 ) ) . GT.O ,000 1 ) GC TO 513 
IFIX.GT.XT) GC TO 533 

T7T=T{IN,J + l} + {TEX(l,J+i)-T(IN,J+l)) XT- X)/(XT-X1(J+1)) 

GC TO £25 
Sie CCNTINUE 

TTT=T ( IN,J+1J+{T(IN-1,J+1)-T(IN,J + 1))*( XT-X)/CX 
GC TO £2£ 

520 CCNTINUE 

X =DV*XWN+e2 

IF(XT.GE.X3{ J- 1) ) GC TO £24 
IF(X.GTcXT) GC TO 523 

T T7=T( I N , J-1) + (T { IN-1 , J- n-T { IN, J-1 > ) * (XT-X)/CX 
GC TO 52£ 

523 TTT=T( IN ,J-1 ) 

GG TO 

524 TTT=TEX{3 , J-1 ) 

XNL=(E2-X3(J-1))± (B2-X3 { J-1 ) >+CY+CY 
XNL=SCRT ( XNL ) 

GC TC 535 

£2£ >NL=(82-X )»(B2- X ) + DY^DY 

XNL=SORT (XNL > 

GC TG 535 
530 TTT = T( I N , J > 

X NL=XT-E2 
GC TC 535 

532 TTT=TE{IN,2) 

X NL- ( E2-XT ) (E2-XT ) + { C2-Y2{ IN ) )=f ( C2-Y2( IN ) ) 

XNL=SGRT{ XNL ) 

GC TO 535 

533 m = T(lN,j+l) 

A 

\ 
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534 XhL=(a2-XT)»{B2->T)+0^r:*D> 

>: ^L=SCRT (XNL) 

535 At-HX{ 1 , J) / >K 

T^EW= (TTT-^AL=tXNL^TGX( 1, J ) )/( 1, + XNL^AL) 
GC TO 545 

54 2 TNEW=( TB( IN, 1 )+TS( IN-1 , 1 ) ) /2 
545 CCMINUE 

TNElA=(i.-CiVE)*T3>;(l ,J)•+C^E4T^EW 
£UN=SUM+ (TNEW-TBX( 1 , J ) ) ^ ( lNEVy-T8 X ( 1 , J) ) 
TEXa ,J)=TNEVii 

IF ( J«GT, N3 ,AKD« J ,LT. N4) GC TC 55C 
GC TC 625 


C FCF J LINES INTERSECTING INNER 8GUNCART-BOUNOAR V POINTS CN THE 
C NEARER SURFACE CF INNER ECLCARY 
C 
C 

55G CCNTINLE 

IM=IFIX( X3{J)yoX+0.C0Cl) + l 

IF<AES(X3(J)-FLCAT{ INI-l)=»DX),LT.C.CCOOi) INI = INI-1 
lAIiilFI X(X4 (J) /OX+O.OOCl )+2 
E2=XE { J ) 

CALL CUR (E2 ,C2 ,B ,C, 1 ) 

>^N■=8 + 2. *C=tB2 
IF(J.EQ,N31) XMN--1.2 
IF(J.EG»2E) X^N=1,2 
XT=FLOAT{ IN I-1)*DX 
X F = XT + CX 
>L=XT-DX 

IF(J.EQ.2£> XMN=l-20 
IF ( A ES( XNN ) ,GE .5 , ) GC TO 565 
IFC (C2-{ T1 ( INIl + VK I M + 1) )/2.) .LE.CY) GCTC 570 
X = CY^=XNN+E2 
IF<X.GT.B2) C-G TC 555 
IF( X .LE .X 1 ( J- 1 ) ) GO TO 5TG 
IF { (X 1 ( J-1 )-XT ) .GT .0 .000 1 ) GO TO 552 
IF( ( XT-X) .GT. 0.0001 I GC TC 562. 

TTT=T ( INI , J-1 )+(T( INI-M ,J - 1 )-T( INI , J- 1) ) *( X-XT) XD X 
GC TC 560 

552 TTT-TBX( i,J-I) + lT(IN+l.J-i)-TeX{l,J-l))*< X-Xl { J-1 3 3 / (X F-X i (J-l 3 3 
GC TC 5&0 
555 >=-CY*x:v^ + E2 

IF(X ,L£ .X H J + l 3 ) GO TO £72 

TTT = T<IM,J-fl) + (T(INI,J-+l)-T(IN’Ii^l,J-rl)3=^'(X-XT)XCX 
GC TO 5c C ' 

562 IF ( <XL-X U J-l ) ) .GT .0 .000 n GO TO £63 

11T = TBX( 1 ,J-1 3 + (T( I NI , J-l 3-TEX (1 ,J-1 ) 3-{XT-X )/(XT-Xl{ J-l ) ) 

GC TO 560 

563 TTT^T CIM - 1 , J-l 3 -KT { IN U J- 1 )-T ( IN I- 1 , J- 1 ) ) =5‘( XT- X 3 XD X 
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OKIGINAI; PAGBlf 
Poor quauth 

5^C >^L=(B£-X) ¥(S2- X )+OYADY 

X^L=SQRT (XNL. J 
CC TC 575 

5 65 1TT=1(IM,J) 

XNL=E2-XT 
GC TG 575 
570 CONTINUE 

IF ( (XL-X 1 ( J ) ) .GT .0 .000 1 ) GO TO 574 
nT={ TB( IM , 1 ) + TEX( i , J ) ) X2 . 

XNL=E2-( XT + X 1 ( J ) )/2. 

GC TC 575 

572 TT T=T6{ INI ,2) 

X NL=<XT-E2 ) 4 (XT-E2) + { Y2( JN I )-C2 ) =♦ ( Y2( IN I )-C2) 

XNL=SCRT ( XNLJ 
GC TO £7£ 

574 TTT=TE UM , 1 ) 

XNL=( E2-XTD SI" { B2-XT) +( C2- Y1 < I M ) ) » ( C2 -Y 1 ( I M J ) 

XNL=SCRT ( XNL ) 

575 7L = HX<3 ,jJ/XK 

TNEW=( TTT + AL=sXNL*TGX(3.J) ) / ( 1 . +X NL=f: AL) 

5 85 TNEVS=C 1 .-CMS )=*T8XI 3, J )+ONE=S'TNEW 

£CN = SLN'+ ( TNE\!i-TBX (3 , J) )* (TNEW-TEX (3 , J ) ) 

TEX{ 3, J ) = 7NEVk 
CC 590 I = M ,N2 
XR = FLOAT! I-l )^DX 

IF( ( XR-X5 < J ) ) .GE .0 .000 1 3 133=1-1 

IF( ( A8S (XR-X3( J) 3 .GE. (0 .S999=?CX ) ) ) 133=133-1 

IF< ( XR-X3{ J ) ) .GE . 0. 00 Cl ) GO TO £93 
590 CONTINUE 

5S3 IF( I33.lt. M l ) I33 = N11 

TCX(3,J )=TG( 133.3) 

IF{J.6G.Nc2) TGXC3.J) =TC(N5i,3) 

C 

C FAFTFEP SURFACE OF INNER BOUNDARY 
B2=X4( J) 

CALL CUR {E2 ,C2,B,C, 1 ) 

XNN = B+2.=fC=?E 2 
IF(J.EQ.N31) XMN=-1.7 
XT=FLCAT( lA I-l T^CX 
XL=XT-CX 

IF(A8S< XI^N) .GE. E.)GC- TC 6C5 
X = -CY5i<XNN + E2 

IF( J.GT. N3. AND. J.LT.N6 ) C-C TC 614 
IF < (X-XT) .GT .0 .CCOl ) GO TO 595 
IF {X.LT.X2 ( J4l ) ) GG TC 597 
5?£ CCNTINUE 

IF < (Y2( I A I >-{ C2+DY ) ) .GT . C.OOO 1 ) GO TO 596 
T TT = TEX( 2 , J + 1 ) 

XNL=(B2-X2( j+1) ) *(B2-X2( J4 1) )+CY^sCY 
X^L=SGRT (XNL ) 
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WOR QUALITY 


G C TO 6 1 C 

59c TTT=T ( lA 1, J 41 ) 

X ^L-( £2-XT ) 4 ( E2-XT ) +DY =^DY 
>^L=SGRT( >NL) 

C-C TO 6 10 
597 CC^^I^UE 

IF{ X2( J4- 1 ) . LE . XT ) GO TO 59E 
TTT^T {I/>I»0+1)4-{T{IAI-1,J41)-T(IAI,J+1J1 =(■( xt-X)/D X 
GC TC eCG 

5 98 TTT = TEX( 2 , J 4 1 ) + ( T { I A J- 1 , J+ 1 )- TB X ( 2 , J+ 1) ) * ( X2 { J+ 1) - X ) / < X2 { J+1 )-XL > 
6CC X^L= (E2-X ) ^ ( E2-X >4-DY=?DY 
>NL=£QRT( >^L) 

GC TO 610 
eC5 m-T(I AI ,JJ 
XNL=XT-B 2 

61C /!L=HX <4 , J )XXK 

T^E\lli=( TTT-4AL=t=XNL*TGX {4,J))y(1.4-XNL=?:AL) 

611 CONTINUE 

612 TNE\»=(1 .-C^E)=FTaX^4 , J ) 4C N E =»T N EW 
SLM=SLM+{ TNEH“TBX(4 , J) ITNEU-TEX (4 , J) 1 
TEX (4 , J )=TNEW 

614 CONTINUE 

CO 6 15 I-Ic3,iN2 
XE = FLCAT (1- 1 )=?CX 

IF ( { XR-X4 ( J) ) ,GE.O .OCCl ) 14=1 

IF(AES(XR-X4{J )} .LE.G.CCCl) 14=1 4+1 
IF ( (XR-X4 ( J ) ) ,GE .0 .000 1 ) C-C TC 62C 

615 CONTINUE 

620 IF { 14 4GT .N2 1 ) I4=N21 

TGX (4 , J) =TG < 14 ,4 ) 

C 

C FAFTFEK BOUNCARY POINTS ON CUTER SURFACE 
C 

625 E2=X2(J3 

CALL CUR (E2»C2,e,C,0> 

>NN = B+2. #C!^e2 
' IFIJ.EQ.15) XMN=-2«a “ 

IF(J.E0.17) XNN=-2,3 

IF(J.EQ«2 ) XNN=6«0 

XT=FUCAT ( lA-1 )=^DX 
IF (AB S { X.VN) . G£ .5 . ) GC TC 645 
>R=XT+DX 
XL=XT-CX 
>=D Y+XN-N + BE 

J F ( J .GT .NS . AND .J .LT .Nc 1 GG TO 650 
IF(J.LE.N31 5 GO TO 635 
IF(IA.GT.N2U GO TC 627 
IF{<C2-Y4( lA )) .GT.DY ) GO TO 627 
T1T=TB( I A ,4 ) 

XNL={ C2-Y4 ( lA) 1 4! ( C2- Y4( I A ) ) 4 ( 8 2- X T ) « ( B2- X T ) 
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> ^L = S CRT <XNL > 

C-C TO €4€ 

627 CCNTINUE 

IF ( (XT-X4 ( J-1 )). GT.O .00001 ) GC TC 629 
IF ( (X4{ J- 1 )-X) .GT.O. OGCl ) GO TO 633 
IF ( {XP-X2 ( J-1 ) ) ,LT .0 .000 1 ) GO TO 628 

TTT = Tax( 4,J-1)-(TEX(4,J-1 )-TEX (2 ^ J-1 ) ) (X-X4 tJ- 13)/tX2(J-l)-X4 

1 (J-in 
GC TC 642 

626 TTT=TBX( 4 , J- 1 )-( TBX( 4 , J-1) -T ( I /S-t-1 , J-1 ) )^ (X-X4 { J-1 ) }X (X R-X4 ( J-1 ) ) 
GC TO 642 

629 CCMINLE 
IF(X.GT.XT) GC TC 64C 

630 CCMINUe 

IF { ( X-X4 C J- 1 ) ) .GT .0 .0001 ) GC TC 632 
621 TTT=T(IA,J-n 

XhL=(E2-XT)^{E2-XT) + (CY^CY ) 

>NL=SGRT( >NL) 

GC TO 646 

632 CCMIM.E 

IF( (XL-X4C J- 1) ) ,GT. 0. GCCl) GC TC 634 

633 TTT=TEX<4,J-1 ) 

>NL= ( E2-X4 ( J-1 ) ) ^ (B2-X4 ( J-1 ) ) +CY^CY 
X^L=SGRT { XNL ) 

GC TC 646 

634 TTT=T(iA-l,J-l)+(T(IA-i,J-l)-T(IA.J-l))*{ X-Xt ) /CX 
IF { (XL-X ) . C-T .0 ,0001 ) GO TO 621 

GC TC 642 

635 CONTINUE 

IF { (C2-Y 1 { lA ) ) ,GT ,DY ) GC TO 636 
T TT=Te{ I A , 1 ) 

XNL=(C2-Y1 ( lA) )5>{C2-YU JA) ) + (B2-XT)^(82-XT) 

X^L=£GRT (XNC ) 

GC TO 646 

636 CCMINUE 

IF< (XT-XK J-1 )) .GT.O. GGOGl ) GC TC 639 
IF { (XKJ-l )-X ) .GT.O.OCC 1 ) GO TC 644 
IF ( (XR-X2 ( J-1 ) ) ,LT.O .000 1 ) GC TO 623 

m=TSX( 1 , J- 1 ) - ( T8XC l.J-l)-TEX(2.J-l))=f=(X-Xi (J-1 ))/(X2(J-l >-X 1 
1 (J-l)) 

GC TC 642 

638 TTT=TBX( 1, J - i)-{ TBX( 1 , J- 1 )- T( I A+ 1 , J-1 ) ) ( X-X 1 ( J- 1 ) )y (XF-Xl (J-1 ) ) 

C-C TC 642 

639 CCMINLE 
IFIX.GT.XT) GO TO 64C 

IF((X-X1 ( J-1) ) .GT.O .0001 ) GC TC 643 
GC TO 621 

640 CONTINUE 

IF ( { XR-X2 { J-1 ) ) . LT. 0, COOl ) -GC TC 641 

TTT=T (IA,J-1) + (TBX{2,J-1)-I(IA.J-1))»CX-XT)/(X2(J-1)-XT> 
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GC TO 642 

e<ll ^T^=^{IA^J-l} + (7UA+l ,J-I)-T(IA,J-1)) « (X-XT ) / CX 

642 X ^L= (E2-X )4 (E2-X )+DY:»CY 
> ^L=SCRT ( X^L) 

C-C TO 646 

643 CCNTINUE 

IF< ( XL-X I ( J-*1 ) ) ,GT. 0. OGO 1) GC TC .63A 

644 TTT=T£X ( 1 ^ J- I ) 

XNL = (H2-X1 (J-1 ) )w(B2-Xl ( J-1 ) ) + CY4CY 
XNL=SQRT( XNL ) 

GC TC 646 

645 TTT=T(IA,J) 

X^L=£2-XT 

646 AL=HX C2 ,«) /XK 

TNEW = (TTT + AL*XNL4:TGX(2«J))/( 1 . + XNL»AL) 

647 T^EA={l.-C^^E)=TeX<2,JJ4□^E4TNEW 
£L^=£l.•M+ {TEX<2 ,J)-TNEVi, )4! <TEX(2 , J )-TNEW) 

TEXC2t J ) = TNEW 

65C CCMINLE 

DC 655 1 = 1. NXl 

CC 655 J~1,NY1 
X=f LCAT( I-l ) *DX 
TE( 1, 1)^TQX( 1 . J) 

Y = FLCAT { J-1 )4CY 

IF(ABS{X-XKJ)).LT.O.OOOC1 . AN C , A ES ( Y-Y 1 ( I ) ) ,LT «0 .0 00 0 1 ) T( I, J )=TB< 

11,1) 

I F ( 46S {X-Xl ( J ) ) .LT .0 .0000 1 . ANC«AES( Y-Y2( I ) ) .LT .0 .OCGCl ) T( I , J ) =T3( 

11,2) 

IF(AES(X*-X2<Jl).LT.0.CG0Cl.AXD.Ae£{V-Y2<I)).LT.C.CCCCl) TII, J )=TB ( 
1 I ,2) 

66E CCMINLE 

CC 657 I=N11,N21 

DC 657 J=^51,^41 
X = FLOAT{ I-l )*OX 

Y = FLCAT { J- 1 )*CY 

I F { AB£< X-X2 ( J) ) , LT. 0. OOOCl , ANC .ABS (Y-Y3 ( I ) ) .LT .0 .0 0 001 > T ( I, J ) = Xa( 

11.3) 

IF(AES(X-X4(J)).LT.O,0 00 01,AKC.AES(Y-Y4(I)) .LT .O.OGCCl) T(I,J) = TB( 

11.4) 

I F (AES (X-X3 ( J )),LT.O,OCOC1.ANC.AE£(Y-V4(I)).LT,O.CCGC1) T(I,J) =T8( 
1 I ,4) 

657 CONTINUE 
C 

C CCNFLETICN CF AN I TER A T I C N- C F E CKS FCF REGUIREC CCNVEFC-ENCE 
C 

c 

IF (ITER. EG. 1) WRITE<6,66C) 

660 FCFI^AKlhl) 

ViFITE (6,6c5 ) ITER, SUM 

665 FCRMAT( 5X , ' I TEfiA TICN= • ,I 3 ,3X ,» EFFCF=* ,E12 .5 ) 
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IF{ ITER/M£:»NTE.EG.ITER> GO TO 67C 
IF(SLN.LE.SLf-'.V) GC TC 67C 
If ( ITER ,G£ 4 IMAX) GO TO 2CCC 
I TEF=1TER + 1 
GC TO 3CC 
670 CONTINUE 

ViPIIE (6 ,700 ) 

700 FCRMATC //, 1CX-, »T-MATRI X, GI VEE ALL. INTERIOR PC I NTS , I N I T J AL SETTING 
1 OF T IS 1600 4 O* /lOX, •CNLV THOSE INSIDE BLADE ARE CHANGED AND WILL 

a ee different frcn^ i 6 CC*cc*) 

W F ITE ( 6 , 7G5 ) 

7C5 FCRMAT(y ,2X ,» J = • ,4X , • 1 • ,7 X , * 2 • , 7X , • 5 * , 7 X , • A • , 7 X , ’ 5 » , 7 X , « £ ’ , 7X , *7 
1*,7X, •a*47x,*9*,ex,*ic«,ex,*ll« ,ex,»12* ,6X,*i3* ,6X,M4’,6X»*15*/ 
22X , * I • ) 

IFtNY.GT.lS) WFITE(6,710) 

710 FCF.MAT( 8 X,« 16*,6X,»17*,€X,' 1£»,£X,*1S* ,6X,*20* ,6X,*21* ,6X,*22* , 6 X 

1 ,*23* ,6X,*24* ,6X,*25* ) 

CO 71£ i= 1 ,NX 

715 WFITE(6,720 ) I, (T ( I , J ) , J= 1 , N Y ) 

720 FORMAT (/ ,2X,I2,lX,15Fa.2 ,2(/SX,15F8 .2) ) 

WF1TE(6, 7£E) 

725 FORMAT {» 1 • ,//15X ,* GAS TEMPERATURE DISTRIBUTION ON THE BLADE INNER 
lANC OOTER SURFACE’,//) 

W F ITE ( 6 , 745 ) NX 1 
VFITE (6 ,720) ( T G ( I , 1 ) , I = 1 , NX 1 ) 

WRITE(6,7EC) NXl 

WFITE(6,72C) < T G ( 1 , 2 > , 1= 1 , NX 1 ) 

IF(M .EQ.N2 .AND4N3.EG.N4 ) GC TC 730 
WFITE (6,755) N11,N21 
VFITE{6,aOC) {T6(I,3),I-M1,N21) 

WRITE(6,7£C) M1,N21 

WF ITE (6,800 3 (TG ( I, 4 ), I-M 1 ,N2 1 ) 

73C CCNTINUE 

WFITE (c,77C > NYl 

1>FTTE (6 ,800 ) ( T GX ( 1 , J ) , J =2 , NY 1 ) 

WRITE! 6, 78 C ) NYl 

ViF ITE (6 , aOO ) ( TGX ( 2 , J ) ,J = a,NY 1 ) 

1F< M .EC. N2. AND, N3. EG, ^4 ) CC TC 735 
WFITE(6,7££) N31.N41 

WFITE (6,800) (TGX(3,J),J = N31,N41) 

WRITE(e,7?G) N31,N41 

WF ITE (6, 800 ) {TGX(4,J>,J=N31,N41) 

WFI TE ( 6 ,795 ) 

735 CCNTINUE 

W FITE (6 ,740 J 

74C FORMAT! • I* ,// ,5X,’ TEMPER ATCRE CISTRIEUTIGN GN THE ECUN CAR Y ’,//// ) 
WF ITE( 6 , 745 ) NXl 

745 F C FM AT (// ,1 OX , • ECUNDARY F C I NTS CN I-LINES FOR LOWER SURFACE OF 
1 CUTER SCUNDRY’/ ,12X, • (START FRCM 1= 2 TC I = ’,I3,» ) ’ ) 

WF ITE (6,720 ) (TE( I, 1 ) , J= 3, NX 1) 
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ViRITE(eW£C) NXl 

750 FCFM;iT(//, lOX, ‘BOLNOARN FCINTS CN I L^^ES FCR UFFEF SURFACE CF 
1 CCTER eCLNCFY’' / ,12X ,• (START FROM 1= 2 TO I=«,I3, ’ )♦) 

\»KITE(6,72C) C T8(l ,2) ,I = 1,^X1) 

If (M .EG .^2 .A^C .^3 .EG .N4 ) GO TO 7€5 
ViRITE (d,7££) M1.N21 

755 FCPMAK//, ICX, »BCLNOAR V POINTS CN I LINES FCR LCViEF SURFACE CF 
1 INNER eCUNCARY* / .1 2X , » ( ST ART FROM I=*,I3,< TO I=*,I3,» )•) 

'ARITE{e.€CC)<TB(l,3).I = Ml.N 21 ) 

VfiR ITE (6 ,760 ) Ml , N21 

76C FCRWAT(/y ,10X ,* BCUNOAFY FCINTS CN I LINES FOR UPPER SURFACE OF 
1 INNER BCUNDARY* y , 1 2X, »( START FRGN' TC I = ',I3,* )’) 

V^RITE(6,800)(TE(1,4),I=M1,N21) 

765 CCNTINLE 

WR ITE (6,770 ) NY 1 

77C FCRNAT(/y ,1 CX , • BCUNDARY FCINTS CN J LINES FOR NEARER SURFACE OF 
1 CUTER BOUNDARY' y,12X, M START FRCNi J= 2 TC J = »,I3,> )•} 

^RITE(6,eOO)(TEX{l,J),J=2,NYl) 

V«FI TE ( 6 , 7£C > NYl 

760 FORMAT (//, 1 CX , 'BOLNOART FCINTS CN J LINES FCR FARTHER SURFACE OF 
1 CCTEfi BCUNCARY* / ,1 2X , » ( ST ART FROM J= 2 TO J= « , I 3 , • )•) 

ViR ITE ( 6 , fiCC ) ( T8X( 2 , J) ,J=2,NY1) 

IF{M.£G.N2.ANC«N3,£G,N‘<|) CO TO 6C5 
ViRI TE( 6 ,7££ } N3 I . N4 1 

765 FCRMAT(/y, ICX, 'BDUNOART FCINTS CN J LINES FOR NEARER SURFACE CF 
lINNER BCUNCARY ' / 1 2 X , • ( ST A RT FROM J='.I3,» TO J= * , I 3 , • )M 
UFITE{6,8CCJ(TBX(3.J),J=r31,N41) 
inR ITE (6 ,790 ) N31.N41 

79C FORMAT! // ,10X ,* BOUNDARY FCINTS CN J LINES FOR FARTHER SURFACE OF 
lINNER BOUNDARY ' y 12X , • ( STAR T FROM J=*,I3.» TC J=»,I3.* )*) 

1ft R ITE (6 ,800 ) (TEX (4, J ), J= N3 1.N4 1 ) 

HRITE(6, 7SS) 

795 FCPMAT(lhl) 

see FORMAT (3 </ ,5X ,15F8 .2 ) ) 

6C5 CONTINUE 

IsRITE (6,610) 

610 FCRMAT( IHl/yy ,15X ,* TEMFE PATURE DISTRIEUTIGN INSIDE THE BLADE AND 0 
XN ITS SURFACE 'yyy2x, * J=* ) 

CC 885 J1=2,NY1 
DC €15 1 = 1, NX 
J=NY+1-J 1 
XR=FLCAT{ 1-1 ) ?CX 

1F( ( XR-X 1 ( J) ) .GE. C. CO C 1 ) 11 = 1 

IF (AES (XR-Xl { J ) ) .LE.C .COG 1 ) 1 1= I 4 1 

IF( { XR-X I ( J) ) .GE. C .COCl ) GC TC 820 
815 CONTINUE 
82C CC 82'5 1 = 11 ,NX 

XR = FLOAT{ I- 1 )^DX 

IF( (XR-X2{J ) ) .GE .O.COOl ) 12=1-1 
IF{.AeS(,^RTX2(. JU, (,0 .999S^,C,X U „ 12,== T2 t,1 ; 
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IF( (XP-X2(J ) ) .GE.C.COCl) GC TC 83C 
825 CCMINLE 
€3C CCMINLE 
J = NY4 1-j 1 

IF( J.GE. N4.CR. J. LE. N3 ) C-C TC 860 
835 CCNTINUE 

CC 840 1 , N2 

J=^Y+ 1-31 

XF = FL.04T < I- 1 )^>DX 

IF{ (XR-X3 ( J) J .GE.O .0001 j 133=1-1 
IF( ( A8£( >R-X3( J ) ) .GE . ( G .cggg*DX ) ) } 133=133-1 

1F< {XR-X3 ( J ) ) .GE.O .000 1 ) GO TO £45 
84C CCMIMJE 
845 CC 850 I=I 33, N2 
XR = FLCAT ( I- 1 )^CX 

IF { (XR-X4( 3) ) .GE.O.OCCl ) 14=1 

IF(4ES(XR-X4{J}) .LE.C.CCCl ) 14=14+1 

IF ( (XR-X4 (J )) ,GE ,0,000 1 1 GG TO 855 
8EC CCNTINLE 
855 CCMINUE 

3= NY+1-31 

WR ITE( 6 ,875 ) 3,TBX(1,J), ri(I,J),I = II ,133) ,TBX{3,3) 

'A RITE (6,88 0) TEX(4,J),{T(I.J),I=I4,I2),T8X(2,J) 

GC TO £65 

86G V«FITE(6,87C ) J , TBX { 1 , J ) , < T ( I , J ) , I = I 1 , 1 2) , TB X { 2 , J ) 

865 CCMINLE 

870 FCfi.MAT ( y S X , I 2 , 4 X , 1 SF £ . 2 ) 

875 FCFf/ AT (/ 2X , 12 , AX , 1SF8 .2 ) 

esc FCRVAK/l cX , • ,3X,12F3,2) 

885 CCNTINUE 

IF(£LN,Le,SLf^M) CC TC 89 0 
I1ER= ITER+1 
CC TO 300 
esc CCNTINUE 

IF{ A8S( TE^PC-TG( N5 ,3)).LE.S) GC TO 910 
C 

c 

C NEVd VALUE OF FILM CCCLANT TE^F IS OBTAINED 
C THEN GCES BACK TC GET FINAL SET OF TEMPS 
C 

TENFC=TG(N5 ,3) 

C£NSC=1 98 ,57/TG( N5 ,3) 

CFC={ S.e7-3£.l/( TGI N5 ,3 ) *•* 0 . 5 ) -2389 . /T G ( N5 ,3 ) +9 . 06 EO 5/ ( TGI I , J ) * + 
12 . ) )/28 ,S7 

VI £CC=7. 77E-07* ( TG( N5 ,3 ) ) *=»= 1 . 5/ (T G ( No ,3)4198.) 

C 

C REASSIGNING GAS TEMPS OCWNSTREAM OF SLOT TO GRIG, VALLES 

C 

C 

CC 906 I=N51,NX1 
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7G( I «2) =TGTEi^ ( I ) 

906 CONTINUE 

CC 908 J=2,N62 

TG><( 2^ J) =TG>TEN{ J) 

908 CCNTINUe 
GC TC 245 
910 CCNTINUE 

C CET4INS ISCTI-EFM/iL LINE LCC/^TIONS 
C 

940 WR ITE{ 6» S<5£ ) 

945 FCFN4T( INI ,28X, » ISQTHERM/JL L I NE LOCAT IONS •/ , SX » • I » , 4 X , * J • . 3X , « T» , 
14>,*T - I*,6X,»T - 2 • ^ e> ,• FRAC* ,/ ) 

C;sLL TEMP 
9SC CCNTINUE 

ViR ITE( 6 . SC2 ) 

992 FCPN AT { iX ♦ 1 H ) 

DC 999 NN = 1 ,NSCLT 
XRITE(6»99£) TSCLT(NN) 

9 95 FCFNATdX, ‘The ISOTHERM AT ’,15./) 

ViRI TE(e,997) 

997 FCFKAT( £X, » X-COORD • , 1 CX , « T-COCRD • ,/) 

NN2=NNY(NN) 

Vl RITE <6, 956) ( {XI£C(NN,NNL> ,YI£G(NN,NNL) ) ,NNL=1,NNZ) 

998 FCF/^ATCSX ,Fe.4,10X, Fa ,4) 

995 CCNTINUE 

CC 1000 I=N51,NX1 

TC- ( I ,2)=TC-TE(^ ( I ) 

ICCC CONTINUE 

CC 10 10 J=2,N62 

T6X (2 , J) =TGXTEf^ ( J) 

1010 CONTINUE 

CC 680 1 = 1 , NXl 

EFFd ,1)=(TG(I,1 )-TB{I,l))/(TG(I .1 )-TCCLN ) 

EFF ( 1 , 2 )-( TC( 1,2 )-T3( I , 2 ) > / ( TG ( I , 2 ) -TC CLN ) 

68C CCNTINUE 

CC ees j= 2 ,NVi 

EFF<1,J )=(TC-X{ 1, J)-TaX( 1,J) >/{ TGXC 1,J)-TC0LN) 

EFF (2 ,J) = (TGX(2,J)-T8X {2 ,J ) )/ (TGX (2 , J )-TCOLN) 

685 CCNTINUE 

Vif ITE (6,687 ) 

687 FCRMAK * 1* , // , lOX EFF IC IE NCY CF C CCL I NG *.//// ) 

WFITE(6,74£) NXl 

Vif ITS (6 ,690 ) ( EFF ( I, 1 ) , 1= 1 ,NX 1 } 

6 SC FCRNATC / ,2X.I2,1X,15F£.4,2(/5X,1£F8,4)) 

F ITE (6 , 750 JNX 1 

UPITE(6,£S0)(EFF(I,2),I=1,NX1) 

WF ITE{6, 7?G ) NYl 

V\FITE(6,e92)(EFF(l,J),J=2,NYl) 

6S2 FGFMA T( 2( / , £X , loFa,4 ) ) 

WF ITE (6 , 780 ) NYl 

ViR t 7E (6 , 69 2 ) (EFF (2,J) ,J = 2,NY1 ) 

2000 CONTINUE 
STCP 
ENC 
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SLBRCITINS C4TA 

dimension 0£LS(3 0,4J tNASC<20,4) ,C(3C,4 3 tCELM30t4) 

CIMENSiON I (30.65) ,T8(30,4) ,TBX(4,6E) , TG ( 3G , 4 ) , TGX ( 4 , 6 5 ) , 

1H(3C,^) ,HX(4,65) ,Y1(30) ,V2{3Q) ,V3(30) ,Y4(30)t)!l(65>,X2(65) , 
2X3(65) , X4( 65 ) ,BX( 200 ) .BY (2C0) tBXK 2GC) ,BYH2C0) . 1P(200) 

CIMENSICN VeLG(30,65) 

PE^L LB.MASC 

CCM-vcN T,TB,TEX,TG,TGX,H,PX.Yi,Y2,Y3,Y4,Xl,X2,X3.X4,BX,3Y,SXI.3YI 
l.IP 

CCMMGN CX.CY .NX.N'Y ,NXG ,NXI .ITEP.M ,N2,N3,K4,NP,N5.N6,N7 
CCMNGN/DAT/ NX 1 , N Y 1 , M 1 , N 2 I , N 3 1 , N4 1 . N5 1 , N52 , N6 1 , N 6 2 , M7l , NSR , 
ICELX.MASC .CELS.XK ,Ct'E,SL^N,IMAX,MS,TCCLN,LE,CP 
CCW,vCN'/TEM/ VE LG, PRKC.T EM PC, CPC, VI see, A LPHA ,FLCC, DEN SC 

READS PHYSICAL PARAMETERS CF THE GPIC SYSTEM 

PEAC(5,5) CX,CY,NXC,NXI 
5 FCPMAT(2F5.3,2I5) 

real (5,10) NXjNYtM ,N 2,N3,N4,N5,N6,N?,NP 
IG FCPMATdOIS) 

READS FPCPERTIES GF THE ELACE AND THE CQCLANT 

REAC(5,15) LB,XK,CF 
15 FCRMAT( 3F8.6) 

READ (.5, 20) PRNO, ALPHA. TCCLN, TEMPO 
2G FORMAT (4F8 o2 ) 

REAC(5,25) CENSC,CPC,VISCC,FLCC 
25 F0RMAT(2F13.10,E16,1G,F8.4) 

reads ITERATION CHECKS 

R£AC(5,30) OME.SUMM ,NTE, IMAX 
30 FCRMAK2F8.3 ,215) 

NX1=NX-1 
NY1=NY-1 
M1 = N1+1 
N21=N2-1 
N31 = N3+'l 
N4 l=N4-i 
N51=N5+1 
N52=N5+2 
N61=N6-1 
N62=N6-2 
N71=N7~1 
NSP=N5-1 

INIATIALLY TEMPS AT ALL FCINTS APE SET TO A CCNVENIENT VALUE 
DC 21 1=1, NX 
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CC 31 

J = 1 ,NY 

31 

T( I , J) = 

2 400. 


DC 32 

1 = 1 , NX 


CC 32 

J = l,4 

32 

TEC I , J3 

= 2400. 


DC 33 

J = 1 ,NV 


CC 33 

1 = 1,4 

■33 

TEX ( I ,J) 

= 2400 


READS IN BGUNCARY PCIMS SET ANC THE CGRR E S PC \C I NG SECUENCE NUMBERS 
fCF THE 8 CUNCAFY PCINTS 


BEAC(5,35) (Yi( I) ,I=1,NX1) 

35 FCFNAT( { 12F6.3} ) 

REAC(5,40) (IP( n ,I = 1,NXI) 

^0 FCPMAT(20I4) 

CC 45 1=1, NXl 
X = FLCAT( I-l )=^CX 
ex( IPU ) ) = X 
45 eY{iPm)=Yi(i) 

REAC(5,35) (Y2{ U-,I = 1,NX1) 
PEAC{S,40) (IP{ I),I=1,NX1} 

DC 50 I=1,NX1 
X=FLCAT{i-i)^CX 

exciPd) )=x 
5C BY(IP(n)=Y2(I) 

PEAC(5,35) (Xl{ J ) ,J = 1,NY1> 
REAC{5,40) (IP{ J),J=1,NY1) 
CC 55 J=1,NY1 

Y = FLCATU-1 )*CY 
8X{ IPlJI ) = XH J) 

55 EYUP{J))=Y 

R5AC{5,35) CX2( J) ,J = 1 ,NY1) 
R£AC{5,40) (IP{J ) ,J=1,NY1) 

CC 60 J=1,NY1 

Y = FLCAT( J-1)=^CV 
EX ( IP{J) )=X2(J ) 

6C B V{ IP{ J) ) =Y 

PEAC<5,35) (Y3(I),I = Mi,N21) 
PEAC(5,40) ( IF( I ),T = Ml,N 2 l} 
CC 65 I=Mi,N2l 
X = FLCAT( i-1 )-CX 
exi UP{ I ) ) = X 
65 6 YI (IPU ) )=Y3 { I ) 

R£AC(5,35) { Y4( I) ,I=M1 ,N21) 
PEAC(5,4D) {IP{n,I = Ml,N2D 
DC 70 I=Mi,N21 
X=FLCAT{ I-1)*CX 
e X I ( I F ( I ) ) = X 
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70 8>I C IP( I) ). = Y4(I ) 

PE.£C(5,35) (X3( J) ,J = N31,N41) 
PEAC(5,40J ( IP{ J) ,J=^31,^41) 
CG 75 J = K-31,N41 
Y = FLC;STU-U=i‘CY 
8X1 { IP( JJ )=X3( J) 

75 EV! (IP{J) )= Y 

RE-4C( 5,35 3 { X4 ( J ) , J= ^3 1 , N’4 I ) 
RE/'C(5,40) ( IP( J) tJ = N31,N41) 
CC 77 J=^2^,^41 
Y^FICAK J-l)#DY 
EXIl IP(J3 )=X4(J ) 

77 S VI { IF( J) )=Y 


QUALITY 


RE^CS IN GAS TENFS ANC TFEN FEAT TRANSFER CCEFFICIENTS 


PEAC(5,90) (TG(ItlJ,I=l,KXl) 
REAC(5t90) [TGI I ,2) ,1^1 , NX 1) 
REAC(5,9G) (TG( I ,3 ) , I = N11,N2I) 
REAC [5,90)ITG( I ,4 ) , I = M1,N21) 
90 FCPPAT{10F8,2 } 

REAC (5,90 ) (TGX( 1, J) , J^2 ,NY1) 
R£AC(5,90) (TGX(2,J) ,J=2,NYU 
R£AC(5,90)(TGX(3,J3 ,J=^31,^41) 
REAC(5,90) [TGX(4,J) ,J=^31,^41J 
PEiC (5,993 (F( I, 1), 1=1, NX!) 

REAC (5,99) (H( I ,2) ,I = 1,NX1) 
REAC(5,S9)(F(I,33,I = Nll,i\;21) 
REAC(5,99)(F(I,4i,I = Ml,N21) 

99 FCRMAT(10F8.4) 

REAC (5,993 (FX(1,J),J=2,NY1> 
REAC( 5,99 3 (l-X(2,J) ,J=2,NY13 
REAC (5,99 3 (FX(3 ,J3 , J = N31 ,N41) 
peAD(5,S9 3 (HX(4,J3 ,J = N31 tN413 

REACS YELCCITIES 


N'R= 1 

^ IF(NP.6G.13 NP=2 
PEAD( 5,1053 (VELG( I ,NP) ,I = N51 ,NX1J 
PEAC( 5 , 10 5 3 (VELG(NF,J ) ,J = 2,N6) 

ICE FORYAK 10F8.2) 


REACS CCCLANT NASS FLCV^ 

REAC(5,11C) (MASC{ I,33,I=N11,N5U 
REAC(5,110) (NASC( 1,4), T = M1,N513 
PEAC (5, 110) (DELN( r ,3 3 , I = M1,N5 ) 
PEAC(5,110)(EELN(I,43,I = M1,N53 
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.OF POOR QUAIiirY' 


C REiCS SURFACE CIST^SNCES 
C 

RE/SC(5,110) (CELS(I»3 3, I = ^11,^5J 
RE4C(5,UO] (DELS (I ,4),I = M1,N5) 
110 FCF.v AT (10F8.4) 


eCFQ CHECK FGR IIMPbT CATA 
I TER=1 

WPITE(6, 1151 

113 FCRPiT (// ,23Xt • INPUT CATAM 
WRITE! 6 ,120) 

120 FORMAT!///, 15X, '8LACE BCUNCARY PCINTS CGCRCINATSS '/// //17X CUTER 
1 eCLNCARY* ,9X, ' INNER 8CUNCARY ',//10X,'PT,NQ.*,5X,’;<',EX,*Y*,6X, 
2'PT.N'0.» ,6X,»X» ,9X, 'V'/) 

NXI1=NX!+1 

URITE(6,125}((I,6X(I),BY{I},I,SXUI),BYI{I3),I=1,NXI) 

WPITE(6, 130) { (I ,eX(I ) ,8Y(I )) ,1 = NX! 1 ,NXC) 

125 FCRMAT(9X,I6,2F10,6, I6,2F10,6) 

130 F0Rf^AT(9X, I6,2F10-6) 

WRITE (6 ,135 ) 

135 FCRMAT( »!• ,//15X, ’INPUT DATA » //15X GAS TENFERATURE DISTRIEUTIGN 
ICN THE PLACE INNER ANC CUTER SURFACE’,//) 

WRITS(6,7A5) M,NX1 
WRITE(6,720) I TG { I , 1 ) , I =1 ,N X 1 ) • 

WRITE(6,750 ) M,NX1 
InFITE(6,720) (TG( !,2),I = 1,NX1) 

WRITE (6,7553 M1,N21 
WRITE (6, 8 00) (TGU ,3 ) , I = M1 ,N21) 

WFITE(6,760J M1,N21 
WPITS(6,8C0) (TGd ,4) ,I=Ml,N2I) 

WRITE (6,770) NYl 
WRITE(6,800) (TGX( l,J) ,J = 2 ,NY1 ) 

WRITE(6,7S0) NYl 

WRIT£(6,800) (TGX(2, J) ,J = 2,NYl) 

WRITE(6,785) N31,N41 
WPITE{6,800) (TGX(3,J),J=A3i,N41) 

WRITE(6,7CC) N31,N41 
WRITE(6,800) (TGX14,J),J=N31,N41) 

WPITE(6,755J 
WRITE(6, 150 ) 

150 FCRMAT( • 1 •, //15X, ' INPUT DATA • // 1 5 X , ’ HEAT CGEFfICIENT DISTRIBUTION 
ICN THE FLACE INNER ANC CUTER SURFACE*,//) 

WRITE(6,745 ) M,NX1 
WR!TE(6,720) ( h ( I , 1 ) , I =1 , N XI 3 
WRITE(6,750 ) N1,NX1 
WRITE(6,720) ( H { I , 2 ) , I = 1 , NX 1 ) 

WRITE (6,755) Nli,N21 
WRITE(6,300) { H( I,3.),I=M1,N21) 
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WRITE{ 6,760) Nli,K21 
WRITE(6,8003 I 1- ( 1 , 4 ) , I=N 1 1 , i\2 1 ) 

WRITE [6,770 

WPITE{6,800J ( hX{ 1 , J ) t J = 2 ,K'rl ) 

WRITE(6,780) NYl 

WPITE(6,8G0) ( l-X(2, J) ,J = 2,NY1) 

WRnE(6,785J 

WRITE{6,8Q0) IHX(2,J) ,J = N31,N4i) 

WRITE(6,790) N31,N4L 
WRITE{6,8C0) ( hX(A,J), J = ^31,■^'4n 
WR!TE{6,7<;5) 

WRITE{6,165) 

165 FQRy^TC « 1 ‘ , //15X, » INPUT T A * / / 1 5X , ' COOL ANT MASS FLOW DISTRIBUTION 
1' ,//J 

WPITE(6,755 ) M1,N51 

WRITE(6,17C) (MASC( I ,2) ,I = M1,N5U 
WPITE{6,76CnNl 1,N51 
WRIie(6,170) (MASCn ,A) ,I = M1,N5U 
170 FCPMAT(3(/,5X,15Ffi,4) ) 

WRIT£(6 ,175) 

175 PQRMAT( *1 ’ ,//15X, ' INPUT DA TA • , / / 1 5 X^ • VELC Cl T I E S AT POINTS WHERE A 
ICIAEATIC WALL TEMP IS TC EE C ET E RM IN E C » , / / ) 

IF(NP-EQ.l) GC TC 178 
WKITE(6,7A5J N51.NX1 
KPITE(6,180) (VELGd ,NR) ,I = N5i,NXl ) 

WRITE{6,770) K6 

WPITE (6,180 } (VELG(NR , J) , J = 2,N6) 

GC TC 182 
17E CCNTINLE 

WRITE(6,750) N51,NJa 

WRITE (6, 180) CVELG( ItNR) ,I = N51,NX1) 

WRITE(6,7eO) N6 

WRITE (6, 180) (VELG(NP , J)^J = 2,N6) 

18C F0RMAT{3(y ,5X,15FS,2)) 

162 CCMINUE 

WRITE (6,185) 

165 FORMAT! » 1» ,//15X , ‘ INPUT DATA') 

WRITE(6,1SG) 

ISO FGPM.AT(//,15X PHYSICAL PARAMETERS CF THE GRIC SYSTEM*) 
WPIT£(6,1S5) CX,DY,NXC,NXI 

1G5 FORMAT!// ,15X,'DX=»,F5.3,5X, 'DY= ' ,F5.3, 5X, «NXQ=' ,I5,5X,*NXI=*,I5) 
WRITE (6,200) NX,NY,M,N2,N2,N4,N5,N6,N7,NP 
2 00 FCRMAT(//,15X,' NX = * , 15 ,5X, *NY^ ' , I 5 , 5X , ' N 1= * , I 5 , 5 X , » M 2 = ' , 1 5 , // 15 X , 
l'N3 = » ,I5,5X,5X, »NA=' ,I5,5X, *N5=',I5,5X, *K6=' ,I5,5X, 'N7 = ' ,15, 
25X,'NF=' , 15) 

WRITE(6,2C5) 

205 FCRMAT(// , 15X, ' PROPERTIES OF THE ELACE AND THE COOLANT') 
WPITE(6,210) L8,XK,CF 

210 FORMAT!// 15 X, 'LE=',F8,6,5X,»XK=»,F8.6,5X,*CP=',F8.6) 

WPITE(6,215) PRNC , ALPHA ,TCCLN , TEMPO 
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215 FGRMiiK//, 15X, « PRNQ= ‘ tFfi.2 ,5X, *ALFHA = » ,F8 . 2 , 5X , » TCCLN = ’ ,F8„2,5X, 
1«7EVFC = SF8.2) 

V,RIT£(6,22Q) DEN SC ,C FC , VI SCC ,FLCC 

220 FGRiVfiT(//, 15Xt ‘CENSC^* ,F13.L0t5X* »CFC= • ♦ F 12* 10 , 5 X , • V I SCC= • ,E16.10 
1,5X , 'FLCC=» ,F8.4) 
i^PITE(6»225) 

225 FCRNAT(//15X, ’ITEP^TICN CFECKSM 
Vi RITE (6,230) CV E , SUP N , N T E , X 

230 FCP''^T(//15X,’0i'^E=' ,F8-3-,5X,»SUPP=*,F8.3,5X,'.NTE = > ,15 , 5X, 

1 * !P4X = ' , 15 ) 

PFIN’T TITLE FGR SUBRCUTINE CURVE OUTPUT 
VPITE(6,235) 

235 FCPPAT( LHl, //,5X, 'SUBROUTINE CURVE -OUTPUT- ' // 5 X , ’ I =0 , CUTE R 3GU-N 
ICART* 1 = 1, INNER B CUN'CA RY '/ 5X , » I F BCUNCARY PCIKT NUMBER') 
^^PITE{6,2^0i 

240 FGPMAK / ,5X, ' { E2,C2 ) IS THE PCINT ON THE BCUNCARY, (61, Cl) AN 

1C (B3,C3) APE THE SURRCUNCING PCINTS CN ELACH. A,B,C, ARE THE COFF 
2S*/5X ,*0F THE FAREGIIC CURVE Y^A-s-E*X+C*X =^X ' , // , 4X , ‘ I » , 3X , • I P ' , 5X , 
2'El',aX,'Cl*,8X,'B2»,8X,'C2',8X,'e3',ax,*C3',5X,*A',9X,'B',9X,*C*, 
h/n 

244 CCNTINUE 
IF(NF,EQ.l) CD TO 245 
N5 = i\5-1 

N5i=N5+l 

N£2=N5+2 

245 CONTINUE 

720 FCRNAT{/,2X,I2,lX,i5Fe.2,2(/5X,l5F8*2) ) 

745 FCRPAT(//,iOX, ' eOUNCARY POINTS ON I-LINES FGP LOWER SURFACE CF 
1 CUTER BCUtNCARY'/,12X, '(START FPCW I=»,I3,* TQ I = »,!3,')M 
750 FQRNAT{//, lOX, 'BOUNDARY PCINTS ON I LINES FOP UPPEP SURFACE CF 
i CLTEP ECUNCAPY»/,12X, '{START FPCR I = ',I3,' TQ I=',I3,»)M 
755 FCRNAT(y/, lOX, 'BCUNCAPV PCINTS CN I LINES FCR LCWEP SURFACE OF 
1 INNER eOUNiOARY »/, 12X, '( START FROM I = ‘,I3,» TO I = *,I3,» )’) 

760 FORMAT(//,lOX,'eCUNCAPY PCINTS CN ! LINES FOR UPPER SURFACE- OF 
1 INNER eOUNCARY»/,12X, '(START FRCP I=»,I3,' TO I=',I3,' )') 

770 FCPiMAT(//,10X,* BCUNCAPY POINTS CN J LINES FOR NEARER SURFACE OF 
1 CUTER BOUNDARY'/, 12X, ’(START FRCP J= 2 TO J=»,I3,* )') 

?£0 FOPyAT(// ,10X,' BCUNCAPY PCINTS CN J LINES FGR FARTHER SURFACE OF 
1 CUTER E0LMCARY»/,12X,« (START FRCN J= 2 TO J=»,I3,' )M 
785 FCRNAT{//,1QX,' BCUNCAPY PCINTS CN J LINES FOR NEARER SUR.FACE CF 
lINNEP BCUNCARY ' /I 2X , ' ( START FROM J=',I3,' TG J=',I3,' )•) 

790 FCRMAT(//,inX,' BCUNCAPY PCINTS CN J LINES FOR FARTHER SURFACE OF 
lINNEP BQUNCARY * / 12X , ' ( STA PT FRCP J=',I3,’ TC J=«,I3,' )') 

795 FCPNATdHlJ 

SOO FCR-'^AT (3{/, 5X ,15F8*2 ) ) 

FE7LPN 

END 
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SL8PCUTINS ^'ESi- (I ,JTT^EK) 

CCNKGN TOO ,65) ,TE (30,4) ,TBX {4,65) ,TG( 30,4) ,TGX(^, 65) , 
1H(3G,4),HX(4,65),Y1(30),Y2{3G),V3{3Q),V4(30),X1(65),X2{65), 
2X3(65) ,X4 (65) ,8X (200) ,8Y(20Q) ,exH20G) ,BV! (200 ), IP (200) 
CCf'NGN CX,CY,NX,NY ,^XC ,NXI ,ITER,M,N2,N3,N4,NP,N5,N6,N7 


SUBROUTINE FINDS SI-1 , S 1-2 ,DcL-U , DEL-2 . FOR EACH IKTERICR PCIKT AND 
GETS NEW TEVPS FOR EVERY INTERICR FCIM 


N5l=N5+l 

PAR = FLCAT(I-l}=i=CX 

IF(NP.SG.l) GC TG 7 

IF{ J.GT.N3 .ANC. J.LE.A6 ) GC TC 6 

GC TG 7 

6 IF( I.GT.N51,ANC.I.LE.N2) GC TO 12 

7 CCNTIME 

IF((PAR-X1{J)).LE.CX)S1=(FAP-X1(J))7CX 
IF( (PAR-Xl (J) l.LE.CX )T1=TEX( 1,J) 

IF ( (PAR-XK J) ) .LE.CX) GC TC 20 
IF( J.GT.N3 ,ANC,J,LT.N4)GG TC 1C 
GC TC 15 

10 IF(PAR,GE.X4{ J) )GC TC 12 
GG TC 15 

12 IF( (PAP-X4( J) ).L£-CX)S1=(FAR-X4{J) )/CX 

IF( (PAR~X4(J) ) .LE.CX)T1=TBX(4, J) 

IF { (PAR-X4I J)) -LE.CX) GC TC 20 
15 Si=1.0 

T1=T ( I-l , J ) 

IF(ABS(X1( J)-PAR + DX) .LT,0.0C001) Tl = TeX( 1 ,JJ 
IF(AeS(X4U)-PAR4CX ) . LT » C. COCO i ) T 1=T8X(4, J ) 
2G CCNTINUE 

IF((X2(J)-PAR).LE.CX)S2={X2(J)-FAP)/CX 
IF{{X2(J)-FAP),LE.CX) T2=TEX(2,J} 

IF( ( X2( J) -PAP), LEoCX) GC TC 35 
IF( J,GT.N3 .ANC.J.LT,N4)GG TC25 
GC TC 30 

25 IF ( X3( J ),G7,PAP )GC TC 27 

GG TC 30 

27 IF( (X3( J)-PAR).LE.CX)S2=(X3(J)-FAP)/CX 

IF( (X3 !J )-PAR) .L£.CX)T3=TBX(3, J ) 

IF( (X3{ J)-PAP ) ,LE,CXJ GG TC 35 
30 S2=l. 

T3=T(I+1,J) 

IF(ABS(X2( J)-PAR-CX) ,LT,Q.G0001) T3=TBX(2,J) 
IF (/IBS (X3 ( J )-PAP-CX ) ,LT,C,OOGCi) T3 = TBX( 3,J) 
35 PAR = FLCAT( J-1)=^'DV 

IFINP.EQ.U GC TC 3 7 

IF( J,GT,N3 .ANC, J,lE,h6) GC TC 36 
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GC TG 37 

6 IF{ I.GT,N51.AND.I .LT .N2) GC TC A2 
1 CCNTIKLE 

IF{ ( n ) ,L£.CV)D1={PAR-Yin ) J /CV 

IF ( (PAR-VU I ) ) .LE.CY)T2=Te(I , 1) 

IF( (PAR-YKI) ) .Lt-CY) GC TC 50 
IF{I.GT.M.AiNC.ULT,N2)GG TC^O 
GC TC 45 

40 IF (PAR.GE. Y4( inCC TC^2 
GC TC 45 

4 2 IF ( (PAR-Y4( I ) ) ,LE.CY)D1=(PAR-Y4a ) )/DY 
IF{ (PAR-Y4 ( I ) ) ,LE.CY3T2=TB{ 1,4) 

IF { (PA.R-Y4( U ) .LE .CY ) GC TC 50 
45 Ci=1.0 

T2=T{!,J-1) 

IF (ABS( Yin )-PAR-i-DY) .LT.O.OOCOl ) T2=TS ( I ,1 ) 
IFtABS{Y4( I )-FAR4CY l.LT.O.OOOO n T2 = TS(! ,4) 
50 CCNTINLS 

IF( (Y2( I )-PARJ ,LE.DY)D2= (Y2(I )-PAR) /CY 
IF{(Y2m-FAF) ,LE=CV }T4=TE( I, 2 ) 
IF({Y2(I)-FAR).LE.CY) GG TC 65 
IFn.GT.M.ANC.I.LT .K2)GC TC55 
GC TC 60 

55 IF{ Y3( I ) -GT.PARIGC TC57 
GC TC 60 

57 IF((Y3m-PAP).LE,CY)C2=(Y3(I)-FAF)/CY 
IF n Y3 ( I )-FAP ) .LE,CY }T4=T5( I, 3 ) 

IF ( (Y3 ( I )-FAB) ,LEoCY ) GC TC 65 
6G C2=l-0 

T4=T{ I , J+1 ) 

Ir{A8S<Y2(! }-PAR-CY) .LT. 0.00001) T4=T8(I ,2) 
IF{A8S(Y3{ I}"PAR-DY).LT,0.QC<}01) T4*TSn ,3) 
65 CCMINLE 

Al=Tl/Si/(Sl+SZ) 

A2=T2/C1/ (C1+C2} 

A3=T3/S2/ (S1+S2) 

A4 = T4/C2/ < C14C2 3 

e=i./Sl/S2+(CX/CY)4 (CX/CYI/C1/G2 
TN£W=( A1 + A34CX4CX/(CY4DY)=^ (A24A4) ) /£ 

RETlP^ 

ENC 
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SUECCUTIN? SLCtXl,Yl»X2,Y2,X3TY2,;^TE,C) 


FINES EQUATICN TFFO TFPEE FGIMS AS Y=A + B=^X+C=i=X>>X 


C4LL ENCFL'rt 

IF<A8S(X1-X2).LE.0.Q01-AKD.AES{X2-X3).LE. 0-001) GC TC 10 

XS3=X3*X3 

>S2=X2^X2 

XS1 = X 1^X1 

C1 = X2X=XS3-X3^XS2 

C2=X1*XS3-X3«XS1 

C2=X14XS2-X2^XS 1 

C=C1-C2+C3 

IF (A85(C) .LE-0,CCOQ01) GC TC IQ 
A1 = Y1*(X2X=XS3-X3=)'XS2 ) 

A2=Y2’^'(X1^XS3-X3X'XS1} 

A3 = Y3=^{Xl-i=XS2-X2=^XS U 
A= (A1-A2 + A3 )/C 
e l = Y2X'XS3-Y3X‘XS2 
B2=Y1*XS3-Y3^XS 1 
B3=Y1*XS2-Y2*XS1 
B = (B l-82+83)/C 
C 1=X2*Y3-X3*Y2 
C2=X1X‘Y3-X3*Y1 
C3 = Xl-'iY2-X2X<Yl 
C= {Cl-C2^■C3 )/C 
RETLRN 
10 A=1G0. 

E=1CQ. 

C=10Q. 

RETURN 

END 
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SUBPCUTINE CUP (EH ,C2 ,B,C, n 

CC.’«'i^GN TOO ,65 J ,TB( 30,4) ,TBX (4,65) ,TG £ 3 0 ,4 ) , TGX ( 4 , 65 ) , ‘ 
Ih £30,4) ,hX(4,65 ),Y1(3G) ,Y2 0G) ,Y3( 30) , V4£ 30) ,>1£65) ,X2{65) , 
2X3 (65 ) ,X4 (65) ,9X(2G0 ) , BY ( 200 ) , EX 1 ( 2CC J , B V I £200) , IP (200) 

CO WON CXrCYfNX ,NY,NXC,NXI tITFP,M,NZ,N3,N^,NP,N5,N6,N7 


SUEFCUTINE TC FINC ;fiCJ/!CENT TWO POINTS FQP THE FOINT(82,C2) iND THEN 
PASSES A CURVE Y=A+EX+CX2 ANC CETERPINES A,E,C. 8,C CCEFFICIENTS 
APE RETURNED TO MAIN PPGGRAP TC GET SLOPE CF CURVED BCUNDARY AT THE 
POINT 


IFd.EC.DC-C TC 25 
DC 10 P=1,NXC 

lF(ABS£P2-eX(N} ).Lc.C.CG01.AND,ABS(C2-BY(Y) )»Le.0»CG01) GG TO 12 
1C CCNTIME 

12 IF(V,EC.l) GC TC 15 
E1=BX(H-1) 

Cl=BY(N-i) 

GC TO 13 

15 ei=ex£Nxo) 

C1=BY£NX0} 

13 CCNTINUH 

iFcy.ec.Nxc) gc tc loo 

63 = eX(M+l) 

C3 = PY(M+1) 

100 CONTINUE 

IF(P«EC«NXC) E3=EX(1J 
rF(H.eC.NXC) C3=BY(1) 

IF £ B3,£C.B2*ANC.C3,E£.C2 )GC TG lA 
GC TG 37 
lA B3 = aX(N+2) 

C3 = eY(Y42) 

GC TC 37 

25 DC 30 H = ItNXI 

IF£ ABS (62-eXI £M )) «LE.G,OG01.AND.ABS (C2-BYI (f') ).Ue. 0.0001 ) GG TC 32 
30 CCNTIME 

32 IF£P.EC.l) GO TG 3A 

ei = exi £Y -1 ) 

C1=3YI (Y-1) 

GC TC 33 

3A E1=EXI(NXI) 

Cl^eVUKXI ) 

33 CCNTINUE 
IF{Y.EG.NXI) GC TC 200 
83=eXI(N + l ) 

C3=BYI (M + 1 ) 
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200 CCMINUe 

IF .EC.N;<1 ) E3=EXI(1) 

IFtf^.EC.NXI ) C3 = SYI(1) 

IF(e3.EG«E2.^KC-C3.EG-C2)GG TG 35 
GC 7C 37 
3 5 B2 = 6 XI(.w + 2) 

C3 = eYHf^ + 2) 

37 CALL SLC(Bl,Cl,E2,C2,E3,C3,Ar8,a 
Xf/N = E + 2-*C^e2 

IF ( ITER. SC. I ) WRITE{&r40)I»iV,Bl,ClfE2,C2»e5,C2TA,B,C,XMN 
50 CCNTINUE 

40 FCFf^AH 3X , I2» 15,9F10.5,3X» *Xf^N=’ ,F10.63 
RETURN 
ENC 
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SLEPCLTINE TEF4 C ( I , J , T^C ) 

CII^ENSICN VELG(30t£5} 

CCN^CN T{30,65),TEO0,4},TBX(^,e5),TG(30,^),T6X(4f65), 
1H(3Q,A) »HX(^,65),V1(3G),Y2(30J ,V3(30) ,Y4(3<1) ,XU65) ,X2(65) , 
2X3 (65 ] f X4{65),BX{200J taY(200),EXI (200i ,eYU2QG) , IP (200) 
CGMYCN OX, CY,NX ,NY,KXC,NXI ,I TER ,M ,N2,N3 ,N4,NF ,N5 ,N6, N7 
CCi''f'C.N/T5f'/ VcLG jFRNOvTENFC , C F C , V I S C C , 4 L P‘r a , F LOG , CEN SC 


SLERCL’TINE TC FIND 4CI4E4TI.C ’rt4LL TENFEFATURE 

IN THE PAIN FRCGRAP GAS TEMPERATURE AT POINTS DOWNSTREAM FROM THE SLG 
E£ REPLACED BY CCRR ES FC NC I NG ACIAEATIC WALL TEPFEPATURE 


N51=N5+1 

REYNC=FLQC/VISCC 
I F(NF,EC.l ) GC TC 6 
IF(I.EC.l) GC TC 10 
XRE = FLGAT (N51-1)X>CX 
XF = FLCAT( I-l )=^«DX 
!F(I.GE.N2J GC TC 5 

CIST=:( ( {Y4(N51}-YA( I ) )4>?2 + ( XRE-XP )7X^C.5)/12. 

GC TC S 

5 CCMINUE 

CIST={ ({Y4(N51}-Yl(I))4=i'2 +(XRS-XP)4*2 )X-^Q.5)/12, 

GC TC 9 

6 CCNTINUE 

IF{ I.Ee.l.CP.I,EG,2) GO TC 10 
XPE = FLQAT(N5 -U=^DX 
XF=FLCAT( !-l)=i'CX 
IF(T.GE.N2) GC TC 8 

DI5T=( ( (V3<N5 )-Y3(in^=^2 + (XF-XPE)=^‘’»2)^=^0.S)/12. 

GC TC 9 

8 DI£T=(UY3(N5 ) - Y2 (I 2+ { XF-XR E ) =>72 ] , 5 ) / 12 . 

9 CCNTINUE 

CPG=(9.87-'3 6.1/ (T6( I , J ) 7*0 , 5 ) - 2389 . /TG ( I , J ) +9 .06 BOS / ( TG ( I,J )7A2. ) ) 
1/28.97 

VISCG = VISCC7 ( (TG { I , J )/TEPFC)7>s=i. 5* (TEMPC+19S. )/(TG(I,J)+198.)) 
CE.\SG = CeNSC7TEMFC/lG ( I , J ) 

GC TC 20 
1C CCNTINUE 

IF(NP.EC.l) GO TO 16 
XR£ = FLCAT(N51-1 )7CX 
YP = FLCAT( J-1)=*CY 
IF(J-Lt.N3) GC TC 15 

CIST=( ( (YA(N51)-YF)=5'*27(XRE-X4( J) 3 *=*2 J 7=«=C .5 ) / 12 . 

GC TC 19 
15 CCNTINLE 
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DIST=( ( (Y^ (N51 )-VF)^^2 + <XPE-Xl{ J) ) «*2 ) 0 . 5 ) / 12 . 

GC TC 19 
le CCNTINie 

XRE=FLQAT(N5 -1)*DX 
VP = FLC/1T ( 

IF(J,LE.K2) GC TC 18 

CIST=( ( ( Y3 (N5 3-YF )>!-*2 + (X3( J)-XFE )^=*2J**G.5) /12. 

GC TC 19 

18 CIST=( ( { Y3(N'5 ) -Y P > *^^2+ ( X2 ( J ) -XF E ) -*^2 3 .5 ) / 12 . 

19 CCNTIKUE 

CPG = {9. 87-36. 1/ (TGX { I , J ) «=>0.5 3 -2389 ./TGX ( I , J 3 +9.06E05/ (TGX tl , J ) =^=i'2 
1.3 3/2S.97 

VISCG=VISCC=«=( (TGXn, J )/TEf<FC)>i'’!'1.5* (TEf^PC+ 19 8. 3 / ( TG X (I , J 3 +198 . 3 3 
CENSG=CENSC^TEMFC/TGX< I , J 3 
2G CIYF = CENSC-=«=VELG C I ,J )=tClST/FLCC 

e£TJ=l.+5S-C4=i=R£VNC’f-VI SCC« S I N ( A L PH A ) / V I SCG 

D^C^' = 1,+0.32S*{ REYKC 3*^ {-0 .2 3>f CPG/CPC=^ (VISCG/V ISCC ) 2=!“ { 0 1 M P 3 -=* ( 

10.83*eETA 

CNLNc = 1.9=?F PNC* ^0.66 666 
FCEFF=DNGf^E/CNC^ 

IF{J.EC.1.CP.J.EC.23 GC TC 3C 
TAC^TGX (I ,J 3 + FCEFF* {TEN-FC-TGX( I, J3 3 
GC TC 40 

30 TAC=TG( I, J3+FCEFF*(TEMPC-TG( I7J) 3 
40 CCN TINGE 

IF(TAC.LT.TEYPC 3 TAC=TFj^PC 
50 CCNTINLE 
PETLPN 
ENC 
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SLEROUTINE TEMP 

CCf-i«Ch T(30»6 5J,Ta(30»4),TBX(4.65)^TG(50,41»TGX(4,6S)» 

1H(3G.4) fHX{ 4,65) ,Vi{30) ,V2(30) ,Y3(30) ,Y4(30),X1<65),X2<65), 
2X3 (65 ),X4(65), ex (£0 0),eY (200,8X1(200), 8 YT(2CQ),IP<2C0) 

CCNMQN DX,OY , NX , NY , NX C , NX I , IT ER , M , N2 ♦ N3 , N4 , NP, iM5, No , N 7 
C 

C 4LCNC I-LINE NIN TEMP STCREC IN • IL • AND MAX TEMP STORED IN «IH* 
C 

N X 1=NX- 1 
CC 987 1=2, NXl 

IF( I .GT, M ./ShO.I <.LT,N2 } GC TC 967 
i;s=IFIX(Yl{ n/DY+O.OOGl) + l 
IC = IFIX(Y2( I)/CY40.000 1 )+2 
I T1=TB( 1,1) 

IT2=TE( 1,2 ) 

IL=IT1 

IF=IT2 

IF ( ITl .GE . IT2 ) IL=IT2 
IF (IT l.GE. IT2) IH=IT1 
CC 955 J= 14 , ID 

IF ( 7 ( I , J) ,LT .FLOAT! -IL ) J IL=IF IX ( T ( I , J ) ) 

IF( T( I , J ) ,G7 ,FLOAT( IH))IF=IFIX(T(1,J)) 

950 CONTINUE 
S55 CCNTINLE 

CC 965 II=IL,IH 
CC 962 J=IA,IC 

IF((FLGAT(» )=i'DY).GT.V2 (I ) ) T ( I , J+1 )=T8 ( 1 , 2 ) 

IF (Y 1 ( I ) .GT , (FLO AT( J- 1 } ) T( I , J)=TB( 1,1) 

T \=FLC AT (II) 

I F ( T( I , J) ,LT .TV. AND. T { J , 4+ 1) , GE, TV) GC TC 960 
IF (T ( I, J ) , GE .TV .AND ,T ( I, 1 ) ,LT. TV ) GO TO 960 
GC TC 962 
960 CONTINUE 

FX=(T(I,J)-TV)/(T(I,J)-T(I,J + 1}) 

CALL I SOTh( I I ,I , J,DX,DT , FX ) 

WRITE (6, 972) I,J, II, T{ I,J),T(I,J+1),RX 
962 CONTINUE 
965 CONTINUE 
GC TC 987 
967 CCNTINUE 

IA= IF IX{ Y 1 ( I )/DY + C. OCC 1 ) + 1 
I E= IFIX( Y3 ( 1 )XDY +C .000 1 ) 42 
IC^IFIXIVA! I )/DY+C. 0001)41 
I C= IF IX { Y2 ( I )/DY4C.0 0 0 1)42 
I 71 = TE( I ,1 ) 

IT £=Ta( 1,5) 

IT3=TE ( 1,4) 

I T4=TB( 1,2) 

IL= IT 1 
IF=IT2 
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972 

S7£ 

977 


96C 


5£2 

9es 

Sfi7 




IL 1=IT2 
11-1 = 1X4 

IF(i Tl.GE. IT2) IL=IT2 
IF( IT 1,GE . I T2) IH=IT1 
IF (IT3.GE . IT4 J 1L1=IT4 
IFC I T3.GE, I T4) IH1=IT3 
CC 977 I 1= IL, Ih 
CC 975 J=I4,IB 
1 \< = FLOAT( I I ) 

IF{Y1 ( I ) ,GT ,{ FLOAT! J* 1 )^DY ) ) T ( I , J ) = TB ( I , 1 ) 
IF((FLCAT(^ ) =!‘DY ) . GT .Y3 < 1 ) ) T ( I , J + 1 )=TS ( I , 3 ) 
IF( T( I,J ) .LT.TV.AND. T( 1, J+1) ,GE. TV) GO TO 970 
I F <T { I , J) . G£ ,TV . AND.T ( I, J-f 1 ) ,LT ,TV j GO TO 97C 
GC TO 975 

P>i={T{I,J)-TV)/(TCI , J )~T( 1,3+1)) 

CALL ISOTH!II ,I , J,DX,OY, FX ) 

V\RITE<e,97£) I,J,II,T(I,J) ,T(I ,J+1) ,RX 
FCFF AT (5X ,3 IS , 3F1 1 .4 ) 

CCMI NLE 
CCNT INUE 

C C 985 I 1= ILl , Ihl 
CO 9£2 J = IC ,ID 
TV=FLCAT (II) 

IF( Y4(I) .GT. (FLOAT! J-1 )=i'CY) ) T( I, J 3=TE( 1,4 ) 
IF( (FLGAT( J-1 ) =D Y) , GT. Y2( I ) ) T( I , J+ 1 ) = TS ( I ,2) 
IF (T ( I , J ) ,LT .TV .AND.T ( 1, J41 ) .C-E ,TV ) GO TO 980 
IF( T( I , J) . GE. TV. AND . T( I , 4+1 ) . LT, TV ) GO TO 980 
C-C TO 98 2 

FX={T(I,J)-TV)/(T(I,J)-T(I,J+1)) 

CALL ISOThC II , I , J ,DX,D Y, F)< ) 

V»FITE(6,97 2) I,J,II,T{I,J),T(J,J + 1),RX 

CCMI NUE 

CCNT INUE 

CCNTINUE 

FETLRN 

ENC 
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SiEFCLTINE ISQTH ( 1 1 , I , J , C X , DY • FX ) 

CIMENSION XISG<^C,30) ,YI£C(30i30) »TSCLT(30) ,NNY(30 ) 
INTcGER=!<4 TSCLT 

CC^^C^yI £Cy XI sc »YI SG » TSCLT -,^SCLT iNKY 


SLBKCLTIhE TC CETER^^INE X-CCCKO AND Y-CCCRD GF POINTS ON A REQUIRED I 
C 
C 

DC IGC NN=1,NSCLT 

IF (TSCLT (NN J- J I .EG .0 ) GC 7G 2CC 
C-C TC 100 

2CC NN2=NNY(NN) +1 

XISC(NN»NNZ) =FLDAT( I- 1 )*DX 
TI£C{NN,NNZ) = (F LCAT { J -1 ) +FX ) ^DY 
NNY( NiN)=NiVV(NN)+ 1 
GC TC 300 
ICG CONTINUE 
300 FETURN 
ENC 

ELCCK DATA 

CIMENSIGN X ISC( 3C,30) ,YI £C( 3C,20) , TSCLK 30 ,N.NV(30) 

INTEGEF^f^A TSCLT 

CCVMQN/ I SC/ XI SO ,YI SO ,TSCLT ,NSCLT,NNY 

DATA NNY,NSCLT,TSCLT/30^C, 2C,1SGO,1S50.1S60»2000 » 20 20 1,20 50,20 30 , 2 
1 ICC, 2150 ,2 200 , 2220, 2 23 5, 22 50,2 26 0 ,2280 ,2300,2320,2350, 2380,2400, 24 
220, 24A0,2AEC, 2480, 2500,2520, 2540 ,255 0 ,2560 ,2 570/ 

ENC 
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IKFOT CiTA 

tlACE eO^>^CARV FCIMS COCBCJNATtS 


CUTER eOUNCARV I^^ER BCIACAPY 


,hC. 

X 

T 

PT.NO. X 

V 

1 

0.0 

1.299959 

1 

0. lOOOCO 

1,312999 

2 

Q.0260Q0 

1.269995 

2 

0.116CQO 

1.295999 

3 

Q.GSOODC 

1.238000 

3 

o.isonuo 

1.287999 

A 

O.ICOOOO 

1.226999 

6 

O.2OC0C0 

1.268999 

c 

0.150000 

1.226000 

5 

0.250000 

1.261C00 

e 

0.200000 

1.216000 

6 

0.265000 

1.269999 

T 

0.250U00 

1.20SOCO 

7 

C.2CCCC0 

1.242999 

e 

0.266000 

1.200000 

0 

0.250000 

1.219000 

<i 

0.20COQO 

1.166995 

5 

0.3E10CC 

1.200000 

10 

0.350000 

1.166000 

IQ 

C.40QCC0 

1,106999 

11 

0.371000 

UlSCOCC 

ll 

0.650000 

1.150000 

12 

0.600000 

1.126CCC 

12 

0.698000 

1,099999 

13 

C.620CQ0 

1 .099959 

13 

0.5000CO 

1,096000 

1< 

0.650000 

1.077959 

14 

0.530000 

1.049999 

13 

0.6 750 JO 

1 .065995 

15 

O.5SOOC0 

1.030999 

u 

C.5C0000 

1.016999 

16 

0.572000 

l.CCOOOO 

17 

0.516000 

l.OCCCCO 

17 

C.6QOOCC 

0.961000 

le 

0.55OU00 

0.950000 

18 

0.607000 

0.950000 

1? 

0.5650CC 

C.9CC000 

19 

0.663Q0Q 

0.900000 

20 

0.60COOO 

0.Q780C0 

23 

0.65C0CC 

G.B90000 

21 

0.621000 

Q. 850000 

21 

0.678000 

C.8EC0CC 

22 

0.65COOC 

0.009000 

22 

0.700000 

o.nionoo 

23 

Q.657CBO 

o.eocooo 

23 

O.712CC0 

o'.eoocGO 

2^ 

0.693000 

0.75QOQ0 

24 

0.745000 

Q.75QQ00 

23 

O.7OC000 

0.728000 

25 

0.7500CC 

0.746000 

26 

0.726000 

0.700000 

26 

0.779000 

0.700000 

27 

Q.75C000 

Q.666GOO 

27 

a.OOOQCO 

0. 668000 

26 

0.763000 

0.650000 

28 

o.aiioco 

0.650000 

2S 

0.796000 

0.600000 

29 

O.8410G0 

0.6CCOQO 

20 

0.800000 

0.5960CO 

30 

0.850000 

0.580000 

31 

0. 833000 

O.SCOOCO 

31 

0.875000 

O.550C00 

32 

0.850000 

0.525000 

32 

0.900000 

0.512000 

22 

0.869CC0 

0.5C0CCC 

33 

D.9C6CO0 

O.5C0OOO 

36 

0.9OOUO0 

0.656000 

26 

0. 93^000 

0.450000 

2 j 

0. 903000 

0.650000 

35 

0.950000 

0.436000 

26 

0.960000 

0.600000 

36 

0-972CCC 

0.400000 

3i 

C.95Q0D0 

0.387000 

37 

1.000000 

0.256C0C 

26 

0.976000 

O.35CO0O 

3D 

1.002999 

0.350000 

36 

1 .000000 

0.316000 

35 

1.032000 

0.300000 

60 

1.013000 

0.300000 

60 

1.049959 

0.275000 

61 

1. 065999 

0.25COGO 

61 

1.Q630CC 

C. 250000 

62 

1.09UIOO 

O.2U0OOO 

42 

4.. 099999 

0.200000 

62 

1.099959 

0.191000 

63 

L.1300C0 

0.150000 

66 

1.136000 

0.150C0Q 

66 

1.1500CO' 

0. 13S000 

65 

1.150000 

0.135000 

6} 

l.ieioco 

O.IOCOOO 
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1.18lC.(iC 

q.ieoooo 


1.183000 

9.100000 


t.20C0OQ 

0.015000 

41 

1.1500GC 

0.150000 

4S 

l;224S55 

c;oso()oo 

48 

1.113995 

0.2CCOCO 

4S 

1.245595 

C.022CCC • 

45 

1.0SS55S 

0.22SOQO 

SO 

1.210000 

d. 050000 

50 

1.065555 

0.250000 

SI 

1.245555 

0.081000 

51 

1.05B599 

0.300000 

52 

1.240555 

C.ICOOOQ 

52 

I.C45555 

0.315000 

ca 

1.210595 

C. 150000 

53 

1.033955 

C.3SC00C 

54 

1.2C0000 

C. 168000 

54 

1.009999 

0.400000 

C £ 

1. laiOQC 

o. 2 cccoa 

55 

l.COOCCU 

0.423000 

5£ 

1.15C00Q 

0.250000 

56 

0.986000 

0.450000 

51 

1.122555 

C. 300000 

51 

0.961CCO 

0.500000 

58 

1.055999 

0.331000 

58 

0.5500C0 

O.518OC0 

S5 

I.C52595 

0.350000 

59 

0.933000 

0.550000 

60 

1.063000 

0.4CC0C0 

60 

O.SOSOCO 

0.600000 

6> 

1.045999 

0.4250CQ 

61 

0. 500000 

0. 605000 

62 

1.031595 

C.45COCO 

62 

0.0760C0 

0.650000 

63 

I.CC5CQ0 

0.5CCCCQ 

63 

0.85CQCU 

0,100000 

64 

l.COOOOQ 

a.5160C0 

64 

0.8220C0 

0.150000 

65 

C. 582000 

0.5S0CCC 

65 

C. 8QQOCO 

0.153000 

66 

0.956000 

0.600000 

66 

0.15£0Cfl 

0.8COOOO 

61 

C.5SOOOC 

0.605000 

61 

0.113000 

0.050000 

68 

0.929000 

0.65COQO 

68 

0.15C0CC 

O.B94000 

65 

0.9010Q0 

0.1GCQ00 

69 

0.146000 

0.500000 

10 

0.5CG0O0 

Q.7C40GO 

70 

0.124000 

0.950000 

11 

C.811C0C 

0.15CC00 

11 

C.UOCCO 

l.OCCOOO 

12 

o.esinoo 

o.eooooo 

72 

0.671000 

1.049595 

12 

G.85CCCC 

0.8C2C00 

13 

0.6S40CO 

1.059999 

14 

0.821000 

0.B500C0 

14 

0.65QCO0 

1.113000 

U 

C.EC3000 

0.500000 

15 

0.633000 

I.ISOOOQ 

16 

C.SCOOOQ 

0.5CEQC0 

76 

0.610000 

1.200000 

11 

o.iaiooo 

0.950000 

11 

0.600000 

1.215559 

18 

0.158000 

l.OGCOOO 

18 

0.585000 

1.249999 

15 

0.150000 

i.oiecoo 

IS 

0.56 lOQO 

1.255999 

£0 

C. 136000 

1 .049999 

80 

0.550000' 

1,3)6000 

81 

0.11300C 

1.055555 

61 

0.52EOC0 

1.349999 

82 

0.130000 

1.131000 

82 

0.5CQCC0 

1.384999 

83 

C. 691000 

1.150000 

63 

0.405000 

1.400000 

84 

C.61C000 

1.2C0000 

84 

0.450CCQ 

1.421999 

65 

C. 650000 

1.245599 

85 

0.413000 

1.4500DC 

£6 

•0.625000 

1.259595 

86 

0.4Q0OO0 

1.455000 

81 

C.6CCQ0C 

1.344555 

81 

0.250008 

1,469000 

£6 

C. 558000 

1,349955 

86 

0.300000 

1.465995 

£5 

0.564CCC- 

1.40C0OC 

05 

0.25CC00 

1.462999 

50 

0.55UOOO 

1.415000 

50 

C. 213000 

1.450000 

51 

C.S230C0 

1.450000 

91 

0.200000 

1.443999 

52 

0.500000 

1.465CC0 

92 

0.150000 

1.412999 

52 

0.453COO 

1.455559 

93 

0.138000 

1,400000 

54 

0.45C000 

1.500595 

94 

0.100000 

1.349999 

55 

0.4COOOO 

1.S21CC0 




56 

C.350UQQ 

1.532000 




51 

0.300000 

1.535555 




58 

0.250000 

1.526555 




55 

C. 200000 

1.5IIOC0 




ICC 

O.IBOCOO 

1.455555 




lot 

0.150000 

1.4e3QC0 




1C2 

0. 104000 

1.45CCCC 




1C2 

o.iooooo 

1.441CC0 




104 

C.C50000 

1 .4000C0 




105 

0.G1300Q 

1.345955 
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IKPtI CATA 


GAS TEPFECATUPE C 1STP IBUTICN' CN THE ELACE INNER AND CUTER SURFACE 


BCUNCAPV POINTS CN I-LINES FCR LChEP SURFACE CF OUTER BCUNCARV 
(START FRCP I» 2 TO I". 261 

»« 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2T60.00 2760.00 2760.00 2760.00 2760.00 
2760. CC 2760.00 2760.00 2760.00 2760. CO 276C.OO 2760.00 2760.00 2760.00 2760.00 


BCUNCARV POINTS CN I USES FOR UPPER SURFACE CF OUTER OOUNCARV 
(START FPCP I« 2 TC 1» 26) 

2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 

2760.00 2760,00 2760.00 2760. CO 2760.00 2760.00 2760. OC 2760.00 2760.00 2760.00 


eCUNCAPV POINTS CN 1 LINES FCR LONER SURFACE OF INNER BOUNOARV 
(START FRCP I» 3 TC I» 26 ) 

1560. CO 1560.00 1560.00 1560.00 1560.00 1560.00 1560.00 1560.00 1560,00 1560.00 1560.00 1560.00 1560.00 1540.00 1540.00 
1560. 6C 1560. CO 1560.00 1560. CO 1560.00 1560.00 1560.00 


I 

BCUNCARV POINTS CN I LINES FCR UPFER SURFACE CF INNER flOUNCARV 
(START FROM l» 3 TO I» 26 > 

1560. CO t56C.C0 1560.00 1560.00 1560.00 1540.00 1560.00 1560.00 1560.00 1560.00 1540.00 1540.00 1540.00 1540.00 1540.00 
1560. CC 1560. GO 1560.00 1560. CO 1560. CO 1560.00 1560.00 


eCUNCARV PCINTS CN J LINES FCR NEARER SURFACE CF CUTER EOUNCARY 
(START FRCP J» 2 TC 31 J 

2TAa.CC 2760.00 2760.00 2160. CO 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.09 2760.00 

2760.00 2760.00 2760.00 2760.00 2760.00 2160. CO 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 


BCUNCARV POINTS CN J LINES FOR FARTtiER SURFACE OF OUTER DOUNOARV 
(START FRCP J* 2 1C J» 31 I 

2760. CQ 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760*00 2760.00 
2160. CO 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 


BCUNCARV POINTS CN J LINES FCR NEARER SURFACE OP INNER BCUNCARV 
(START fPCP J« 3 TC J« 20 ) 

1560. OC 15^0. CO IS'iO.OO 1540.00 1540.00 1540.00 1540.00 1540.00 1540.00 1540.00 1540.00 1540.00 1540.00 1540.00 1540.00 
1540. OC 1540.00 1540.00 1540.00 1540.00 1340.00 1540.00 1540.00 1540.00 1340.00 1540.00 1340.00 1540.00 


ECUNCARV POINTS CN J LINES FDR FARTHER SURFACE CF INNER BCUNCARV 
(S1APT FRCP J« 3 TC J“ 30 ) 

1540.00 1540.00 1540.00 1540. QO 1340.00 1540.00 1540.00 1540.100 1540.00 1540.00 1540.00 1540.00 1540.00 1340.00 1540.00 
154G.CC I540.«a 1540.00 t540.<S0 IS40.00 1540.00 1340.00 1540.00 1540.00 1540.00 1540.00 1540.00-1540.00 


jffiEinvflb %)od idu 



tKPn DAT* 

FEAT CCEFTIUEMT CISTIM0UTION Q« THE BLADE INNEft AND CUTER 


eCLKCABV PCIMS CN I-LIRES FCP tCHEB SURFACE OF OUTER BOUNDARY 
(START FROM I« 2 TO I» 261 




4.51 2.40 2.12 

1,99 

1.90 

1.04 

1.70 

1.74 1.69 



2.3E 2,53 2. SI 

3.27 ■ 

3.SS 

4.31 

4.69 

4.94 5.17 



eCUNCAPY POINTS ON I 

LINES FCR 

UPPER 

SURFACE 

OF 

CUTER BOUHOARV 



(START FRCH I« 2 

TO I« 26) 




• 



4.51 4.06 4,58 

5.14 

5.63 

5.94 

6.11 

6.22 .6.31 



• J.TE 5.67 5.57 

5.47 

5.40 

S.33 

5.28 

S.24 5.21 



flCUNCAPV PCINTS CN t 

LINES FOR 

LQViER 

SURFACE 

OF 

INNER DOUNDARV 



(S1AFT FPCH I« 3 

TC I« 24 

) 






9.03 B.19 5.56 

5.42 

5.49 

5.56 

5.42 

5.36 9.21 



1.91 2.00 2.36 

. 2.57 

2. 78 

2,76 

2.78 




ECUNCAPY PCINTS CN I 

LINES FOR 

UPPER 

SURFACE 

OF 

INNER. BOUNDARY 



IS1APT FRCH 1* 3 

TC 1- 24 

) 




H 

to 


6.89 5.56 5.42 

5-49 

5,56 

5.42 

5.56 

5.42 5.56 

cn 


1.72 1.91 . 2.12 

2.49 

3.99 

4.62 

4,94 



BCChCARY PCINTS CN 4 LINES FCR NEARER SURFACE OF OUTER OOUNOARY 
(START FROM J» 2 TO J- JT I 

S.24 5.10 4.B6 4.61 4.31 3.99 3.63 3.06 2.81 

l.£E 1.T8 1.72 1.6t 1.67 1.67 1.71 1.77 1.88 


CCUNCAflV POINTS OH J LINES FOR FARTHER SURFACE OF OUTER QOUNOARV 
(START FPCH J« 2 TC 4« 31 I 

S.21 5.21 5.21 5.22 5.24 5.26 5.28 5.31 3.33 

5.66 5.72 5.78 5.82 3.88 S.S4 6.01 6.11 6.18 


8CUNCAPY PCINTS CN J LINES FCR NEARER SURFACE OF INNER BOUNDARY 
(START FROM J» 3 TC J» 30 1 

2.78 2.76 2.70 2.78 2.70 2.64 2.49 2.29 2.10 

1.39 1-39 uoa 2.78 3.89 5.21 5,56 ■ 5.49 5,56 


6CUNCART POINTS CN J LINES FOR FARTHER SURFACE CF INNER SCUNOARY 
(START FPCH J- 3 TO J« 30 t 

5.10 4.86 4.69 4.40 4.17 2.64 2.33 2.40 2.22 

1.39 1.39 1.39 1.E8 2.,43 -• 3.19 *3.96 S.i-7 3o42 


1.67 

5.20 

1.60 

i:?4 

1.04 

2.06 

2.26 

6.37 

5.21 

6.39 

6.35 

' 6.26 

6.09 

5.92 

3.47 

2.08 

1.39 

1,39 

1.53 

1.72 


3.49 3.21 3.33 1.80 1.39 1*51 


2.57 

6.2B 

2.67 

6.67 

2.37 

6.30 

2.26 

4.51 

2.13 

4.06 

2.04 

5.00 

5.38 

6.26 

5.64 

6,32 

5,47 

6.39 

5.51 

6.39 

5.56 

6.39 

5.60 

6.31 




1.90 1.08 1.79 1.61 1.31 1.39 

e.68 0.89 6.94 5.42 


1.61 


2.03 1.94 t.Ql 

3.56- - -9.42 - • 9.36 


1.70 

9.S6 


1.49 
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IhFUT CiT* 

CCCLANT MASS FLOW DIS1B!BU^IC^ 

8CUKCAPY PCIKTS ON I LINES FOR lOUER SURFACE OF INNER BOUNDARY 
ISTAPT FROM I« 3 TC !« 22 J 

0.0 A. 4100 6.3700 6.3700 8.3300 8.3300 10.2900 10.2900 12.2500 12.2500 26.4600 26.4600 26.4600 26.4600 26.4600 

26.460C 26.460C 26.4600 26.4600 26.4600 

BCUNCAPY PCUTS CH I LINES FCR UPPER SURFACE OF INNER BOUNCARY 
(ISTAPT FROM I« 3 TO !• 22 » 

0.0 4.4100 6.3700 6.3700 0.3300 6.3300 10.2900 10.2900 12.2300 14.2100 26.4600 26.4600 26.4600 26.4600 26.4600 

26.4600 26.4600 26.4600 26.4600 26.4600 


INPLT DATA 

VELCCITIES AT PCINTS WHERE ADIAOATIC HALL TEMP IS TO BE OETERHINEO 

etUNCAPV PCINTS CN I-LINES FCP ICWEP SURFACE OF OUTER BOUNDARY 
ISTAPT FRCP I« 22 TC I" 261 

1641.66 1635. CO 2030.00 2266.00 2416.00 

BCUNCAPY POINTS CN J LINES FOB REARER SURFACE Of OUTER BOUNDARY 
ISTAPT FBCP J< 2 1C J« a 1 

2350.00 2150.00 1975.00 1796.66 1641.66 1408.00 1216.33 
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IhPUT CATi 


PHYSICAL PARAHEIERS CF THE GRID SVSTEP 

DSOC.CSO CV-O.OSO ^XC« 105 hill- 94 

27 32 M" 2 N2= 25 

N3> 2 K4« 31 fkS« 21 H6o Q N7« 13 NPn 

PPCP6RTIES OF THE BLADE AND THE COOLAKT 

LB«0. 335000 XK><!.537SOO CP>>0,2T02aO 

PRACo 0.12 /LPHAa 0.0 TCOLOi° 1340.00 TENPC<« 1540.00 
D£NSC- 0.4513199925 CPC« 0.26B5999ST0 VISCCaO. 2690299080E-04. FLOC^ 
I1EFAT1CN CHECKS 

CPEo 1.000 SUHN» 0.020 hTEu lOO |HAH« 300 


0.2602 


DSIGINAIl PAGE IS 
OE POOR QUALITY 
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JUBBOLtlKE CLPVE -CUTfLT“ 

l«Q,O.TEP DCU^D«PY. l«ltI^^EP ECUNCflPV 
IP 6CLKCAPV FCIKT NLPCEP 

<E2rC21 IS THE POINT CN THE ECUNDAPV* 
Cf THE F/PECLIC CUPVE Y«#4E« > + C’»>« > 


(BlfCn MO IB3,C3I «RE THE SURROONOUG POINTS ON ElACE. a»B,Cf ARE 


I 


0 

0 

0 

1 

1 

0 

0 

I 

I 

0 

0 

1 

1 

0 

0 

I 

I 

0 

0 

I 

1 

0 

0 

L 

1 

0 

0 

1 

1 

0 

0 

1 

1 

0 

0 

I 

1 

0 

0 

t 

I 

0 

0 

1 

1 

0 

0 

1 

1 

0 


IP 


ICA 

A 

1 

103 

5 

3 

«2 

ICl 

6 

4 
SI 
SS 


Bt 


0.02400 
C. 10000 
C.OECCO 

a.ioojo 
C.13ECC 
0.10400 
0.10000 
0. 1 1400 
C. 20000 

c.ucco 

C. 15000 
0.15000 
C.213CC 
0.25000 


Cl 


1.25000 

1.44700 

1.23800 

1.350C0 

1.4COOO 

1.45000 

1.22700 

1.30000 

1.44400 

1.50000 

1.224CO 

1.28800 

1.45000 

1,52700 


P2 


0.05000 
0.05000 
C. 10000 
0.10000 
0. lOOOO 
C.IOOOO 
C. 15000 
C. 15000 
C. 15000 
0.15000 
C. 20000 
0.20000 
C. 20000 
C. 20000 


C2 


1.23000 

1.40000 

1.22700 

1.31400 

1.35000 

1.44700 

1.22500 

1.28000 

1.41300 

1.48300 

1.21000 

1.260CG 

1.44400 

1.51100 


7 

0.20000 

1.21900 

0.25000 

1.20500 

5 

C, 20000 

1.26900 

C. 25000 

1.26100 

E9 

0.300CO 

1.470CO 

0.25000 

1.46300 

58 

C.3C0CQ 

1.53600 

0,25000 

1.52700 

9 

C.264C0 

1.20000 

C. 30000 

1.1E7C0 

7 

0.28500 

1.Z500U 

0.30000 

1.2440C 

E8 

C.35QCO 

1.46900 

0.30000 

1. 47000 

97 

■ C. 35000 

1.53300 

0.30000 

1.53600 

to 

0.30000 

1. 16700 

0.35000 

1.16400 

a 

C. 30000 

1.2440C 

C. 35000 

1.21900 

t7 

C .40000 

1,4£SC0 

0.35C00 

1.46900 

56 

C.4C0C0 

1.52100 

0.35000 

1.53300 

12 

Q.3T1C0 

1.15000 

C. 40000 

1.13600 

10 

0.36100 

1.20000 

C. 40000 

1.1610C 

E6 

0.41300 

1.45000 

0.40000 

1.45500 

95 

C. 45000 

1.50100 

C. 40000 

1.521CC 

14 

0.42900 

1.10000 

0.45000 

1.07000 

11 

C. 40000 

1.18700 

C. 45000 

1.15000 

B4 

0.48500 

1.40000 

0.45000 

1.42600 

54 

C. 45300 

1. 50000 

0.45000 

1.5O1C0 

16 

0.475C0 

1.C5000 

C. 50000 

1.01900 

13 

y .45600 

i.inono 

0.50000 

1. 05800 

(2 

C.E26CC 

1.35000 

0.50000 

1.38500 

92 

0.52300 

1.45QOC 

C. 50000 

1. 465CG 

18 

0.51400 

l.UOUOO 

O.55OG0 

0.95000 

15 

C. 53600 

1.05000 

0.55000 

1.03100 

£0 

C.561CC 

1.3COQO 

C. 55000 

1.31600 

90 

0.564C0 

1.40000 

0.55000 

1.4190C 

20 

0.56SCC 

0.90000 

C. 60000 

0.0760C 

17 

0.57200 

l.CCOOO 

0.600CO 

0. 96100 

n 

C.61C00 

1.2C000 

0.60000 

1.22000 

67 

C. 62500 

1.30000 

0.60000 

1.34500 

22 

0.62100 

0.85000 

0.65000 

0.80900 

20 

0.64300 

0.90000 

C.65C00 

0.69000 

74 

0.654C0 

1.10000 

0.65000 

1.11300 

B5 

C. 67000 

1-20000 

0.65000 

1.25000 


83 

03 

A 

8 

0.10000 

1.22700 

1.26489 

-0.69665 

0.01300 

1.35COO 

1.32936 

1,64925 

C. 15000 

1.22600 

1.25900 

-0.52005 

0.11400 

1.30C0O 

100.00000 

loo.Donoo 

0.18000 

1.31400 

100.00000 

100.00000 

0.05000 

1.40000 

■ 1.33539 

1.46772 

0.20000 

1.21900 

1.21100 

0.28003 

0.20000 

1.26900 

1.32872 

-0.19011 

0.13800 

1.40000 

1.09580 

3,23474 

0.10400 

1.45000 

1.34445 

1.22121 

0.25000 

1.2D500 

1.20500 

0.34999 

0.25000 

1.26100 

1.41100 

-1.15004 

0.15000 

1.41300 

1.24458 

1.49955 

0.18000 

1.50000 

1.28272 

1.79873 

0.26400 

1.20000 

1.21475 

0.26279 

0.28500 

1.25000 

1.21022 

■ 0.65697 

0.21300 

1.45000 

1.24581 

1,47579 

0.20000 

1,51100 

1. 37700 

0.95017 

0.35000 

1.16400 

1.20429 

0.28724 

0.35000 

1.21900 

1.23242 

0.49882 

0.25000 

1*46300 

1.30799 

1.01994 

0.25000 

1.52700 

1.30201 

1.49907 

0.37100 

1.15000 

1.01939 

1.43274 

0.30100 

1.20000 

1.24T6I 

0.40591 

0.30000 

1.47000 

1,20299 

1,67014 

C. 30000 

1.53600 

1.36496 

1,11009 

0.42900 

I.IOOOO 

-0.61070 

9,59536 

0.45000 

1. 15000 

1.33760 

-0.05368 

0.35000 

1.46900 

1.33443 

0.96731 

0. 35000 

1.53300 

1.39298 

0.96017 

0.47500 

1.05000 

1.24523 

0.33498 

0.49000 

1.10000 

0.92805 

1.87601 

0.41300 

1.45000 

1.16536 

1.86765 

0.40000 

1.52100 

1.90804 

-1.47368 

0.5140Q 

1.00000 

0.92663 

1.67391 

0.53600 

1.05000 

-0.04059 

5.56863 

0.48500 

1.40000 

0.47525 

4. 72607 

0,45300 

1.50000 

1.25980 

1.60035 

O.SBSOO 

0.90000 

1.55580 

-0.79414 

0.57200 

1.00000 

3.41223 

-7.14569 

0.52800 

1.35000 

2.97744 

-4.54135 

0.52300 

1.45000 

0.50309 

4,31923 

0.62100 

0,05000 

3.04908 

-5.83246 

0.60700 

0.95000 

0.04511 

4.59130 

0.58500 

1.25000 

2.40039 

-1.98438 

0.59800 

1.35000 

12.46875 

-34.54544 

0.65700 

0.80000 

3.17572 

-5,96721 

0.67800 

0.65000 

1.31919 

-1.41739 

0.63300 

1.15000 

-18.47290 

63.29166 

0.62500 

1.30000 

-1.96834 

12. 1803T 


3.17883 
“4.72301 
1.00997 
100.00000 
100.00009 
“3.51010 
“1.19994 
“0.54251 
-7.47174 
-1.98394 
-1.39997 
2.20009 
“2.51620 
-3.28668 
-1.20549 
“1.81558 
-2.42940 
-1. 40001 
-1.14973 
“1.54046 
-1.59976 
-2,39964 
-2.9U05 
-1.39434 
-2,60023 
-1.80008 
-13.07059 
-0,80803 
-1.66220 
-1.60012 
-1,66895 
-3.07883 
-2.05358 
1.26316 
-2.96696 
-6.58823 
-5.01518 
-2^37983 
-0.55896 
5.11920 
2.76692 
-5.09615 
3.68586 
-5.10435 
-0.01563 
26.72726 
3.57377 
-0.00870 
-51.04166 
-11-12732 


XHN« -0. 
XHN* 1. 
XHN« -0. 
XMN=»119. 
XHN»119, 
XMN« 0 
XHN* -0 
XHN* -0 
X«N- 0 
X“N=> 0 

XHN=. -0 
XHN- -0 
XKN» 0 
XMN« 0 
XHN’ -0 
XHN» -0 
XMNa 0 
X^N^ 0 
XHN= -0 
XHN- -0 
XHN» 0 
XMN=> 0 
XHH= -0 
XKNa -0 
XHN« -0 
-0 

XHH* -0 
XKN« -0 
XHN= -0 
XKN= -0 
XHN=- -I 
XHN« -0 
XHN« -0 
XHN» -0 
XMN=< -1 
XHM» -I 
XHH= -I 
XHN- -0 
XHN» -1 
XHN» -1 
XMNa -1 
X^N" -1 
XHH» -1 
XHN« -I 
XHN« -2 
X«N» -3 
XHN<« -1 
XHN.: -1 
XHN- -2 
KHK« -3 


THE COEFS 


3TB962 
176443 
120056 
999985 
999985 
765696 
079949 
352866 
.993222 
,626027 
.209995 
.270001 
.493033 
.404061 
.339958 
.250813 
.261087 
.250165 
.402599 
.425460 
.O6OO80 
.060081 
.604996 
.570126 
.150015 
.149967 
.661109 
,700739 
,362450 
.319927 
.077091 
.894144 
.700575 
.336842 
.313042 
.019607 
.089108 
.771473 
.408993 
.514571 
.497745 
.28653T 
.409425 
.533912 
1,003124 
>.472733 
1.321313 
t. 428695 
1.062485 
'.205142 
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^ER4•tIC^» 
n£P/STICK« 
tt£p«uc^» 
lTEP.^^lC^ = 
iiEPancNa 
nEPAHCN» 
ITEPAT ICN« 
nePAiKP* 
I1EPATICN* 
HEFattCh:' 
iiEpancp- 
tTEPa T ICSa 
ITFPanCPa 
^£PA1Ic^» 
JTEP4T IC^« 
HEf aiICF» 
^EPATIC^“ 
IlcPAT ICP« 
nE'ilICF« 
HEPAl ICA' 
ITEP41 tC'* 
I7EPA11CN* 
tTEPAT ICA:> 
I7EFA7 ICA* 
UFRAIICA’ 
17 £fA 7 ICA= 
I7EPATICA* 
ITEPA7ICA* 
r7E3A7!CA= 
nERA7ICK« 
nFPA7iCK« 
llEPAnCA* 
nEBA1ICA« 
!TEPA7ICA« 
17EPA71CA* 
nEPA71CA = 
IIEPATICA* 
I7ERA7ICK- 
ITEPAUCA» 
HEPATICA* 
IIEPAIICA* 
17EFA71CA" 
nFFATICA* 
I 7EPA7ICA* 
7TEAA7ICA» 
tlEFAIlCA* 
I7EPA7 ICA* 
17EPA7ICA- 
I7EPA7ICA* 
17EP47ICA* 
I7EPA7 ICA* 
nEPAIlCA* 
ITEBA7ICA* 
IIEPAIICA* 
tXEPA71CA» 
I7EPA71CA* 
tlPPA7ICA» 
IIEPA7ICA* 
nEPATICAt. 
tTEFAllCA- 


1 

ERPCB» 

2 

EPPCR= 

3 

EPRCP* 

4 

ERRCRa 

5 

EPflCR» 

4 

BPRCRx 

7 

FRBCRx 

B 

EPRCR* 

9 

FRRrp, 

10 

FRRCRx 

11 

EPnCPx 

12 

EPRCR* 

13 

EBRCRx 

14 

EBRCPx 

IS 

ERPCRx 

16 

EPRCRx 

17 

EPPCPa 

18 

EPRCR= 

19 

EPRCP- 

20 

EPRCRx 

21 

ERRCRx 

22 

EPRCRx 

23 

ERRCRx 

24 

EflPCpx 

2S 

EPRCPx 

26 

ERRCRx 

27 

EPPCP= 

28 

EBRCPx 

29 

EPRORx 

30 

EBPCP* 

31 

EPRCR= 

32 

EPfiCRx 

33 

EPRCPx 

34 

ERPCRx 

3S 

EPPCPx 

36 

ERROR" 

37 

ERRORS 

38 

ERROR* 

39 

ERROR* 

40 

ERROR* 

Al 

ERROR* 

42 

FRRCR* 

43 

ERROR* 

44 

ERROR* 

45 

ERROR* 

46 

ERROR* 

47 

ERROR* 

48 

ERROR* 

49 

ERROR* 

SO 

ERROR* 

51 

ERROR* 

S2 

ERROR* 

53 

ERROR* 

54 

ERROR* 

55 

ERROR* 

56 

ERROR* 

57 

ERROR* 

56 

ERROR- 

59 

ERROR* 

60 

ERROR- 


O.IEEEAE 07 
0.69ECAE 04 
0.23E59E 04 
0.12mE 04 
0.7S727E OS 
O.S'iSSAE OS 
O.AU24E OS 
0.222;2E OS 
0.2SA37E OS 
0.20USE OS 
0.1SS92E OS 
O.liXSlE 05 
O.IOISSE OS 
O.E17A1E 04 
0.4S44SE 04 
0.52e54E 04 
0.42414E 04 
C.244CCE 04 
0.27ei9E 04 
0.22517E 04 
O.Ia242E 04 
0.14TU7E 04 
0.11S91E 04 
0.S7247E 03 
0.78fl5eE 03 
0.62S’5E 03 
O.S160BE 03 
0.41944E 03 
0.33S24E 03 
0.27432E 03 
0.221S1E 03 
0.17E74E 03 
0.144C4E 03 
0.US4SE 03 
O.S3221E 02 
G.74E41E 02 
O.4071SE 02 
0.4S216E 02 
0.3E44IE 02 
0.30S44E 02 
0.247t2F 02 
O.ISECIE 02 
0.15827E 02 
0.12447E 02 
0. tOt07E 02 
0.e048SE 01 
0.44427E 01 
O.S1499E 01 
0.41U3E 01 
0.3277IE 01 
0.26174E 01 
0.2C928E 01 
O.J6717E 01 
0. 1334SE 01 
0.10456E 01 
0.8SC39E 00 
O.4ei03£ 00 
0.S4222E 00 
0.43450E 00 
0.34733E 00 



T-NATBIX, GIVES Alt IMEPICB PClATS.IMTIil SETTING OF T tS 1600.0 
CMV THOSE INSIDE ELADE APE CHANGED AND PILL OE DIFFERENT FROM 1600.00 

J » 1 2 3 6 5 6 T a 9 10 U 12 13 14 15 

t 

16 17 18 IT 20 21 22 23 24 2S 

1 2400. OC 24C0.D0 2400.00 24CO.OO 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 

Z4CO.OO 2400.00 2400.00 2400.00 24S0.C0 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2333.60 2400.00 2400.00 2400.00 

24C0.CC 24CC.OO 

2 2400.00 2400.00 2400.00 2400.00 24CQ.OO 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 

24CC.CC 24C0.0C 24CC.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 22SS.86 2181.92 2145.64 20Q9.64 2400.00 

24C0.CC Z4CC.00 

3 24C0.OG 24C0.00 2400.00 2400.00 2400.00 24CO.OO 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 24Q0.00 2400.00 2400.00 

24CC.0C 2400.00 2400. OC 2400.00 2400.00 24CC.00 2400.00 2400.00 2400.00 2400.00 2101. T4 1992.52 1875.01 1996.63 2400.00 

24CO.CO 24CO.00 

4 2400.00 24C0.0C 2400.00 2400. CC 2400.00 2400.00 2400.00 2400.00 2400.00 2400,00 2400.00 2400.00 2400.00 ^400.00 2400.00 

24CO.OO 24CO.OO 2400.00 2400.00 2400.00 24C0.C0 2400.00 2400.00 2400-00 2400.00 2002-B6 2400.00 2400.00 2400.00 2095.87 

24CC.CC 24CC.C0 

5 2400. OC 24C0.00 2400.00 24CO.OO 24C0.C0 24CC.Q0 24C0.Q0 24D0.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 

2400. OC 24C0.0Q 2400.00 2400.00 2400.00 2400.00 24C0.00 2400,00 2400.00 2400.00 1958-21 2400.00 2400.00 2400.00 2087.57 

H 2235.26 24(0. CC 

LO 

H 6 2400. OO 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400. 00 2400.00 2409.00 2400.00 2400.00 2400.00 2400.00 2400.00 

24CC.CC 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 1931.79 2400.00 2400.00 2400.00 2400.00 
<22S,E6 24CC.00 

7 2400.00 24Ca.0C 2400.00 2400.00 2400.00 24CO.OO 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 

24CC.CC 24C0.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 1953.31 2400.00 2400.00 2400.00 2400.00 2400.00 

2221.77 24CC.O0 

a 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2408.00 2400.00 2400.00 2409.00 2400.00 2400.00 2400.00 2400.00 2400.00 

24C0.CC 2400.00 2400.00 2400.00 2400.00 24C0.00 2400.00 2400.00 2400.00 1912.25 2400.00 2400.00 24QQ.OO 2400.00 2400.00 

1253.41 24CC.CC 

9 2400.00 24C0,00 2400.00 2400.00 24C0.0a 2400,00 2400.00 2400. CO 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 

24C0.OQ 24C0.oa 2400.00 2400.00 2400.00 2400.00 24CO.OC 2400.00 1906.22 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 

2290.63 24CC.0C 

to 2400.00 2400.00 2400.00 2400.00 2400.00 2400. 00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2409.00 2400.00 2400.00 

24CC.C0 24CO.OO 2400.00 2400.00 2400.00 2400.00 2400.00 1903. S4 1BS3.21 2400.00 2400.00 2400.00 2400.00 2400.00 2226.25 

2342.61 24CC.CC 

11 2A0Q.0C 24GO.OO 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 

24CO.OC 24CO.OO 2400.00 2400. 00 2400.00 24CC.OO L93U61 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2196.80 2310.24 

240a-CC 2400.00 

12 2400.00 2400.00 2400.00 2400.00 24Q0.00 2400,00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 

2400.00 2400.00 2400.00 2400.00 206C.61 2006.41 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2199.46 2308.93 2400.00 

24C0.CC 24CC.OO 

13 2400.00 24C0.no 2400.00 2400.00 24C0.C0 24CC.C0 2400.00 2400.00 2400.00 2400. OQ 2400.00 2400.00 2400.00 2400.00 2400.00 

2400. OC 24C0.00 2400.00 2109.41 2090.76 2400.00 24C0.Q0 2400.00 2400.00 2400. 00 2158.02 2240.16 2408.00 2400.00 2400.00 

24CC.CC 24CC.C0 

14 24CQ.0G 24C0.OQ 2400.00 2400.00 2400. CO 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 2400.00 

24CC.00 24G0.O0 2180.76 2400.00 2400.00 2400.00 2400.00 2400.00 2296.12 2296.96 2296,84 2400.00 2400.00 2400.00 2400.00 

24C0.CC 24CC.C0 


ORIOINAr; PAGE li 
m POOR QUAIOT 
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G^S tehpepature oistributich cn the claoe inner and cuter surface 


eCUKCARY POIKTS CN l-LIRES FCft LCHER SURFACE OF OUTER BOUNCftY 
(START FROR I» 2 TC 1= 26 ) 

2760. GO 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 
2760. CO 2760.00 2760.00 276C.00 2760.00 1765.57 1765.57 1765.57 1865.15 1027.63 


OCUACAPY FCUTS CN I LINES FCB UPPER SURFACE OF CUTER BCUNCRY 
(START FRCP 1= 2 TO 1= 26 1 

2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 
2760. CO 276C.C0 2760.00 2760.00 276C.C0 2760.00 2760.00 2760.00 2760.00 2760.00 


BCL'NCARY POINTS CN ! LINES FOR LONER SURFACE OF INNER BOUNDARY 
(START FRCP 1= 2 TC I» 2R J 

1560.00 16C1.52 1608.67 1626.05 1617.66 1626.92 1623.05 1633.28 1626.99 1635.81 1671-22 1668.05 1702.36 1712.81 1721-07 
1731.61 1761.95 1753,36 1765.55 1560.00 1S6C.OO 1560.00 


OCUNCARY PC(NTS CN I LINES FOR UPPER SURFACE OF INNER BOUNOARY 
(START FRCP I* 3 TC 1= 26 1 

1560. CO 16C5.92 1609.86 1636.13 1633.60 1656. 66 1650.36 1669.06 1667.67 1669.75 1671.22 1683.05 1702.36 1712.61 1721.87 
1731.61 1761-95 1753.36 1765.55 1165.57 1765.51 1765.57 


DCUNCAPY PCINTS CN J LINES FOP NEARER SURFACE OF CUTER BOUNDARY 
(START FRCP J= 2 TC 8* 31 I 

1900.16 1B27.17 1765.57 1765.57 IT6S.57 1765.57 1765.57 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 
2760. CO 2760.00 2760.00 27AC.00 2160.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 


BOUNDARY POINTS CN J LINES FOP FARTHER SURFACE OF OUTER BOUNDARY 
(START FRCP J» 2 TC J= 31 1 

2160. CO 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2750.00 2760.00 2760.00 2760.00 2760.00 2760.00 

2760.00 2760.00 2760.00 2760.00 2160.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 2760.00 


2 S 

HTjg 

2 

& > 
:1a 


DCUNCAPY PCINTS CN J LINES FOR NEARER SURFACE OF INNER BOUNDARY 
(START PROP 3 TC J= 30 1 

1560.00 1560.00 1560.00 1560.00 1765. SS 1765.55 1753.36 1761.’95 1761,95 1731.61 1721.87 1721.87 1712.81 1702.36 1702.36 
1608.05 1671.22 1671.22 1625.01 1626.99 1633.28 1623.05 1629.92 1626.05 1560.00 1560.00 1560.00 1608.67 


BCUNCARY PCINTS CN J LINES FOP FARTFER SURFACE CF INNER BOUNOARY 
(START PROP J- 3 TC J= 30 > 

1765.57 1765.57 1765.57 1765.57 1765.57 1765.57 1765.57 1765.55 1765.55 1753.36 1753,36 1761.95 1761.95 1731.61 1721.87 

1721.01 1712.81 1712.81 1712.81 1702.36 1702.36 1680.05 1608.05 1671.22 1671.22 1669.75 1667.67 1669.86 
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TEf'PEPATUPE CISTRIBUTION CN THE 0CUNCARY 


eCUNCARY POINTS ON I-LINES FOR LCHEP SURFACE OF 
(START FPCH I>= 2 TC I« 26 I 

** 2271.97 2133. 6A 2060.33 20C5.50 1999.09 1970.72 l95fi.S5 
2295. SS 2313.64 2392.89 23*9.85 2760-01 2400-00 2400.00 


eCUNCAPY POINTS CN I LINES FOR UPPER SURFACE OF 
(START FPCH I* 2 TC I» 26 ) 

*+ 2099.64 2099.64 2199-12 2267.69 2306.02 2330.80 2342-52 
2548.49 2562.91 2534.69 2488.94 2446.68 2379.60 2318.13 


BCUNCARV POINTS CN I LINES FCR tCUER SURFACE GF 
(START FPCH 1= 3 TO 1= 24 1 

1964.56 1917.97 1925.76 1914.34 1092.72 1087.71 1071.10 
22E9.;96 2320.23 2327.36 2662.64 24CC.C0 24CC.C0 2400,00 


OCUNCARV POINTS ON I LINES FCR LPFEP SURFACE GF 
(STAFT FRCH I« 3 TC I« 24 1 

1875.01 1982.40 2070.26 2122.10 2141.22 2164.90 2162.32 
2460.25 2413.16 2365.13 2319.93 2221.64 2222.48 22CB.61 


BCUNCARY PCINTS CN J LINES FOR NEARER SLRFACE OF 
(STAFT FPCH J= 2 TC 0* 31 » 

2281.42 2400.00 2400.00 2400.00 2400.00 2400.00 2564.93 
2236.10 2182.47 2104.04 2104.04 2019.87 1949.46 1924.02 


BCUNCARY POINTS CN J LINES FCR FARTHER SURFACE OF 
(START FRCH J= 2 TC 31 ) 

2367.41 2354.66 2369.99 2345.94 2350.45 2360.00 2349.90 
2571-14 2545.72 2577.34 2539.35 2564.75 2527.55 2503.89 


CUTER BCUNCRY 

1920.24 1927-aO 1971.12 2068.61 2139.47 2225.48 2246.73 2295.20 
2400.00 2246.49 2287.32 

OUfER BOUNORY 

2346.23 2345.40 2352.12 2342.27 2278-67 2296.84 ‘2354.52 2523.31 
2332.20 2341.65 2376.58 

INNER eOUNOARY 

1853.27 1874.44 1962.43 2077.11 2148.73 2210.82 2228.76 2256.72 
INNER aOUNCARY 

2173.06 2167.65 2119.41 2124.46 2261.97 2436.80 2474.10 2485.95 
OUTER BQUNCARV 

2381.37 2381.37 2353.27 2304.79 2295.57 2295.57 2270.96 2236.10 

1935.40 1984.91 2340.34 2333.60 2352.32 2099.64 2099.58 2237.61 

OUTER BOUNDARY 

2408.46 2412.40 2457.69 2447.26 2496.53 2492.20 2530.22 2536.21 

2393.78 2386.13 2302.56 2322.78 230i;92 2345.98 2350.44 2345.13 


BCUNCARY POINTS ON J LINES FOR NEARER SLRFACE OF INNER DCUNOARY 
(START FPCP J* 3 TC J» 30 J 

2400.00 24QO.OO 2400.00 2400.00 2400.00 2760.01 2333.97 2334,32 234 1. 52 2271.05 2249-23 2265.81 2245.57 2202.35 2215.18 
2173.39 2090.37 2094.58 2000-93 1902.20 1854.12 1844.94 1887.03 1918.06 1936.84 1875.01 1925.65 2084.42 


BCUNCAPY PCINTS CN J LINES FOP FARTHER SURFACE OF INNER BOUNDARY 
(START FRCH J= 3 TC J* 30 1 

2206.44 2210.31 2212.65 2195.57 2226.05 2299.52 2325. 7T 2358.35 2369.68 2408.16 2416.14 2458.68 2462.01 2491.47 2462.07 
2498.15 2448.19 2414.57 2437.95 2341.32 2273.00 2187.54 2142.20 2100.27 2131.96 2155.52 2150.29 2166.33 
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TEMPERATURE CISTMBUTION INSIDE THE BLADE AND ON ITS SURFACE 
Ja 

31 2237.61 2235.28 2229.68 2231.77 2253.91 2290.63 2342*81 2349*19 

30 • 2C99.se 2C95.B7 2087.57 2064.42 

««««« 2166.33 2226.25 2310.24 2350.44 

29 2099.64 1996.63 1926.65 

2156.29 2196. CO 2308.93 2349.90 

26 2352.32 2145.69 1675.01 

thVItt 2155.52 2199.46 2301.92 

27 2233.60 2181.92 1992.52 1936.84 ‘ 

thhlitl 2131.96 2240.16 2322.76 

26 2340.34 2255.86 2101.74 2002*86 1958.21 1931.19 1916*08 

*«««« 2108.27 2158.02 2302.96 

25 19E4.91 1953.31 1912.25 1887.03 

UIIIIHH 2 142.20 2298.96 2388.13 

24 1939.40 1906.22 1644.94 

2187.54 2256.12 2393.78 
23 1924.02 1903.84 1854.12 

2273.00 2432.48 2303.89 
22 1949.46 1931.67 1902.20 

*»««« 2341.32 2419.07 2527. 5S 

21 2Ct9.87 2006.41 2000.93 

#««/)« 243T.95 2538.22 2S64.7S 

20 2104.04 2C90.76 2C94.58 

uitvtut 2474.57 2502.74 2539.39 

19 2104.04 2109.41 2090.37 

2448.19 2456.76 2563.36 2377.34 
18 2162.47 2180.76 2173.39 

2490.15 2320.93 2349.72 
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tSOTI*ERH>L 1. INE LOCATIONS 
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0 .7«ea 
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2 

28 
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3 

26 
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2145*6931 

2099.6396 

0.7J16 


2 

26 

21 18 
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2055*6396 

C. 7C59 


2 
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21 14 

2145*6921 

2099.6256 

C.6E02 


i 

26 

2 1 IS 
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2099*0396 
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0.2000 
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0.1000 
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CCORC 
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1.2611’ 

1,2626 

1.2265 

1.2027 

1.1552 

1.0266 
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EFF1CIE^CV OF CODLING 


ECU^CAPV POINTS ON l-LINES FCR tCUEP 
(START FRCP I» 2 TC !* tk ) 

SURFACE 

OF 

CUTER BCUNCRV 







O.AOQO 0.S134 0.5899 0.6189 0.6237 
0.3804 0.36S9 0.3U09 0.3198 -C.OOOO 

0.6469 

0.2951 

0.6569 

0.2991 

0.6683 

0.2951 

0.6021 

0.4209 

0.6466 

0.3874 

0.5667 

0.9086 

0.4381 

0.4207 

0.3810 

eCUNCAPV POINTS CN I LINES FOB UPPER 
(STAFT FRCP 1= 2 TC I» 26 1 

SLPFACE 

OF 

CUTER BCUNCRV 







0.3218 0.5413 0.4597 0.4035 0.3721 
0.1734 0.1615 0.1845 0.2222 0.2568 

0.3518 

0.3116 

0.3422 

0.3622 

0.3392 

0.3507 

0.3398 

0.3425 

0.3343 

0.3143 

0.3424 

0.3945 

0.3796 

0.3324 

0.1940 

ECUNCAPV POINTS CN J LINES FCip NEARER SUPFACE OF 
(SIAPT FfiCP 2 TC • 31 1 

OUTER BOUNUARV 







0.3923 0.2551 0.2951 0.2951 0.2951 
0.4294 0.4734 0.5377 0.5377 0.6067 

0.2951 
0.66 44 

0.1599 

0.6852 

0.3103 

0.6726 

0.3103 

0.6353 

0.3334 

0.3440 

0.3731 

0.3495 

0.3007 

0.3342 

0.360T 

0.5413 

0.4009 

0.5413 

0.4294 

0.4282 

OCLNCARV POINTS ON J LINES FOP FARTHER SUPFACE OF 
(START FRCP 2 TC J= 31 ) 

OUTER BOUNDARY 



. • 




0.3218 0.3322 0.3197 0.3394 0.3357 
0.1548 0.1756 0.1497 0.1809 0.16C0 

0.3279 

0.1905 

0.3361 

0.2099 

0.2881 

0.3002 

0.2049 

0.3048 

0.2478 

0.3749 

0.2563 

0.3584 

0.2160 

0.37S5 

0.2195 

0.3394 

0.1818 

0.33S7 

0.1034 

0.3401 
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APPENDIX B 


HEAT BALANCE EQUATIONS TO OBTAIN COOLANT TEMPERATURES 


In this appendix, the equations governing the coolant tem- 
perature distribution will be derived. A one-dimensional flow 
analysis is performed down each inter tube-vane passage, account 
being taken of the area changes, mass addition due to additional 
impingement from the flanks of the tube and heat addition from 
the vane walls. The temperature of the coolant inside the 
impingement tube is assumed to be constant. 

Referring to Figure 16, I=N1 and I=N7 are the I mesh lines 
enclosing the impingement tube without intersecting it and 
N71 = N7-1. To derive equations governing the temperature 
distributions, consider the control volumes in the coolant flow 
path. A typical control volume is made up of two consecutive 
I mesh lines and two surfaces . The four types of control 
volumes (CVl, CV2, CVS and CV4) which are considered to cal- 
culate the coolant temperature rise are shown in Figure 17. 

The sketches illustrating the m.ass balance and heat balance 
for control volumes are shown in Figures 18 and 19 respectively. 

Consider first the control volume CVl, shown in Figure 19. 
The heat balance equation for this system can be expressed 
semantically as : 

Influx of 
enthalpy at + 

I mesh line 


or mathematically as 


Influx of Rate of heat 

enthalpy from + inflow from 
the impingement the vane wall 

tube 


Efflux of enthalpy 
at Cl+1) mesh 
line 
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• • 

m(i,3)C TGCi,3) + Am{i,3)C TCOLN 

*r 

+ h(i,3) [TB(i,3) - TG(i,3)'] As(i,3) LB 


where mCi,3) 

' TG(i,3) 

Am (i, 3) 

TCOLN 
h(i,3) 
TB{i,3) 
As {i ,3) 

LB 

3) 

TG(i+l,3) 


m(i+l,3)C TG{i+l,3) (b.l) 

IT 

- total mass of coolant entering the control 
volume at l mesh line per unit time (see 
Figure 18) . 

= specific heat of coolant at constant pressure. 

= coolant temperature at I mesh line intersection 
point on the inner lower surface. 

= coolant mass added to the control volume between , 
I and I+l mesh lines from the impingement tube, 
per unit time. 

= temperature of the coolant in the impingement 
tube . 

= heat transfer coefficient at I mesh line inter- 
section point on the blade inner lower surface. 

= blade temperature at I mesh line intersection 
point on the blade inner lower surface. 

= surface distance between consecutive points on 
inner lower surface intersected by I mesh lines, 
(see Figure 18) . 

= length of the blade. 

= total mass of coolant leaving the control voliame 
at (I+l) mesh line per unit time. 

= coolant temperature at (I+l) mesh line inter- 
section point on the inner lower surface. 
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The equation (b.l) can be written in a more convenient form 
as follows : 


TG{i+l,3) = 


m(i,3)C TG(i,3) + Am(i,3)C TCOLN + Q (i , 3) As (i , 3) LB 

E — E 

m{i+l,3)C 


(b.2) 


where Q(i., 3) = h(i/3) [TB(i,3) - TG{i,3)] (b.3) 

Following the same procedure described above, a similar 
equation can be obtained for the control volimie CV2 shown in 
Figure 19 as follows: 


m{i,4)C TG{i,4) + Am(i,4)C TCOLN + Q {i , 4) As (i , 4) LB 

TG(i+l,4) = 2 E 

m(i+l,4)C„ 


•(b.4) 


where the second subscript * four * for the variables appearing 
in equation (b- 4) refers to the conditions at a point on the 
upper inner surface. 

Referring to Figure 19, the blade inner upper and lower 
surfaces and the mesh lines I=N7 and I=N71 define the boundaries 
for the control volume CVS. The two streams of coolant coming 
from the lower and upper surfaces of the impingement tube merge 
at I=N71 line. The heat balance equation for this system 
(CVS) can be written as follows: 


• « 

m(i,3)C TG(i,3) + m(i,4)C TG(i,4) + Q(i,3) As(i,3)LB 

P P . 

+ Q(i,4) As(i,4)LB 

= mCitl,3)C TG(i+l,3) 

P 

= mCi+l,4)C TG(i+l,4) (b.5) 

hr 
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Since/ for I ^ N71, mCi+1,3) = mCi+1/4)/ equation (b.5) can 
be rewritten as follov/s : 


TG(i+l,3 = TG(i+l,4) 


[m(i/3) C TG(i,3)+m(i/4)C TG(i,4).+Q(i/3) as (i , 3) LB+Q (i , 4 ) As (i/4)LB] 

■ . _P P . ■ . - 


m(i+l/ 3) C 


(b.6) 


Finally, the heat balance equation for the system CV4 shown 
in Figure 19 can be written as follows; 


« 

m(i,3)C TG(i,3) + Q{i,3) As(i,3)LB + Q(i,4) As(i,4)LB 

tr 

= m(i+l,3)C TG(i+l,3) (b.7) 

P 

• • 

Since, for I > N71, m(i+l,3) = m(i+l,4), equation (b.7) can 
be rewritten as follows : 


TG(i+l,3) = TG{i+l,4) = 


m(i,3)C TG(i,3) + Q(i,3) As(i,3)LB + Q(i,4) As(i,4)LB 

m(i+l,3)C 


(b.8) 


The equations (b.2), (b.4), (b.6) and (b.8) are incorporated 
in the main program to calculate the coolant temperature rise 
over an incremental length As. 
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